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XIX. BAND ERSTES HEFT 1951 


Der zeitliche Verlauf des Aufspulvorganges einer instabilen 
Unstetigkeitsflache von endlicher Breite. 


Von W. Kaufmann. 


1. Einfithrung. In einer bemerkenswerten Abhandlung hat H. Kaden den Aufspulvorgang 
einer ebenen, instabilen Unstetigkeitsflache untersucht, die sich einseitig ins Unendliche erstreckt!. 
Die wesentliche Kigenschaft einer solchen ist, daB der Vorgang infolge ahnlicher VergréBerung 
nach einem vorgegebenen Zeitgesetz erfolgt, dergestalt, daB sich die Zeiten t und t, zu den Mab- 
staben r und r, der ahnlichen Strémungsbilder verhalten wie 


alee (1 


Dieses Gesetz gilt jedoch nicht mehr, wenn es sich um Unstetigkeitsflachen von endlicher 
Breite handelt, wie sie insbesondere hinter Auf- oder Quertrieb erzeugenden Fliigeln oder Schau- 
feln bei Flugzeugen und Strémungsmaschinen auftreten. Es soll deshalb hier der Versuch unter- 
nommen werden, etwas iiber den zeitlichen Verlauf der 
Aufwickelung einen ebenen Unstetigkeitsflache von 
endlicher Breite und vorgegebener, zur Flachenmitte 
symmetrischer Zirkulationsyerteilung auszusagen. Die 
strémende Fliissigkeit wird dabei durchweg als’ un- 
zusammendriickbar angesehen. 

Wie aus der Tragfliigeltheorie bekannt ist®, kann 
man sich eine derartige Unstetigkeitsflache durch ein 
System paralleler Wirbelfaden ersetzt denken, die sich, 
ausgehend von der Hinterkante des Fliigels, nach riick- 
warts bis ins Unendliche erstrecken. Abb. 1a zeigt eine 
Draufsicht auf das System dieser Wirbelfaden im nicht 
deformierten Zustand, sowie die Zirkulationsvertei- 
lung J’, langs der Fliigelspannweite 6. Der Anderung 
von /’. um dl’, entspricht der an der Stelle « nach 
riickwarts abgehende Elementarwirbel mit der Zirku- 
lation — (dI",/dx)dx. Diese Wirbelschicht ist in hohem 
MaBe instabil. Sie wickelt sich von den seitlichen 
Enden her nach der Mitte zu spiralférmig auf und geht Abb. 1. Schema des Wirbelsystems: 
schlieBlich in zwei Einzelwirbel mit der Zirkulation J’, ) Pra errant sean: 
iiber. Bezeichnet t, die Zeit, die bis zur vollstandigen 
Aufspulung in die beiden Einzelwirbel verstreicht, gerechnet vom Augenblick der Erzeugung 
der Wirbelschicht, so hat der Querschnitt durch die zur Zeit t= 0 nicht aufgewickelte Wirbel- 
schicht zur Zeit 0<t <t) etwa die aus Abb. 1b ersichtliche Form. Die zu lésende Aufgabe 
besteht nun darin, zu jeder Zeit t die Form und GréBe der beiden Spiralen anzugeben, zu denen 
sich die Wirbelschicht bereits aufgewickelt hat, wahrend der mittlere Teil — abgesehen von einer 
Dehnung — im wesentlichen noch unverformt ist’. 


2. Das Zeitgesetz. Als Koordinaten eines beliebigen Punktes P der Spirale zur Zeit ¢ sollen 
der vom Pol 0 der Spirale nach P gezogene Strahl O P= r und der zugehérige Polarwinkel p 


1 H. Kaden, Ing.-Arch. 2 (1931), S. 140. 
2 4, Betz, Handbuch der Physik. Bd. VII. 5. 241 ff., Berlin 1927. 
3 Vel. hierzu die von H. Muttray hinter einem Rechteckfliigel im Windkanal ausgefithrten Messungen. 


Dtsch. Luftfahrt-Forsch., Forschungsbericht Nr. 824/3. 
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eingefithrt werden, letzterer von der Verbindungslinie der beiden Spiralenpole aus gerechnet 
(Abb. 1b). Man kann also die oben gestellte Frage nach der Form und Gréfe der beiden Spiralen 


so beantworten, daf man qy als Funktion der beiden unabhangigen Verinderlichen r und t be- 
| 


p= P(r, 4) « (2) 
Wie aus zahlreichen Strémungsaufnahmen und auch aus systematischen Messungen’ bekannt — 
ist, liegen die einzelnen Spiralenwindungen — zumindest in der Umgebung der Pole — sehr nahe — 
beieinander und haben angendhert Kreisform. Das Linienintegral der Geschwindigkeit langs 
eines solchen Kreises vom Radius r, d.h. die Zirkulation J’, ist gleich der Zirkulation I’, des | 
Stiickes der urspriinglichen (ebenen) Unstetigkeitsflache von x bis b/2 (Abb. 1a), die sich zur Zeit t 
innerhalb des Kreises vom Halbmesser r aufgewickelt hat?. 

Vernachlassigt man nun mit Riicksicht auf die hier als Idealisierung angenommene (an- 
geniherte) Kreisform eines Spiralenumlaufes alle radialen Geschwindigkeitskomponenten, dann 
kann sich die Zirkulation J’, lings eines Kreises vom Halbmesser r mit der Zeit nicht andern, 
d: h. J", ist fiir ein bestimmtes r unabhangig von t. Als mittlere Umfangsgeschwindigkeit langs 
der Kreisperipherie sei jetzt der Wert 


stimmt, also 


I, 
ory (3) 
definiert. Dann ist die mittlere Winkelgeschwindigkeit eines Fliissigkeitselementes im Abstand r 
vom Spiralenmittelpunkt 


Ot or 2art° (4) | 
Da es sich bei dem Aufspulvorgang um sehr viele, sehr nahe beieinander liegende Spiralenwin- | 
dungen handelt, soll in erster Naherung von der Veranderlichkeit der Umfangsgeschwindigkeit _ 
lings eines Kreisumfanges abgesehen und somit dy/dt = dq/dt gesetzt werden, so daB aus (4) | 
durch Integration in bezug auf t folgt | 
(PAG 


20 7? 


f(r). (5) | 


YP = 
Die zunachst unbekannte Funktion f(r) hat offenbar bestimmte Bedingungen zu erfiillen, die mit 
der Form der Spirale im Zusammenhang stehen. Es muf deshalb zunachst versucht werden, 
dariiber Genaueres auszusagen. 

Zu diesem Zwecke sei gemaB Abb. 2 die Flache F betrachtet, welche durch eine Spiralen- 
windung begrenzt wird und durch die sehr kleine radiale Verbindungsstrecke Ar, d. h. die Linge, | 
um die sich r bei einem Umlauf andert. Die zeitliche Anderung dF'/dt der Flache Fist gleich dem 
Flissigkeitsvolumen uMdr~ uAr, das — bezogen auf die Tiefe ,,eins‘¢ — in der Zeiteinheit durch 
Ar hindurchtritt. Es muB also sein 3 


dF | 
tATr—=—. | 
dt (6) 
Da aber unter Beachtung der oben gemachten Annahmen 
Abb, 2. Durch eine A 
Spiralenwindung ae = or 
begrenzte Flache F. on! dp 


gesetzt werden kann, so folgt aus (6) mit Riicksicht auf (3) 


nel OL eae 
Tandy edt 
oder 
Oy 7 eee | 
dr dFidt’ (7) 


Die zeitliche Anderung dF/dt der Flache Fist ebenfalls eine Funktion von rund t. Zur Abkiirzung 
soll dafiir zunichst 


1 H. Muttray, a. a. O. 


2 A. Betz, Verhalten von Wirbelsystemen. Z. angew. Math. Mech. 12 (1932), S. 164 ff 
3 H. Kaden, a. a. O. 
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dF 
a Yr) (8) 
gesetzt werden, also 
op ~—SsS, 
erry 
Differentiiert man diesen Ausdruck partiell nach t, so wird 
Go. J, op 
on ee 


Weiter folgt aus (4) mit dp/at = dq/dt durch partielle Differentiation nach r 
inky Smee aaa be 
Ordt 2x 4r (73) , 

Nach den Regeln der Differentialrechnung miissen die Ausdriicke (9) und (10) einander gleich 


sein, also 


(10) 


DO Oa ert. od tee 
oGrar ()= ik a) 
Nun ist 
: LE: r—2T,r 
@) (Le or 
ar\r2 rt =e ae 
womit (11) itbergeht in 
1 (oe 2 1 dy 
Qa\or rT, 7) w Ot” 
Setzt man noch 
i ain, Al 2 
al or rl. aa g(r). (12) 
so wird 
4 Oy 
g(r) ot = ype? 
und durch Integration folgt daraus 
1 
Sem cise 
ys) ¢ ”) 


wo T(r) eine noch unbekannte Funktion von r ist, welche die Dimension der Zeit besitzt. Nach 
aufgelést, erhalt man daraus wegen (8) 

dF 1 

dt g@k— 7] a 
Mit diesem Ausdruck liefert (7) 

og 1; 

28 — F gr) ft —2(n)). (14) 


Einen zweiten Ausdruck @q/@r erhalt man aus (5) durch partielle Differentiation nach r, 
namlich 


a a eae 
aes et Le) (14a) 
Setzt man beide Ausdriicke einander gleich, so wird 
he ip GIR. . 
Set) Eat T)) 5 zap bE (r) (15) 


Da aber unter Beachtung von (lla) und (12) 
Pea pe) is ‘ 
a Nee pale 


ist, so folgt aus (15) 
f= = Far) er (16) 
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als alleinige Funktion von r. Integriert man diesen Ausdruck nach r, so erhalt man wegen (12) 


Oe ies 1 OD) «(r) dr + € 


NM r? 2r2 Or. 


und damit aus (5) 


qe es este 1 “en dr +6. (17) 


207 a}\r3 2r? Or 
Die Zirkulationsverteilung J. = J’.(r) innerhalb der Spirale kann berechnet werden, sobald die 


Zirkulationsverteilung J”, iiber die Fliigelspannweite 6 gegeben ist!. Mit sehr guter Annaherung 
kann man dafiir setzen 


r.=T, t as ae (18) 


wobei r, den Halbmesser der Kinzelwirbel bezeich- 
net, zu denen sich die Unstetigkeitsfliche nach Ab- 
lauf der Zeit t, aufgewickelt hat (Abb. 3), und 1) 
die maximale Fliigelzirkulation. Der Exponent v 
und sein reziproker Wert 1/y hangen von dem Ver- 
haltnis o = 2a/b ab, in dem 


Abb. 3. Wirbelkerne nach Beendigung des 1 


Aufspulvorganges, CC, — | Te dx (19) 


den halben Schwerpunktsabstand der beiden Einzelwirbel vom Halbmesser ry darstellt. Die 
GréBe vy kann aus folgender Tabelle entnommen werden ?: 


p = 0 | 1e25 | 18 _ 1,75 | 2,0 | 220 | 2,5 | 2,15 
(nn 


o =0,5 | 0,6119 | 0,6888 | 0,7472 | 2/4 | 0,8245 | 0,8507 | 0,8716 
» = 3,0 ers | 4.0 | 4.5 | 5,0 
o = 0,8896 | 0,9166 | 0,9318 | 0,9474 | 0,9548 


Das Integral (17) enthalt aber noch die zunachst unbekannte Funktion t(r), durch welche erst 
die zeitliche Anderung der Spiralenflache gemaB® (13) bestimmt ist. Es darf wohl vermutet wer- 
den, da dieser Ausdruck einem der Minimalprinzipe der allgemeinen Mechanik gehorchen muB, 
indessen kann dariiber zur Zeit noch nichts Bestimmtes ausgesagt werden. Nachstehend soll 
deshalb versucht werden, einen plausiblen Ansatz fiir t(r) zu finden, der einerseits dimensions- 
richtig ist und andererseits die zur Verfiigung stehenden Randbedingungen befriedigt. 

Zunachst soll als ,, Wirbelkern‘: derjenige Bereich von Abb. 1b definiert werden, der durch den 
Winkel p = x/2 begrenzt ist. Setzt man diesen Winkel in (17) ein, so kann daraus r= r(t) als 
Kernhalbmesser zur Zeit t berechnet werden, sobald die Integration ausgefiihrt ist. Zur Errei- 
chung dieses Zieles sei zunachst ein Zustand betrachtet, in dem der Aufspulvorgang bereits be- 
gonnen hat, r aber noch einen beliebig kleinen Wert besitzt. In diesem Zustand kann man sich 
die Unstetigkeitsflache mit groBer Annaherung einseitig ins Unendliche erstreckt denken und die 
Ergebnisse der Kadenschen Theorie* zum Vergleich heranziehen. Nach dieser gilt fiir kleine r 


dF 4 fF? 
io —Y) 
und 
Oy OME 
Or 2 mrdi2” (21) 
wo 
ede 
— Sif (21a) 


Legt man zunachst, um zu konkreten Vorstellungen zu gelangen, eine alliptische Zirkulations- 


1 W. Kaufmann, S.-B. math.-naturwiss. Abt. bayr. Akad. Wiss. 1946, S. 109, und Ing.-Arch. 17 
(1949), S. 191, 

2 W. Kaufmann, a.a. O. 

3H. Kaden, a. a. O. 
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verteilung iiber die Fliigelspannweite zugrunde, so folgt aus (19) als halber Schwerpunktsabstand 
der Einzelwirbel a = bz/8, weshalb o = 2a/b = a/4. Fiir diesen Wert liefert die obige Tabelle 
den Exponenten » = 2 in (18), weshalb 


F=T,\/27—(y, (22) 
d.h. im Kern ist die Zirkulation elliptisch iiber r verteilt. La8t man nun r—> 0 gehen, so wird 
ah 
(2, peasy — i \ Hh ? (22a) 
‘ae 
weiter wegen (12) 
o 
8(T) +0 = 4 7g 1? (23) 
und wegen (13) 
dF A st r? 
peo 


Damit dieser Wert fiir sehr kleine r mit (20) iitbereinstimmt, mu 1(r) mit r schneller gegen Null 
gehen als ¢. Fiir die sich einseitig ins Unendliche erstreckende Unstetigkeitsflache wird nach (1) 


i—— Ti konst: (24a) 


Die gesuchte Funktion t(r) muB also mit r” gegen Null gehen, wo n> 3/2 ist (erste Bedingung 
fiir t(r)). 
Aus Gleichung (14a) folgt unter Beachtung von (22a) 


ay Dol oes a aE . 
(ar),.0= ee ee a ee) A 


DQ 7 27» pay 


Durch die Verbindung von (21a) und (22a) findet man 
Ly 
x= —.,, 
V2.7 


so daB 
Oy B93 vy 
uaa aett 


r>0 


Damit dieser Wert fiir sehr kleine r mit (21) ttbereinstimmt, muf 


f'(r)r+0 = 0 
sein. 
Nun folgt aus (16) wegen (22a) und (23) 
/ a 
(t= 


Soll also f’(r),.9 = 0 werden, so muf t(r) mit r” gegen Null gehen, wobei n’ > 5/2 ist (zweite 
Bedingung fiir t(r)). Nach Ablauf der Zeit t = t, ist der Aufspulvorgang beendet; der Wirbelkern 
hat dann den gréBten Halbmesser ry erreicht, und die gesamte Zirkulation /y ist jetzt innerhalb 


des Kernes enthalten (Abb. 3). In diesem Augenblick muf} also 


ie = 0 
dt T=T9 
sein (dritte Bedingung fiir t(r) . Ee . 
Aus Gleichung (13) folgt, da® 1(r) die Dimension der Zeit haben muf, aber eine reine Funktion 
vor rsein soll. Die Zeit kann also nur als Konstante darin enthalten sein — etwa ty — wahrend 
r dimensionslos in T(r) eingeht (vierte Bedingung fiir t(r)}. ve 
Die vorstehenden vier Bedingungen seien jetzt in folgendem Ansatz vereinigt: 


fy UL (25) 
> = Gera 
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wo c eine zundchst unbekannte Zahl bezeichnet. Man iiberzeugt sich leicht, daB dieser Ausdruck 
die oben gestellten Forderungen erfiillt. 
Setzt man nun (25) in Gleichung (17) ein, so wird 


het Cy [tle LO aay ee C 26 
e= reat eile 2r? ar} ; = or nee 30) 
/ 5} 
/ To 
y1-@) 
Die weitere Rechnung soll auch hier fiir den Fall elliptischer Zirkulationsverteilung durchgefithrt 
werden, obwohl dies grundsatzlich auch fiir jede andere J’,-Verteilung méglich ware. Zunadchst 


wird das in (26) auftretende Integral behandelt. Unter Beachtung von (22) erhalt man dafiir 
nach EHinfithrung des dimensionslosen Wertes 


ow (27) 


das Integral 


ea - Tees oT 
ri Irth QP 


UL Vo7 — Te ie (Weyg r7/2 ed Slt ke 7 
Is [( ii z 


ay I 5 - i r . ; T = = 
=7s/@ r) Sees dr. 
Der Ausdruck (3 — r) p= kann mit sehr guter Naherung durch die Parabel 2 |r ersetat werden, 


so daB man erhalt 


Viel, (pear = r2 
Mit diesem Wert des Integrals liefert (26) schlieBlich 


Tit Lot /7 = ! Gs (28) 


=f r dr— Ig _ d(7?) es oy 72. 


r2 
Un 


oes Doo heats TUT 6 
Setzt man nun fiir den Kernrand y= 7/2 (s. oben), so kann r als Funktion von t angegeben 
werden, sobald die drei noch offenen Konstanten c, C und t) bestimmt sind. 

Fiir sehr kleine Werte von r gilt das Zeitgesetz (24a), worin die Konstante zunachst allerdings 
nicht bekannt ist. Aus dieser Gleichung folgt t= 0 fiir r—0 und (dt/dr),,)9 = 0. Diese Be- 
dingungen sollen auch hier gefordert werden. Weiter muB fiir r= ry, d. h. bei Beendigung des 
Aufspulvorganges, t= t, und (dr/dt),_,, = 0, bzw. (dt/dr),_,,== 00 werden. Diese vier Be- 


dingungen werden — wovon man sich leicht itberzeugt — durch den Ansatz 
our | lye ay ais) (29) 
| ok 


befriedigt. Entwickelt man ihn in eine binomische Reihe, so wird 
eo Saha a 
3 e io) ro | é 


(t)r+0 = = ty (zal , (30) 


tee i, 


woraus fiir kleine r folgt 
7o 


Der obige Ausdruck (29) soll nur zur Bestimmung der Konstanten in Gleichung (24a) dienen, die 


lediglich fiir kleine r gilt. Das wirkliche Zeitgesetz wird dagegen erst aus Gleichung (28) ab- 
geleitet. 


Man erhalt also aus (28) mit m = 7/2 fir kleine r wegen (22a) 
es (2/SV OP tareletae ; ; 
oF Pari mr? aie CS) 
Fir r= ry ist t= t) und J;= 1). Setzt man diese GréBen in (28) ein, so wird 


mI t 


2 Qur2 ante (32) 
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oder 


(pe ee 


2: Diccmom 
SchlieBlich liefert die Verbindung von (31) und (32) sofort 
1 


Weel 
c= 3 /2—>=— 0,029. 


Mit g = 2/2 und den Konstanten ¢ und C folgt dann aus (28), wenn noch die dimensionslose 


GroBe 


t = 
= (33) 


to 
eingefiihrt wird, unter Beachtung von (22) 


O= = /27—F +0,058 YI — 7 — 1 


oder 
ee ae (34) 
Wie man sieht, wird 
i); Goce Ris Ray 
ec (120,058) mo Se, 


entsprechend der Bedingung (30), und t= 1 fiir r= 1. Differenziert man in (34) ¢ nach 7, so 
stellt man fest, da®B in der Tat (dt/dr),.9 0 und (dt/dr),_,, = 00 wird, wie oben gefordert 
wurde. 

Der gréBte Kernhalbmesser kann berechnet werden, wenn der halbe Schwerpunktsabstand a 
der beiden Einzelwirbel vom Halbmesser ry aus (19) bestimmt ist. Man erhalt dafiir 1 


i= ay/3/4 : (35) 

In dem Zeitgesetz (34) ist somit nur noch die Aufspulzeit tj unbekannt. Bevor an deren Ermitt- 
lung herangegangen wird, soll zunachst als Kontrolle der Gleichung (25) noch ein anderer Ansatz 
fiir die Funktion t(r) gemacht werden, namlich 
SROs 

jl—r ; 
der ebenfalls die oben fiir r(r) aufgestellten Bedingungen befriedigt. Fiihrt man in Gleichung (17) 
wieder die Zirkulation J’, aus (22) ein und beachtet die Abkiirzung (27), dann lautet jetzt das i? 
(17) auftretende Integral 


SSS aay aris 
ja Bue I Le ioe ro di 


(36) 


a 7 rj 27 rs yor—r/ /1—r 
Be Piatoic ras. dr = 74 [@ )|/ ees oe 
2are} V2r—P VI—F 2a 2—ry1—rF 
Setzt man wieder wie oben 
a 
rene r = 
Sends TANT 
so wird 
— Lyte aR aN Ty t ¢ [Js 3 4) ae __ 9s 
sia mre | \r—?? oie mT? li fae ia ata le 27). 
womit Gleichung (17) lautet 
Iyt Ly ty € ae : = ; 
°= sain an [2 jr — 7? +arc sin (1 — 27)] +C. (37) 


Zur Bestimmung der Konstanten c und C dienen die gleichen Bedingungen wie oben. Fir kleine 
r gilt Gleichung (30), auBerdem ist t= 1, fiir rr) und [}=— J. Mit p= m/2 folgt dann aus 
(37) unter Beachtung von (22a) fiir kleine r 


1 W. Kaufmann, a. a. O. 
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Gia to tose a 
ae rire Qur2 2 ! 
UNG s cur r=" Te 


wT Lot) , Lote m | C 

Oe To gta Tee | ee 

Durch Verbindung beider Gleichungen erhalt man 
2—1 


Tt 


— 0,0182 , 


— & Dy to 
Ga ae 0,971 eT 


Damit folgt aus (37) mit p= 7/2 
o= 4! + 0,0182 [2 ¥F— # + are sin (1 — 27)| — 0,971 
leone 


oder unter Beachtung von (22) 
72 —= 5 
i= a (0,971 — 0,0182 |2 yr — r? + are sin (1 — 27)|). (38) 
27— 7 


Nachstehend sind die Werte t, die sich aus (34) und (38) ergeben, miteinander verglichen. Man 
erkennt daraus, da® praktisch vollkommene Ubereinstimmung besteht (Rechenschieberable- 
sung). Damit liegt das Zeitgesetz t =t(r) fest. sobald r) nach (35) berechnet und die Aufspulzeit t) 
bekannt ist. 


7F=|0/| 0,1 0.2 SI 053 0,4 0,5 0,6 0,7 | 0,8 0.9 | 1,0 
t= |0 | 0,022 | 0,063 | 0,119 | 0,189 | 0,275 | 0,375 | 0,493 | 0,630 | 0,796 | 1,0 (34) 
t=|0]| 0,022 | 0,063 | 0,119 | 0,190 | 0,275 | 0,376 | 0,495 | 0,630 | 0,796 | 1,0 (38) 


| 


3. Die Aufspulzeit t). Die in Ziff. 2 durchgefiihrten Untersuchungen sind rein kinematischer 
Art und dienen lediglich zur Bestimmung der Form des Spiralenkerns; iiber die GréBe der Auf- 
spulzeit vermégen sie jedoch nichts auszusagen. Letztere hangt offenbar von den Kraften ab, 
die zur Aufwicklung der ursprimglich ebenen Unstetigkeitsflache fiithren. Die ins Einzelne 
gehende Untersuchung dieses Vorganges ist recht schwierig und bisher nur fiir das erste Stadium 
gelungen?. Es kann deshalb hier lediglich eine mehr summarische Betrachtung dariiber angestellt 
werden, vorbehaltlich einer spateren Nachpriifung durch den Versuch. 

Im ganzen gesehen sind die Krafte, die beim Aufspulvorgang eine Rolle spielen, der Auf- oder 
Quertrieb 


A= OU, [T.dx 
(6) 


und der induzierte Widerstand 


0 = GPs Giles 
Wee Tae ( Sis 59" 


(6) (6) 
die gemaB vorstehenden Gleichungen von der speziellen Form der Zitkulationsverteilung itber 
die Spannweite b (Breite der urspriinglichen Unstetigkeitsflache) abhingen. Fiir den Fall, da 
dI’,./dx = 0 fiir alle xist (ebenes Problem), gibt es keine Wirbelschicht (vgl. Ziff. 1) und auch kei- 
nen induzierten Widerstand. Bekanntlich kénnen A und W; auch wie folgt geschrieben werden: 


At cagl , W;= Cwi GF , 


wo c, die Auftriebs- und c,,; die induzierte Widerstandsziffer bezeichnen, wahrend q= (0/2) v2, 
den Staudruck und F die ,,Fliigelflache** darstellen. Im Falle elliptischer Zirkulationsverteilung 


c, F ; 
TU ist *¢,,, = = ree Bildet man nun den Quotienten 
J 


A Ca a Ob 


Wi: = Cw; = Ca F 4 (39) 


or Rn By Gly Op 
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so erkennt man, daf dieser Quotient bei gegebener Fliigelstreckung b2/F mit kleiner werdendem 
c,immer groBer wird und schlieBlich mit c, > 0 gegen unendlich geht. Bei nicht elliptischer, sich 
mehr dem Rechteck nahernder /’,-Verteilung ist 


CLF 
Daa ae © 
wo w> 11, weshalb 
: A zw 6 
Wi Wa F . (39a) 


Der Faktor y wird um so gréer je gréBer die Fliigelstreckung ist. AuBer diesen beiden Kraften 
spielen fiir den Aufspulvorgang zwei charakteristische Lingen eine Rolle: Der Abstand 1,, wel- 
cher das Ende der Aufwickelung, gerechnet von der Fliigelhinterkante aus, bestimmt und der 


Abstand a’ = 2 b— a des Schwerpunktes S der beiden Einzelwirbel von den Fliigelenden am 


Ende der Aufspulzeit t) (Abb. 3). Uber die Lange J, ist zur Zeit nichts Genaueres bekannt. An 
einem elliptischen Fliige] von der Streckung b?/F = 5 hat Muttray? die Lage der Randwirbelkerne 
fiir c, = 1,07 ausgemessen und kommt zu dem Schluf$, daf der Aufrollvorgang vermutlich bei 
einer Linge |) = 6 beendet ist. Auch aus Luftaufnahmen® sind die beiden aufgewickelten Ein- 
zelwirbel bereits in kurzer Entfernung hinter dem Fliigel erkennbar. Bei gegebener Zirkulations- 
verteilung J’, ist a’ eine ganz bestimmte Lange, die nur von der Form dieser Verteilung, nicht 
aber vor der absoluten GréBe von I’, abhaingt. So ist z. B. bei elliptischer /',-Verteilung wegen 


(19) a= bz/8 und somit — b(1 —a/4). Entsprechende Werte kénnen fiir andere J/’,- 


Verteilungen angegeben werden. Das Verhaltnis |,/a’ wird sich also bei einer bestimmten Flii- 
gelform nur in dem Mae andern kénnen, in dem sich c, andert. Geht c, gegen Null, dann muf 
1,» co gehen. Denkt man sich als andern Grenzfall einen ,,Hufeisenwirbel*‘, so ist a’ = 0 und 
I, = 0, denn in diesem Falle gehen die beiden ,,Randwirbel* direkt von den Fliigelenden nach 
riickwarts als Fortsetzung des ..tragenden Wirbels‘‘, der den Auftrieb erzeugt. Setzt man nun 
die beiden Verhaltniswerte A/W; und |,/a’ einander proportional, so erhalt man wegen (39a) 


% Ob 
Seer baae 0 
I re see (4 ) 


und erfillt damit jedenfalls die beiden Grenzbedingungen |, = 0 fiir a’ = 0 (Hufeisenwirbel) und 
1, = co fiir c, = 0 bei vorgegebener [’, — Verteilung. Da aber 


iie=tly Ue 
ist, wo v, die als gleichfé1mig angenommene Fortschrittsgeschwindigkeit des Fliigels bezeichnet, 
so ergibt sich aus (40) mit a’ = : b—a als Aufspulzeit 
b/2—a a B 
ae mi ce Fe 


to (41) 
Dieser lineare Ansatz hat zunachst nur hypothetische Bedeutung und bedarf hinsichtlich seiner 
Brauchbarkeit einer systematischen Nachpritfung durch das Experiment. Insbesondere kann 
iiber die Konstante x zunachst keine Angabe gemacht werden. 


4. Vergleich der vorstehenden Theorie mit Versuchsergebnissen. Kine Priifung der Theorie 
an Hand von Versuchen bereitet insofern Schwierigkeiten als bisher weder die Aufspullange /, 
noch die Gré®e der Wirbelkerne zu einer beliebigen Zeit 0< t <t) genau gemessen werden konn- 
ten, da sich solchen Messungen erfahrungsgema8 erhebliche Schwierigkeiten entgegenstellen. 
Eine Vergleichsméglichkeit bieten jedoch die sehr sorgfaltigen Druckmessungen, die Muttray * 
an einem rechteckigen Fliigel mit der Fliigelstreckung b?/F = 5,36, der Auftriebsziffer ¢, = 0,505, 
der Fliigelflache F = 0,1439 [m?] und der Halbspannweite 6/, = 0,439 [m] bei der Anstrémungs- 
geschwindigkeit v= 32,9[m/sec] durchgefiithrt hat. An diesen Untersuchungen interessiert 
fiir das vorliegende Problem besonders die Lage der Kernschwerpunkte in verschiedenen Ab- 


standen hinter dem Fliigel. Fat man — wie oben geschehen — die Kernrander zu einer belie- 
1 Vgl. hierzu: J. Hueber, Z. Flugtechn. Motorluftschiff. 1933, 5. 249—251 und 269—272. 
2 H. Muttray, Luftfahrt-Forsch., 1935, S$. 28—37. 
3 H. Lueder, Jb. 1942 dtsch. Luftfahrtforsch. SI 103/105. 
4 


H. Muttray, Dtsch. Luftfahrtforsch., Forschungsbericht Nr. 824/3. 
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bigen Zeit angendhert als Kreise auf (Ziff. 2), so fallen die Schwerpunkte S; dieser Kerne nicht 
mit den Mittelpunkten M der Kernrander zusammen, sondern liegen um eine gewisse Strecke e 
exzentrisch zu diesen (Abb. 4).1 Bezeichnet nun x, den Abstand des Kernschwerpunktes S, von 
der Symmetrieachse und x die Abszisse der Zirkulation I.=T,, die sich im Kreise vom Halb- 
messer rzur Zeit t bereits aufgespult hat (vgl. auch Abb. 1a), so gilt fiir x, mit der Abkiirzung 


(x, — =)? = f(x) (42) 


folgende Beziehung * 


272 = f(x) — a +a Ja? +2 f(x). 


Die Auflésung dieser in f(x) quadratischen Gleichung liefert 


fla) = 2a +9) (2 FPP Ae He, (43) 


so daB wegen (42) fe 
n= = +O) fees 


; | 


Abb, 4. Wirbelkerne zur Zeit t< to. Abb. 5. Verlauf der Fliigel- und Kernzirkulation. 


fiir jedes beliebige r berechnet werden kann, wenn J",(x) gegeben ist. Nachstehend soll diese | 


Rechnung am Beispiel des Muttrayschen Versuchsfliigels erlautert werden. 
Zuniachst werden die Fliigelzirkulation J’, und die Kernzirkulation J’, durch folgende Ansatze 
dargestellt? [vgl. auch (18)]: 


eee Z| 
To 


Die GréBe o = 2a/b (Ziff. 2) kann einer Arbeit von Hueber® entnommen werden, welcher die 


und 
l/y 


jeeaaar 


(45) 


Zirkulationsverteilungen fiir eine Reihe von doppeltrapezférmigen Fliigeln bei verschiedenen | 
Fliigelstreckungen berechnet hat. Danach ergibt sich fiir einen Rechteckfliigel vom Seitenver- | 
haltnis 1:5 der Wert o = 0,864, dem nach der Zahlentafel auf S.4 der Exponent » = 2,66 


entspricht. Da das Seitenverhaltnis des Muttrayschen Fliigels nahezu das gleiche ist (1:5,36), 


soll hier mit dem vorstehenden Wert von » gerechnet werden. Der damit entstehende Verlauf | 


der Funktionen [°,/J7, und J/I’) ist aus Abb.5 ersichtlich. 
Aus Gleichung (40) folgt nun als Aufspullange mit den oben angegebenen Zahlenwerten und 
py = 1,043 


b Te. 5,36 S 
Cree mean aos = oe ae 
oder 
I, 
aa 2,18 x 


Beachtet man, dais Muttray (s. oben) fiir eine Auftriebsziffer c, = 1,07 bei einem elliptischen Flii- 
gel von der Streckung b?/F' = 5 die Zahl 1,/b = 1 angibt, so darf wohl aus Gleichung (46), welche 
we W. Kaufmann, a. a. O. 


2 W. Kaufmann, a. a. O. 
3 J. Hueber, a. a. O. 
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— 0,505 (also nur halb so gro) enthalt, gefolgert werden, daB die Konstante x sich nicht viel 
von ,,eins‘‘ unterscheiden wird. Fiir die weitere Rechnung soll deshalb vorlaufig x — 1 gesetzt 
werden. Damit wird fiir den vorliegenden Fall wegen b = 87,8 [cm] 


lp 87,8 2,18'== 192 [em] 
und 


1 
i= id = 0,058 [sec] . 
Weiter ist 


a= 0 ) = 0,864 - 43,9 = 38 [em] 


und wegen (35) 
ry = 35,4[em]. 

Fiir eine beliebige Linge 0 </ < J, gilt das Verhiltnis 

bone 2 

ee Rant 
Tragt man nun nach (34) r— r(t) auf, so kann mit Hilfe dieser Figur der zu jeder Lange I ge- 
hérige Verhaltniswert r= r/r, bestinmt werden, womit auch der zul bzw. t gehérige Kernradius r 
bekannt ist. Streng genommen gilt Glei- 
chung (34) nur fiir elliptische J',-Vertei- 
lung. Da sie jedoch nur die VerhAltnis- 
gréBen r bzw. t enthalt, wird man wohl 
keinen groBen Fehler begehen, wenn man 
sie auch fiir die hier maBgebende nicht- 
elliptische J’,-Verteilung anwendet. 

Um nun den gesuchten Schwerpunkts- 
-abstand x, des zur Zeit tim Abstand | hin- 
ter dem Fliigel bereits aufgewickelten Wir- 
belkerns zu erhalten, beachte man (Abb. 5), 
daB die dem Verhaltniswert r entsprechen- 
de relative Zirkulation J°/I) gleich I’,/I, 
sein mu (s. oben). Die zu J’,/Iy gehirige 
Abszisse ist daher 2 x/b = 1— T/T; womitx Abb. 6. Abszisse des Kernschwerpunktes als Funktion 
und somit nach (43a) auch x, berechnet der Aufspullange. 
werden kann. 

Die Muttrayschen Druckmessungen liefern die Lage dieser Schwerpunkte fiir drei verschiedene 
Langen | hinter dem Fliigel, und zwar ist 


U 


= 


1 1 


1 es ! | 1 == ee 
0 20 40 60 60 0 120 WO 0 ocCm 200 


fiir I= 4,9 20,2 61,0 [em] 
x, = 41,9 41,0 39,8 ,, 
xt = 42,8 40,8 20,20 


Die entsprechenden Werte, welche sich nach der vorstehenden Theorie ergeben, sind als x, 
darunter geschrieben. Abb.6 veranschaulicht schlieBlich die vollstandigen Ergebnisse der Rech- 
nung. 

Der Aufspulvorgang ist im Abstand /, = 192 [cm] hinter dem Fliigel beendet, wozu nur die 
geringe Zeit von t, = 0,058 [sec] benétigt wird '. 


(Eingegangen am 8. Juni 1950.) 
Anschrift des Verfassers: Professor Dr. W. Kaufmann, Miinchen 2, Walther v. Dyck-Platz 1. 


1 Bei der Durchfiihrung der Zahlenrechnungen hat mir Herr Dipl.-Ing. R. Wurzbach wertvolle Hilfe 
geleistet. 
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| 
Theorie des Hochdruckmefrohres mit ausmittiger Bohrung. | 
Von W. Wuest. | 

| 


1. Einfiihrung. Dickwandige Rohre mit konzentrischer Bohrung werden durch einen inneren} 


Uberdruck einem rotationssymmetrischen Spannungszustand unterworfen. Bei ausmittiger’ 


Bohrung geht diese Symmetrie verloren. Der Angriffspunkt der Axialkraft (0) in Abb. 1) fallt: 
nicht mehr mit dem Schwerpunkt S der Querschnittsfliche zusammen, so daB ein Biegemoment | 
auftritt, das eine Kriimmung des urspriinglich geraden Rohres bewirkt. O. Hoffmann und 
W. Francke’ haben diese Eigenschaft wohl erstmalig erkannt und zur 
Grundlage eines HochdruckmeBverfahrens gemacht, bei dem die Aus- 
lenkung eines geraden, exzentrisch gebohrten Rohres zur Druckanzeige 
benutzt wird. Gegeniiber der itblicherweise angewandten kreisférmig ge- 
kriimmten Réhrenfeder von ovalem Querschnitt (,,Bourdonfeder‘‘) hat 
das neue Mefiverfahren bei sehr hohen Driicken (iiber 1000 at) den Vor- 
teil, daB die Beanspruchungsspitzen niedriger sind und ein Sicherheits- 
schutz gegen Zerknall der Hochdruckfeder leicht anzubringen ist, indem | 


4Z 


mierten Fliissigkeiten eine erhebliche potentielle Energie besitzen. 
Eine elementare Theorie des Vorganges wiirde davon ausgehen, 
daB die Kriimmung des Rohres den Wert 

were Ti 

| ear. 
hat. Da der Abstand des Schwerpunktes S vom Kraftangriffspunkt 0, 


gegeben ist durch s= r? e/(r? — r?), wo e die Exzentrizitat bedeutet. 


(1) 


0 
erhalt man fiir das Moment M 
Abb. 1. Zustandekommen des ir TPG : 
Kriimmungseffektes bei ex- Ses = er bi (2) 
zentrisch gebohrten Rohren fe ho 


unter Innendruck. S Schwer- es : 5 
punkt derQuerschnittsflache. Das auf den Schwerpunkt bezogene Tragheitsmoment Js ist 
Die Pfeile deuten an, daB auf 
der schwicheren Seite die Js= 
Tangentialspannungen und 2 ae 
damit auch die axiale Quer- ; A ; 
kontraktion stiirker als auf Damit liefert die elementare Rechnung fiir die Kriimmung j, 

der gegeniiberliegenden See 

Seite ist. : 4e To TT P 


Jo (=A) (Ht =) — 4 rer? E (4) 


Versuche, iiber die noch berichtet wird, haben gezeigt, daB die Kriimmung tatsachlich nur etwa 


Te GU CL uAT 


die Halfte der nach Formel (4) berechneten betriagt. In der Tat zeigt eine genauere Beschreibung _ 


des Vorganges, daB auf der schwacheren Seite, wie in Abb. 1 angedeutet ist, die Tangential- 
spannungen gréBer als auf der gegeniitberliegenden Seite sind. Damit ist auf der schwacheren 
Seite auch die axiale Querkontraktion gréBer. Hierdurch wird also der urspriingliche Kriimmungs- 


effekt abgeschwacht. 


iiber das MeBrohr ein Schutzrohr gestiilpt wird. Ein solcher Sicherheits- 
schutz ist beim Experimentieren mit hohen Driicken unerlaBlich, weil | 
die Me&feder an der Grenze ihrer Festigkeit arbeitet und die kompri- | 


— t _y4 | 
4 ahd) pat (3) 


Bei der genaueren theoretischen Behandlung dieses Problems kann man davon ausgehen, 


daB offenbar abgesehen von der unmittelbaren Nachbarschaft der Rohrenden ein ebener Form- 
dinderungszustand vorliegt, so das die Lésung in bekannter Weise von einer Spannungsfunktion 
abgeleitet werden kann, welche die Bipotentialgleichung befriedigen muB. Unabhangig von- 
einander sind von M. Bricas* und E. Weinel® Formeln zur Berechnung des exzentrischen Kreis- 
ringes unter innerer oder auBerer Druckbelastung aufgestellt worden, die aber voneinander 
abweichende Ergebnisse liefern. Die von M. Bricas angegebene Lésung, die sich als nicht ganz 


* DRP. 743 613 (1937), Franz. Pat. 846 870 (1938). K. Palzer, Armaturen fiir den Hochdruckeeriite- 
bau. Maschinenelemente-Tagung Diisseldorf 1938. VDI-Verlag 1940. : 
* M. Bricas, La théorie de V élasticité bidimensionelle. Athen 1936, 
° E. Weinel, Z. angew. Math. Mech. 17 (1937), S. 284. 


XIX. Band 1951. Wuest: Theorie des HochdruckmefBrohres mit ausmittiger Bohrung. 13 


richtig erweisen wird, ist neuerdings auch in den 3. Band des Lehrbuches von L. Féppl* auf- 
genommen worden. 


2. Spannungszustand in einem exzentrischen Rohr unter Innendruck. Bei einem ebenen 
Formanderungszustand kénnen die Spannungen o,, o, und t,, von einer Spannungsfunktion F’ 
abgeleitet werden, wobei 
__ @F arr 


ar 
Oy = On ’ Oo; = oy? p) Txz — dy Oz 


ist. Die Funktion F muf die Bipotentialgleichung A’ A’ F = 0 erfiillen, wobei A’ = 02/0 y? + 02/az2 
gesetzt ist. 
Die axiale Spannung ist gegeben durch 


On = + (a, +0,) + E(k—jy). (5) 


-Hierbei driickt k die axiale Dehnung und j die Kriimmung aus, deren Vorzeichen wir von vorne- 
herein so wahlen, da jim vorliegenden Problem positiv wird. Wenn ¢, und q, beliebige Potential- 
funktionen sind, so daB A’p= 0, ist 


F= (y? +2) (1+ %2 (6) 
eine Lésung der Bipotentialgleichung. 

Fir die Berechnung des exzentrischen 
_Rohres empfiehlt sich der Ubergang zu bipo- 
laren Koordinaten. Die bipolaren Koordina- 
ten 7/,¢ sind mit den kartesischen Koordinaten 


-y,2 durch folgende Beziehung verkniipft?: 
1 1 4 
y+iz=eStg) (nti. (7) 


Durch Trennung in Real- und Imaginarteil 


vergibt sich 


AZ 


sm ENG La 
t= Coj 7 — cos € 7 (8a) 
sin € 
Sas Cof 7 — cos ° oP 
Ferner ist 
Col + cos € 
Die D eee OS 
pate of 7 — cos * (9) 


Die Linien 7, = konst sind Kreise mit dem ee Se ee ee eee 
Halbmesser r—c/Giny und dem Abstand Bare Sei eicoe eeerce ae Oe ere 
a=ctgy7 des Mittelpunktes vom Koordi- 

‘natennullpunkt (y= z= 0). Das exzentrische Rohr ist begrenzt durch die Kreise ry = ¢/Cin yp 
und r,=c¢/Giny,, und die Exzentrizitat e betragt 


e= ¢ (tg 4, — Ctg yo) = Jr? fe [r+e, (10) 


-woraus man findet 
Cipolla) PV 2 
Ave? e® == (r2 —-r? — e*) 4r2r 


a Ht) 
: ( 
An den Randern 7 = 7) und 77, sind folgende Bedingungen zu erfiillen: 


N=? T= 0 Nila Age 

== 4: Tat == 0 Og a 
Da es sich um eine zweifach zusammenhangende Flache handelt, sind noch die drei Eindeutig- 
keitsrelationen zu erfiillen, die physikalisch den Sinn haben, da bei einem Umlauf auf der 


Ringflache sowohl die Verdrehung y,z< als auch die Verschiebungen v, und v; eindeutig sein 
missen. 


1 L. Féppl, Drang und Zwang. 3. Bd., S$. 145. Miinchen 1947, ; 
2 M. Bricas und L. Féppl wablen statt dessen die Transformation y + iz = ¢ Ig 2 (7 +716), was 


lediglich auf eine andere Nummerierung der ¢-Kreise hinaus kommt. 
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Die Spannungssumme s’ = 0, -|-o; = A’F stellt eine Potentialfunktion dar. Wir bezeichne | 
mit t die zu s’ konjugierte harmonische Funktion. Wenn wir beriicksichtigen, daB die an den 
beiden Randkurven angreifenden Krafte im vorliegenden Fall je fiir sich ein Gleichzewichtdl 
system bilden, liefern die Eindeutigkeitsrelationen folgende Bedingungen?: 


j 
fdt= 0, (12)) 
f(s’ it) d(y +iz)=0. (13)) 
Die Gleichung (13) liefert fiir Real- und Imagindarteil je eine Aussage. Der Integrationsweg auf! 
der Ringfliche ist dabei so zu wahlen, daB die Rohrseele einmal umschlungen wird. 
M. Bricas? hat folgende Liésung des Problems angegeben, die von L. Féppl® tibernommen) 
wurde: 
6, = A— B Gin? n, 
oc= A+B Gin? y +2 B(1— cos oj n), (14))) 
t,o 0. | 
Die Spannungssumme s’ = 6, -+o;¢ wird also 
s'=2A+2 B(1—cos€ Co} x). 


Die hierzu konjugierte Funktion lautet t= — 2 Bsin( Giny. Die Grobe t ist also eindeutig; 
und damit die erste Bedingung (12) erfiillt. Die Bedingung (13) geht iiber in: 
f(s’ fit) d(y Liz) = G2 (A+ B— BGof (y +i) d(y tia. 

Fiihrt man «= (y-+iz)/e ein, so wird oj (y-+i6) = («2 -+1)/(a2— 1). Bei der Integration) 
in der komplexen a-Ebene ist der singulare Punkt x1 zu umschlingen. Das Umlaufintegral| 
fi (a? +-1)/(a? — 1) da verschwindet aber nicht, da a= 1 ein Pol erster Ordnung des Integranden) 
ist. Die Lésung (14) ist also zu verwerfen. Pi 

Es ist fiir das Folgende zweckmafig, zu einer komplexen Spannungsfunktion F* — F +-1F 
iiberzugehen, wobei an Stelle der Beziehung (6) folgende tritt: 


Fe = (y2 42%) gt+ gi (6a) 


Pi=GPitiv, und Go—gHtigQg. 


Damit wird die Spannungssumme mit ihrer komplexen Erganzung 


Pq. 2 / es 5 d yy 
s+it= A P*—49)+4(¥ FS Gane (15)) 
oder unter Beriicksichtigung von (7) j | 
Hoge 06 , soy UR LC | 
sit 4g%—4Gin(y +id) fe (15a) 


Die von E. Weinel* angegebene Lésung legt es nahe, fiir die Potentialfunktionen yj und gs! 
folgende Ansatze zu wahlen: 


Pi = A, + 4g (y +i) + As Cof (+16) + Ay Gin (+76), (16a}) 
| 


Py = & (By (4 +1) + B; of (7+ iL) + By Gin (+7). (16b) 
Die A; und B; sind hierbei reelle Konstante. 


Durch Einsetzen von (16a) in (15a) erhalt man fiir die Spannungssumme mit ihrer komplexen | 
Erganzung 


fit (44, +245) £44, (y +-i0) + 4g [4 Cof (9 +52) 
— 2 Go} 2(7 4-1 0)]+ A, [4 Sin (ny +-¢) — 2 Sin 2 (y+ i2)] (17) 
Aus der ersten Umlaufbedingung (12) folgt sofort A,—0, da sich sonst t bei einem Umlauf um. 
8 A, andern wiirde. Die Kindeutigkeitsrelation (13) wird von den Koeffizienten von A, und A,| 
einze!n fiir sich erfiillt, wie nachstehend gezeigt wird. Man fiihrt hierzu wiederum die komplexe 


Veranderliche a = (y --1C)/cein, so daB Gin (7 + i) = 2 a/(a? — 1) und Gof (y + if) = (02 + 1)} 


‘ C. Pate le! und R. Grammel, Technische Dynamik. S. 102. Berlin 1939. 
a. a. O. 

SSa,ae S 

SasanOwowl2o. 
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(«?— 1). Das komplexe Umlaufintegral (13) lautet dann fiir die beiden Koeffizienten von A, 


und A, 
A (at — 1) — 16 x? — f |o® —2—30 a® + 2 — 34) 
eam = bP ean tea |=. 
8 a (a* — 1) — 8a (a? + 1) A A] 
p (a? — 1} = Sle y Cait eens 


In beiden Fallen verschwinden die Umlaufintegrale um die Stelle x—1, weil der Integrand keine 
Pole erster Ordnung enthalt. Von dem Verschwinden der Pole erster Ordnung iiberzeugt man 
_ sich leicht, wenn man 4 = «—1 im Integranden einfiihrt. Damit ist also erwiesen, daB die 
Ansitze (6a), (16a) und (16b) mit 4,—0 die Eindeutigkeitsbedingungen erfiillen, und es ist 
noch zu priifen, ob mit diesen Ansatzen auch die vier Randbedingungen erfiillt werden kénnen. 
Nach Biezeno-Grammel sind die Spannungen fiir bipolare Koordinaten gegeben durch 


tt=— (Coj 4 — cos €) (Go) % — cos C) a + sin ¢ + Gin oF ; (18a) 
C= = (Coj 7 — cos) Goi 7” — cos C) a Sin 7 a + sin € oF ; (18b) 
Oo. = * Go} 74 — cos) (Goi 1 — cos€) a + Gin 7 i sin € a : (18c) 


Wenn man F, d.h. den Realteil von F* hierin einsetzt, erhalt man zunidchst fiir die Schub- 
spannung 
t = — (of 7 — cos C) [— 2 (Coj 7 +-cosl) Gin y cos f A, — | 
— (Coj 27 +2 oj 7 cos C) sin € A,+ B, sin ¢ — 2 (of 7 — cos) Sinn sing B, , (19) 
— (Gof 2 4 — 2 Woj 7 cos €) sin € By). | 
In der eckigen Klammer kommen also Glieder vor, die mit sin€ und sin€ cos¢ behaftet sind 
und an den Randern fiir sich verschwinden miissen. Aus dem Verschwinden des Gliedes mit 


sin¢ cos¢ an den beiden Randern findet man A,— B, und A,—B,. Das Verschwinden des 
Gliedes mit sin¢ an den beiden Randern liefert die Beziehungen 


ae) = 1 Gin (% + m) 9 
eae Ce ey 
Lge Ep = 1 Go} (m +m) p (20b) 


260f %o—m) 


In gleicher Weise findet man fiir die beiden iibrigen Spannungskomponenten o, und o;, wenn 
man A, B, und A,= B, setzt, 


O, = (2A, + A;) — Sin y (Coj 7 cos C) B, 
— 2 Gin y cos 6 (Gin 2 7 A, + Cof 2 A,) + (Coj 27 As+ Gin 2 7 A,). (21) 
Oc = 0, +2 (Cof 7 — cos) [Gin B, +2 cos C(Cof 27 A, + Gin 2 y Ay). (22) 
Setzt man zur Abkiirzung noch 24,14A,—= A und B,= B und fihrt nach Weinel 
Mo =(o + 7y)/2 und fy =(No— N)/2 sowie = — My 


ein, so erhalt man fiir die drei Spannungskomponenten 


ae Gin 2 1 eg : 
a = B E (Sin i Sane Sinn pale ; (23a) 
; in 2 
oc = 0, +2 B (Coj y — cos¢) in +e cose j (23b) 
Cof 2m 
T t= — B (Coj yn — cos¢) sing (1 — os 3) ; (23¢) 
Durch die Randbedingungen o,, = — p fiir 4= 1% (bzw. w=") und o,= 0 fiir n=, (bzw. 
js = — fy) werden die Konstanten A und B bestimmt zu 


Els His (0), fe) ORS, 
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1 1+ Ctg 2m, Sin2n, 
, ers i 24 
A=24A,+ 4, ' 2 Ginty + Sin?y, *’ ee 


Ctq 2 my 
Pag # eee 20 
pee! Sin? 49 + sin? y, ~~ 2) 


Damit nimmt die Tangentialspannung o; am inneren Rand 7) und am duBeren Rand 7, folgende 
Werte an: 

Innenrand 7 = 7%: 

Cig 2 wm, 
Sin?) + Gin? 7, 


OG == Peep 


AuBenrand 4 = 7,: 

Ctg 2 
Pees) 1 
: P Sine No + Sin? n, 
Die Héchstwerte oj und o, hiervon ergeben sich an der engsten Stelle €—ma. Macht man davon 
Gebrauch, daf 


(0) 


(Coj 4, — cos C) (Sin Hy, — Tq 2p, cos). (26b) 


Mo(0,€, 
4 0(0.) a9 ees r? —r? —e? 
SUM Yo = 7% Co} y= oun ro = 
c ’ rj —r? + e 
Sm Vy = eal V4 = oar = 


Ctg 2 u,= i epi het 


DENG 


ist, so erhalt man schlieBlich folgende For- 
meln fiir die gréBte Tangentialspannung innen 
und aufBen; 


1 % ie 2% +e 


r2 


o,= 2p, —3, (1 tee) eT) 
L il 0 


eer 


Die Formeln (27a) und (27b) stimmen mit den 
von Weinel abgeleiteten iiberein, wahrend die 
Formeln (26a) und (26b) bei Weinel Druck- 
Abb. 3. Beiwertfunktion m, (0, ¢) fiir die Berechnung der gréSten ‘ ; 
Beanspruchung eines exzentrischen Rohres unter Innendruck in ehler enthalten. Um auch den Fall einer 
Abhangigkeit vom Radienverhiltnis 9 und vom F . 
8 Mae ectahiaat auBer en Druckbelastung zu erhalten, braucht 
man lediglich einen hydrostatischen Druck p, 
zu iiberlagern. Die sich daraus fiir die Tangentialspannung innen und au®en ergebenden all- 
gemeineren Formeln sind von Weinel angegeben. 
Es ist zweckmaBig, dem Vorgang von Weinel folgend, an Stelle der Abmessungen To tae 
die dimensionslosen GréBen 


ut 


Care (RadienverhAltnis) und ier (ExzentrizitatsmaB) 
z 1 LO 
einzufiihren, so daB 
ao tO (Q. é) 
ee msnun =a (28) 


wird. Der Beiwert my ist in Abb.3 in Abhangigkeit von 9 und ¢ aufgetragen. 


Nach Untersuchungen an kaltgereckten, dickwandigen Rohren! wird die Beanspruchung 
am besten durch die Gestaltaénderungstheorie von Huber-Nadai wiedergegeben. In unserem 
Fall ist also fiir die gréBte Beanspruchung zu setzen 

es 7 


1 Bap oases a ; 


Onp=~, JV ik 
Tepe Aaa | “EH 


(29) 


1 H, Klein, Untersuchungen an kaltgereckten dickwandigen Rohren unter besonderer Beriicksichtigung 
der Veranderungen der Werkstoffeigenschaften. Diss. Aachen 1929. Mitt. KWI Eisenforschg. Diisseldorf 
Abh. 138 (1929), Bd. 11, Lief. 20. 


(Co} 79 — cos ¢) (Gin Hy + Tq 2 uv, cos -P). (26a) | 


—— a 
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3. Berechnung der Kriimmung. Fiir die axiale Spannung gilt nach (5) 


1 ES 
Oo —— ~ (Oy Oc) E(k Ja) 


m 


Hierbei wird Wy = ¢ Gtg 7) vom Mittelpunkt des Innenkreises aus gerechnet. Die Kon- 
stanten k und j sind so zu bestimmen, daB 


1. die itber die Querschnittsflache integrierten Axialspannungen gleich der vom Fliissigkeits- 
druck auf den Innenkreis ausgeiibten Kraft ist; 


2. das Moment der Axialspannungen bezogen auf den Mittelpunkt des Innenkreises ver- 
schwindet. 


Der rechnerische Ausdruck fiir diese Bedingungen lautet 


[ox oe == Te -o Dp. (30) 
F 
fo,ydF=0. (31) 
rs 


Die Spannungssumme s’ =o, +o; ist aber nach (17) oder (21) und (22) 


s' = 2 (24, + A;) +4 (Coj7 cos 6A, + Ein cos C A,) — 2 (Cof2ycos20A, | 39 
+ Gin 24 cos 2¢ A,). i323) 
Bei der Integration ist fiir das Flachenelement dF zu setzen: 
ce 2 
Die Durchfithrung der Integrationen ergibt (vgl. Anhang) 
[s'dF= is a 
F 
[s'ydF=0. 
F 
Damit gehen die Bedingungen (30) und (31) iiber in 
2 9 D) 9 9 = 9 
Sera tg 15) Peele to p, 
oy ee eel 
Eire ek— (5 (ham 0h eee ri] == {le 
J 
woraus man die Kriimmung 7 berechnet zu 
sony 2° 4 Aerory p 
1a (1 ae 7) it) Ae rer ) 


/ 


Ein Vergleich mit der in der Einfithrung auf elementare Weise abgeleiteten Formel (4) ergibt 
das einfache Ergebnis 


j=(1— si. (36) 


Es sei an dieser Stelle bemerkt, daB die von M. Bricas und L. Féppl benutzten Lésungen (14) 
zum gleichen Ergebnis fithren. Bei der Integration heben sich also die Unterschiede zwischen 
beiden Ansatzen gerade wieder heraus. 

Es empfiehlt sich auch hier wieder zu dimensionslosen GréfSen iiberzugehen. Mit dem Radien- 
verhaltnis 9 und dem ExzentrizitatsmaB ¢ erhalt man 


aa | ‘Ie (1 — 0) fe Bal owe all 
(l—e) 4g) . 

Es ist besonders interessant, die auf diese Weise fiir K erhaltenen Werte in Abhangigkeit von 
der Beanspruchung o, nach (29) aufzutragen. Abb. 4 zeigt, da sich die fir verschiedene ¢@ 
berechneten Kurven durch eine Hiillkurve einschlieBen lassen, welche fiir eine gegebene Bean- 
spruchung Héchstwerte fiir den Federweg liefert. Man erkennt daraus, daf® das Exzentrizitats- 
ma zweckmaBigerweise ¢—0,6 zu wahlen ist. In Abb. 4 ist auch angegeben, wie das Radien- 
verhaltnis zu wahlen ist, um bei gegebener Beanspruchung Bestwerte der Kriimmung zu erzielen. 
2 
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Bei der Anwendung dieser Ergebnisse auf die Praxis ist jedoch zu beachten, da man die | 
dickwandigen MeBrohre fiir sehr hohe Driicke einer Autofrettage (Kaltreckung) unterzieht, wobei | 


die inneren Zonen des Rohres zum plastischen FlieBen gebracht werden. Hierdurch werden dem 
MeBrohr nach der Entlastung innere Spannungen eingepragt, die der spdteren Betriebsbean- 
spruchung entgegenwirken. Gleichzeitig tritt hierbei eine Erhéhung der Streckgrenze des Werk- 
stoffes ein. Wenn man das exzentrische Rohr als DruckmeBglied verwenden will, mu die Be- 
triebsbeanspruchung im elastischen Gebiet bleiben, so daB die fiir die Kriimmung abgeleitete 


Formel auch bei vorheriger Autofrettage weiterhin giiltig bleibt. Jedoch ist es bei sehr dick- 
wandigen Federn miéglich, viel héhere 


K(@,€) Druckwerte zu erreichen, als es der nach © 
¢ . + 5p 5 Formel (29) errechneten Beanspruchung 
sae marleiarlk = + und der gemessenen Streckgrenze des 
4 Werkstoffes entspricht. Im auBersten Fall 
wird man etwa eine Verdopplung der 
Druckgrenze erreichen kénnen. Obwohl 
1 Versuche iitber 8000 kg/cm? bisher nicht 
bos vorliegen, wird man doch erwarten kénnen, 
We daB& bei geeigneter Autofrettage diese 
O4¥ Grenze noch erheblich iiberschritten wer- 
+ den kann. 
ge 4, Federweg des exzentrischen Rohres. 
Die Winkeldrehung des Federendes sowie 
ar der Federweg ergibt sich aus der Kriim- 
008 o=02 = mung in einfacher Weise. Die Winkel- 
086 if : P ung Arehung ist nadmlich in Bogenmaf 
Beanspruchung 
OOH i _» Druck ea Saetl 2 p 
3 0 Ag=lj mae m) KE . (38) 


Den Federweg f berechnen wir aus Abb.5 zu 


08 1 i 
f= R(1—cos Aq) = — RAG?=— Pj 
06 2 2 
oder 
i P 2 
fee oT eae’ Le 
02 [> 3 ( =) % Be (39) 
\Beanspruchung 1 
a eae TST TRS aT TE DE 20 30 Druck wenn d, = 2 r, der AuBendurchmesser ist. 


Um - 
Abb. 4. Kriimmung eines exzentrischen Rohres unter Innendruck in : erope Federweg ee erhalten, pore 
Abhingigkeit von der gréBten Beanspruchung. Die Bestwerte sind fiehlt sich also, die Querschnittsabmessun- 
durch eine Hiillkurve gegeben. 0: ¢ = 0,2; +:e=0,4; A:e=0,6; en so . hal 1 di 
@®:«=0,8. Im unteren Schaubild ist jeweils angegeben, wie das 8 geting Zu Halten, wie es die Sos 
Radienverhaltnis 9 zweckmaBigerweise zu wihlen ist. forderte Richtkraft und die Fertigungs- 
méglichkeiten erlauben. 

Meist wird der Federweg aber nicht am theoretischen Federende abgenommen, sondern an 
einem mit dem Federende fest verbundenen Teil. Abb. 6 zeigt zwei charakteristische Ausfiihrungs- 
arten. Im ersten Fall wird der Federweg an einem Verlangerungsstiick abgenommen, im zweiten 
Fall wird er zur Einspannstelle zuriickgefiihrt. Die Riickfithrung des Federweges durch ein 
itbergestiilptes Rohr’ bringt meBtechnische Vorteile. AuSerdem kann dieses Rohr gleichzeitig 
als Schutzrohr bei einem etwaigen Bruch der Hochdruckfeder verwendet werden. 


Bei einer Verlangerung nach oben erhalt man (vgl. Abb. 5): 


a 
fh=s(l +2 4) (40a) 
bei einer Riickfithrung nach unten dagegen 


fr=s(i— 29). (40b) 


1 DRP. ang. Sch. 119 817/42d. 


XIX. Band 1951. Wuest: Theorie des Hochdruckmefrohres mit ausmittiger Bohrung. 19 


5. Vergleich mit Messungen. Im Februar 1939 wurden von den Herren 0. Schumacher und 
H. Bremer auf Anregung der Firma F. Uhde, Dortmund in der MeBgerateabteilung der Firma 
Schaffer & Budenberg, Magdeburg-Buckau mehrere Messungen an exzentrisch gebohrten Hoch- 
druckmefBrohren durchgefiihrt. Bei der Auswertung dieser Messungen wollen wir so vorgehen, 
da wir zunachst aus den Formeln (40a) und (40b) auf den reduzierten Federweg f zuriick- 
rechnen und dann aus (38) j bestimmen. Unter Zugrundelegung eines Elastizitatsmoduls 
E = 2,1- 10% kg/cm? berechnen wir dann j, aus Formel (4). Aus dem Verhiiltnis j/j, = 1 — 2/m 
kann dann die Querkontraktionszahl m bestimmt werden und mit bekannten Werten verglichen 
werden. 


Messungen an exzentrisch gebohrten Rohren. 


ah Pp Ete al ele ee l a firs | joe jn 108 4 
kg/em? mm mm mm | mm mm mm 1/mm 1/mm A) 
| 

1 Mod pase gs i 300 85 4,08 58,0 121,5 0,478 
1 2000 4,5 6.5 1 300 85 4,26 60,4 121,5 0,498 
3 4000* 250 7,25 | 1 300 85 0,77 10,9 20.0 0,545 
4, 4000* 2,0 7.8 300 85 0,30 4,25 8,89 0,478 
5 2000 4.5 7,0 Hal 424 —A15 —3,45 36,7 ALI 0,516 
6 2000 Be We TTS 1 424, —415 —1,76 18,75 37.9 0,495 


* Diese Federn sind fiir einen Nenndruck von 8000 kg/em*® bestimmt. Aus mefstechnischen Griinden 


erfolgte die Messung bei 4000 kg/cm?. 


Man erhalt als Mittelwert fiir das Verhaltnis j/j) = 1 — 2/mden Wert 0,5. Diesem entspricht 
eine Querkontraktionszahl m=4, wie sie oft bei Festigkeitsrechnungen zugrunde gelegt wird. 


6. Zusammenfassung. Ausmittig ge- 
bohrte Rohre haben die Eigenschaft, dai 
sie sich unter der Wirkung eines Innen- 
druckes nach der starkeren Seite der Wand 
durchkriimmen. Diese Eigenschaft kann 
meBtechnisch insbesondere zur Messung 
hoher Driicke ausgenutzt werden. Der 


y 
y 
g 
g 


SS 


Ss SS SOO MMA Mi? WS 
Abb. 5. Berechnung des Federweges bei einem Ver- Abb. 6. Schematische Darstellung zweier Ausfiihrungsbeispiele 
langerungsstiick oder bei Rickfiihrung der Bewegung. yon HochdruckmefSrohren mit Ablesung des Federweges 


durch eine MeBuhr. 


Spannungszustand in exzentrischen Rohren ist von verschiedener Seite mit z. T, voneinander 


abweichenden Ergebnissen behandelt worden. Die Beziehungen fiir den Spannungszustand 


werden daher erneut abgeleitet und auf die Berechnung der Kriimmung angewandt. 


Anhang. 


Zur Berechnung der Integrale (30) und (31). Es ist niitzlich, zunachst einige Umlauf- 


_integrale zusammenzustellen, die fiir die Berechnung dieser Integrale gebraucht werden. Da 


solche Umlaufintegrale bei Rechnungen mit bipolaren Koordinaten vielleicht auch anderswo 
_ gebraucht werden, ist diese Liste noch etwas vervollstandigt worden. 


Die 
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iY dé ayn A = ; Al) | 
otis eoicny Sin 7 7 ( | 
0 ! 
i ei? Co} AQ 
| ay eat = a ieanae (A2) 


PAgA 


Col? 7 1 1 wok 
7a). = 27 3 oa I A3 
i (Cof n- — cos fe m( Sin 1) 1 2 Gin’ 7 ( ) | 
2a ue ae 
a (= ae een (A4) | 
| — cos C)* Sin? 7 2 Sin’ 7 
ra cos € dt cosC dt y = ast 
fe ere = 2% (Oto — 1), (45), 
€ cos € dt teh 1 A6 
ic oj 7 — cos C)? ere (A6) | 
2a 5 . 
; cos Cd oa Sof ae 
(oj 7 —cosl)> ONT j (A7) 


(Cvj 7 —- cos £)” 


2n 
| sin € d¢ 
Qn 


ik woe = 2m (btg? 7 — 3Ctg 7 +2), (A9) 
| Gin=esor = 8 Sara” (A10) 
ee 

erp auc (AL) 

9 : 
| (Go 9 a oe => (A12) 
0 


Die einzelnen Glieder der Integrale (30) und (31) nehmen damit folgende Form an: 


ny 2x 

| | (Go) =a So oe a( Sit h af Al , (A13) 
ny 0 

a Co neosCdyndg il 1 

| | (Gof 4 — cos ¢)? fer nm Sie H : (A14) 
m 

a in 7 cos € dy dé 5 ; ito 

| | (Cof 1 — cos ¢)* x (Gtg 7, — Cty r79). (A15) 


ny 0 
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"0 2a 
Cof 2 7 cos 26 dy dé , 
| | (Co} 7 — cos £)2 ~ 2x (Co} 2 of ca Gj : No — Sin 2 Hi + Sin 2 No) 


} 
m 0 | l l | (Al6) 
2 Sin? m 2 Sine an)? | 


"0 20 


~ ( Sin 2 9 cos 26 dy de 
| | (Co} " — cos C)? 


2 7% «(Cot 2 Iho — Goj 2 41) — (Gin 2 49 — Sin 2 77,), | (A17) 
[ 
) 


a +2 (Gtg 71 — Ctg ro)]. 

yo 27 

[Sing dyde (Gof, __ Gof 

| | (Cof 7 — cos 6)8 = | Sin? 1, Sin’ | “ (A18) 
2 7 

Ny 2a 

“ ( Gin y Coj 9 cos € dy dé F ; 

| | (Goin —eos B27 Wtg? 7, — Cty? 9). (A19) 
"1 0 

No 2% 

| “Gin? ycosSdy do 3 1 1 
I 3 (Coj n —cos fs — 27 ( Sint, aria a (A20) 
iy 0 


No 27 
| | Sin 7) Co} 2 7 cos 26 dy dé 


(Coj 7 — cos 6) ~ —= mt (tg? 4, — Ctg? n) +3 Gtg 7,— 3 tg), (A21) 


"4 0 
10 2 
Sin n Sin 2 cos 26 dy dé ] : 
=3 Te : 
| | (Go} 7 — cos ¢)8 fers Ny Sin? “al (A22) 
n, 0 


Insgesamt wird also das Integral (30) 


[s'aF = 2 (A+ Ay) x (r} — v3) — 4c? (A, — Ay) [Gof 2, — Gin 24, — (Gof 2 yy — Sin 2 np) 


2 


r? —r2 Tei . : 
i onwatuniieualt ake 
und das Integral (31) 
1 é 
[sydF = 2A A{ Sine y, O18 ie one Ctg Tela 2s (Gtg? 7, —Ctg? ny 
F 


— 3 tg 7, +3 Ctg no) mc? — c&tg ny [sd F 
EF 


=[er? + (ri — 15) eGtg mo] A Ur} 15 +1911 Go} (jo — m1)] et As 
—2cign uri p= 0. 


(Eingegangen am 10. Juni 1950.) 
Anschrift des Verfassers: Dr. W. Wuest, (20b) Gottingen (Hann.), Zeppelinstr. 7. 
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Beitriige zur Theorie der durch gleichférmigen Schub 
beanspruchten Platte. 


Ill. Mitteilung. 
Von W. Kucharski.. 


1. Die grundlegende Differentialgleichung als Wellengleichung. In der folgenden Mitteilung 
soll unter Hinweis auf die Grundlagen der in Bd. 18, Heft 6 dieser Zeitschrift zum gleichen 
Thema verdffentlichten Mitteilung I die grundlegende Differentialgleichung 
vom Standpunkt der Wellenlehre behandelt werden'. Durch 


w= f(xcos & + y sin a — ct) (2) 
wird eine Welle mit geradliniger Front dargestellt, die normal zu sich selbst mit der Wellen- 
geschwindigkeit c fortschreitet. ¢ bildet mit der x-Achse den Winkel a. Der Sinn dieser Aussage 
ist der, daB w fiir konstante Werte des Argumentes 


S;= xcosa-+y sina — ct (3) 
konstant ist. s, = konst. bedeutet aber 
ct 
y= — ctg ax iar ++ konst. (4) 


Dies sind fiir jeweils feste Werte von t gerade Linien, die zueinander parallel sind und mit 


der x-Achse den Winkel 
raf : 


einschlieBen (Abb. 2). Bei wachsendem t wird offenbar der gleiche Wert von s und damit von w 
fiir entsprechend vergréBerte Werte von y und x erreicht; diese liegen auf einer zu der gezeichneten 
parallelen Geraden, die fiir t-+-At in der Normalrichtung um cA t fortgeschritten ist. 
Fir 
Sg—= * C08 & + ¥ sin & -- ct (6) 


gilt Entsprechendes, mit dem einzigen Unterschied, dai die Richtung von ¢ gegeniiber der in 
Abb. 2 eingetragenen umgekehrt ist. 


Wenn Striche Differentiationen von f(s) nach dem Argument s bedeuten, so wird 


, Ib 
Wis —— Cf © eae Ziti Wry = sin a cos af" = | sin 2 af”; 


setzt man dies in (1) ein, so erhalt man 
weet 1 - Mt 
f silties 59 sin 2a)f aii: (7) 


2. Fortschreitende Wellen. Der Schubspannungswert @ wird hier als konstant an- 
genommen: 


& = konst. (8) 


Damit stellt dann (7) eine gewdhnliche lineare homogene Differentialgleichung fiir f dar, mit 
der allgemeinen Lésung 


f=A+Bs+C, sin Ks+C, cos Ks (9) 
mit A, B, C,, C, als Integrationskonstanten und 


1 
Ke 2 eas 
c 5 Osin2a. (10) 


1 In dieser Mitteilung wird It. Abb. 1 in Ubereinstimmung mit der sonstigen Literatur tT entgegen- 
gesetzt zu der Festsetzung in den Mitteilungen I und II als positiv angenommen. 
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Wenn K? >0, so sind C, und C, in (9) reell, und sin bzw. cos haben reelles Argument. Fiir 
K?< 0 treten die Hyperbelfunktionen mit reellem Argument an die Stelle der Kreisfunktionen, 
und C,, C, sind wieder reell. 

Die Hyperbelfunktionen und ebenso das lineare Glied in (9) sind hier von sekundarem Imteresse, 


da sie standig wachsende Betrage von w ergeben. Daher wird die Untersuchung im wesentlichen 
beschrankt auf 


f= C,sin Ks +-C,cos Ks , Ke 0m (11) 
Wie aus (9) hervorgeht, hat die Platte (ebenso wie der elastische Stab) die Eigenschaft, daB 
nur ganz bestimmte.einzelne Wellenformen méglich sind, die (abgesehen von dem 
linearen Glied) durch den betreffenden Wert von K gekennzeichnet sind, und da® die Fort- 
pflanzungsgeschwindigkeit c fiir die verschiedenen Formen verschieden ist. Betrachtet man 
in (11) der Einfachheit wegen das Sinusglied, so erkennt 

7] man leicht, dai 


20 2% 


K=- r — (12) 


M=Sin #S 


Abb. 1. Hauptbezeichnungen. Abb. 2. Wellenfront mit Angabe der Winkel. Abb. 3. Eine einfache stehende Welle. 


gilt, worin | die ,,Wellenlinge“ fiir eine Vollperiode bedeutet (Abb. 3). 
Hiermit ergibt (10) den Zusammenhang zwischen der Wellengeschwindigkeit und 
der Wellenlange durch die Gleichung 


2 
oa iT 1 2 sin2a. (13) 


Also ist ¢ von | abhangig; die Wellen haben Dispersion; d.h. eine Wellenform, die zur Zeit t 
auf einem gewissen Gebiet der Platte durch Ubereinanderlagerung endlich oder unendlich vieler 
Partialwellen von verschiedenem K, d.h. mit verschiedener Wellenlange, entstanden ist, kann 
nicht bestehen bleiben, da die verschiedenen Anteile mit verschiedenen Geschwindigkeiten 
weiterschreiten. 

Die Verhaltnisse liegen bei der Platte (und beim Stab) vollkommen anders als bei der stark- 
gespannten Saite. Bei dieser sind beliebige Formen méglich, und cist unabhangig von der Lange 
der harmonischen Partialwellen. 

In dem hier untersuchten Fall der durch Schub vorbeanspruchten Platte tritt 
nun nach (13) ein dritter Unterschied gegeniiber der Saite und ein Unterschied gegeniiber 
dem elastischen Stab und der schubfreien Platte auf: 

Die Wellengeschwindigkeit ist auBerdem von der Richtung abhangig, um 
so starker, je gréBer bei gleichem K der Schubspannungswert @ ist. 

Die Verhaltnisse werden anschaulich, wenn man beachtet, daf der Spannungszustand mit 
Ty = Ty, = 0in Schnitten parallel zu den Koordinatenachsen fiir die dazu unter 45° liegenden 
Schnittrichtungen Zug- bzw. Druckspannungen vom gleichen Betrage liefert, wie in Abb. | 
durch die Pfeile angedeutet. Die Druckspannungen wirken im Sinne einer Vergréferung, die 
Zugspannungen im Sinne einer Verkleinerung der Auslenkungen w. Hiermit steht die Vergrofe- 
rung von cin der Zug-, die Verkleinerung in der Druckrichtung unmittelbar in Zusammenhang. 

Offenbar ist stets bei Betrachtung von c (a) 


1 


Cue a= Ke + 5 On Cmin == K?— 9 a, (14) 
und fiir jedes K? gibt es einen kritischen Wert von 0 
Oy == 2 K?, (15) 


fiir welchen 
Cnin = Oe 


Dies ist offenbar derjenige Wert, bei welchem die Platte in Parallelstreifen unter 45° mit der 
Wellenfront normal zur Druckrichtung ausbeult, 


Man erhalt. eine einfache Darstellung der Verhaltnisse, wenn man z. B. die kritische | 


Wellenlange [x einfiihrt, fiir welche jeweils Cin = 9 wird, die Platte also gerade an der Aus- 
beulungsgrenze sich befindet. Unter Beriicksichtigung von K—2q/l erhalt man aus (14) 
finnce. 0) 


fee 


Hiermit erhalt (13) fiir die Wellengeschwindigkeit die Form 
c'2 = (Ix/l)? + sin 2 0, 


wenn man mit 


eine auf die jeweils mabgebende konstante # bzw. lx bezogene GroBe einfiihrt. Diese Darstellung 
hat den Nachteil, daB c’ mit 0? — 0 nach | 


Unendlich geht. 


Fiir eine anschauliche Kurvendarstel- 
lung ist der Wert c- 1 geeignet; man erhalt 


42 =T = 1 +(7,) sin 2a 


Dieser Wert ist in Abb. 4 mit y als Radius 
vom Nullpunkt und « als Winkel eines 
Polar- Koordinatensystems aufgetragen. 
7 Fir 1= 0 ergibt sich der Kreis um den 
Mittelpunkt mit dem Radius Eins; fiir den 
maximal zuladssigen Wert | =1[x erhalt man 


zwei Kreise vom Radius y1/2, die sich im 
Zentrum beriihren und deren Mittelpunkte 
auf der mit 45° gegen die x-Achse im posi- 
tiven Sinne geneigten Geraden liegen. 
y (a) ist diejenige Funktion, die beim An- 
wachsen des Wertes 9 von Null bis zum 
jeweils zulassigen Maximalwert fiir (I/lx) 
= 1 die zuletzt genannten Kreise in den 
einen zentralen Kreis von gleichem Flachen- 
inhalt iiberfiihrt. 


Abb. 4, Die Abhangigkeit der Wellengeschwindigkeit yom Zentriwinkel. 


3. Stehende Wellen in Parallelstreifen mit staindig ruhenden Randern. Diese sind leicht herzu- 
stellen durch den Ansatz 


w= 5 sin K (h— ct) 4-5 sin K(h +61), (16) 


worin zur Abkiirzung 
h= xcosa +y sina 
gesetzt ist. In bekannter Weise lat sich dies umformen zu 


2ah Qirnie 


w= sin Khcos Kct = 7 COs ihe (17) 


das bedeutet eine stehende Schwingung in der durch sin Kh bestimmten Form mit der 
Kreisfrequenz 


a rc 
o= Kes (18) 
Offenbar ist 

w—0 fir alle t 


fiir sy nesses lca (k = 0, 1, 2, 3) 
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oder fiir 
l 


(= —ctgax tk. 
J 8 _” 2 sina 


(20) 


Dies sind fiir die aufeinanderfolgenden Werte von k jeweils Gerade, deren Normale mit der 
x-Achse den Winkel « einschlieBt, deren Abstand 1/2 ist und die auf der y-Achse Strecken von 
der Linge 1/2 sina abschneiden. In Abb. 5 ist dies angedeutet; in ihrem rechten Teil ist w fiir 
einen Wert von t schematisch dargestellt durch Umklappung des Diagramms in die Zeichenebene. 


S y 
Af 
as 
Wee / 
sw RY 
he tite: 
<< go » 
Ve: ésinoc\| / 
x POEL SS 
he 
aS pana 
— vA 
™ 
eS 
/ S 


Abb. 5. Stehende Wellen in unendlich langen 
Parallelstreifen. 


Abb. 6. Plattenstreiten mit schriger Randbelastung. 

Je zwei aufeinanderfolgende Gerade wo schlieBen einen unendlich langen Parallel- 
streifen von der Breite 1/2 ein, der herausgeschnitten werden kann und in unveridnderter 
Weise schwingt, wenn er an den Randern entsprechend den Navierschen Bedingungen gelagert 
ist und wenn an seinen Randern in 
der Plattenebene geeignete Krafte 
angebracht werden. 

Diese ergeben sich durch eine ein- 
fache Benutzung der Formeln fiir den 
Spannungstensor derart, daBbihr Betrag 
gleich # ist und ihre Richtung mit der 
x-Achse den Winkel -f bildet. In Abb. 6 
und 7 sind zwei derartige Falle an- 
gedeutet. 

Naturgema kénnen in dem Ansatz 
(16) oder in der daraus hervorgehenden 
Form (17) auch die Sinus- anstatt der 
Kosinusfunktionen oder umgekehrt ge- 
wahlt werden (unabhangig voneinander), ohne daB im Ergebnis etwas Wesentliches geandert 
wird, da hierdurch die Parallelstreifen lediglich parallel zu sich gegeniiber den Koordinaten- 
achsen verschoben werden. 


Abb. 7. Plattenstreifen mit schrager Randbelastung. 


4, Stehende Wellen in gewissen Scharen von Parallelogrammen. Diese erhalt man durch 
Ubereinanderlagerung von zwei stehenden Wellen nach (17), deren Normalen mit der x-Achse 
die Winkel «, und a, bilden mégen. Mit 


h, = x cos 0, --y sin 41, 


) 
hi, == 4 C08 0, 4 Si. 5, | 
i eee = kee ee, { (21) 
I, L | 
si 20 : Se RELY 
Ue i A Ws = Pie ss | 


sel angesetzt 


w= w,+wW, 


mit 
— J it t == ey eye S -h, cos Wot (22) 
w= — = 00s K,h,cos@t, UW,= 5 ples ate 
Dies kann eine stehende Welle nur dann werden, wenn die zeitlichen beiden Frequenzen 


gleich sind; die zu wahlenden Daten unterliegen also, wenn man stehende Wellen fordert, der 
Beschrankung 


2 2 
K,¢c, = ae: =0,= K,c¢, = 7, % =o? =w | 
| (23) 


oder 


Die zweite Form dieser Gleichung ist unmittelbar anschaulich: Die Zeiten miissen gleich 
sein, in welchen jede der Welle um ihre Lange vorriickt. Hiermit entsteht aus (22) 


w= — S (cos K,h, — cos K,h,) cos wt, 


was sich wieder umformen laBt zu 


w= sin (k,h, + K,h,) sin (K,h,— K,h,) cosot. (24) 


Um in diesem Ergebnis die leicht itbhersehbaren sin-Funktionen zu erhalten, waren vorher in (22) | 


die Kosinus mit den entsprechenden Vorzeichen gewahlt. 
Aus (24) ergibt sich dann sofort die Umrandung der gesuchten Parallelogramme mit 


w= 0 
fiir 
K,h, + K,h,= tha (e002 eee) (25) 
und fiir 
K,h,— K,V=+hkh,a (eg 10 sa 1257) 


oder ausgeschrieben in x und y 


x (K, cos «, + K, cos «,) + y (K, sin «, + K, sin «,) = k, 2 


und 
x (K, cos a,— K, cos «,) + y (K, sin «,— K, sin a.) = kaa 

oder 

5: K, cos «, + K, cos a It i 

=— ~ ; : x = = : 

: K, sina, + K,sina,° ~— K,sino,+ K,sina, 1 | 

und (26) 
K, cos «, — K, cos a a k | 
x K, sino, — K, sina, ~~ K,sina,—K,sina, * 


Dies sind zwei Scharen von Parallelen, deren Winkel 6, bzw. 8, mit der x-Achse sich aus | 
dem betreffenden Faktor von x ergeben und deren Abstande aus den Abschnitten auf der y-Achse | 


folgen, die aus dem jeweils letzten Gliede in (26) abzulesen sind. 
Durch einfache Rechnung erhalt man fiir den Winkel 6, — 6, den die Geraden der einen 


Schar mit denjenigen der anderen einschlieBen, unter Einfiihrung der Wellenlangen /, und I,: | 


eo dade 


tg (By — Bi) = B—E sin (%:— 0). (27) | 


In (26) und (27) diirfen die vorkommenden GréBen nicht unabhangig vonein- , 


ander gewahlt werden; es besteht die Frequenzbedingung (23), in welcher c, und ¢, 


nach Gl]. (13) durch 1, bzw. 1, ersetzt werden kénnen. Man erhalt so die Gleichung zwischen — 


lhe orp ee van 


1 Uganiaions ib ol #sin 2a, 1 
Ene Tg yet a 28) 


wobei dann 
= 7 


+5 ap Silt 2/0 5 | 
und _* (29) 
is thee @ sin 2 X. J 
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Die Lésungen sind nur periodisch fiir reelle c, also fiir 


c? = 0) sund teal UBS (30) 
Man kann dies mit (28) in einer iibersichtlichen Form zusammenfassen, etwa durch 
il ff @ sin 2 X 7. 1 # sin 2 x 
lee pe) \ eed eee 9 
i: (1 8 m2 |= li (1+ mou. : (28a) 


Jetzt kann man beurteilen, fiir welche Parallelogramme die Lésung der 
Schwingungsgleichung (1) durch den einfachen Produktenansatz (24) gegeben 


eke, wird. In diesen sind offenbar die Wellen- 

; fronten die Diagonalen, 6,— f, der eine, 
x — (h,— B,) der andere Eckenwinkel 
und «,—«, der eine der Diagonalen- 
winkel. In Abb. 8 ist dies schematisch 
angegeben; hierin sind die Fronten der 
erzeugenden Wellen gestrichelt, die ent- 
stehenden Parallelogramme stark aus- 
gezogen. Die Parallelogrammseiten ver- 


SS 


Abb. 8. Parallelogramme fiir stehende Wellen. Abb. 9. Quadrat unter Zug und Druck. 


binden Schnittpunkte der im Abstande 1,/2 bzw. 1,/2 gezogenen Wellenfronten fiir w— 0. In 
der Abbildung ist auf (28) noch keine Riicksicht genommen. 

Es seien nun einige interessant erscheinende Fragen aufgeworfen. Gehéren Rechtecke zu 
diesen Parallelogrammen? Hierfiir miiBte 6, — 6, = 2/2, tg (Bb, — B,) = © sein. 

Nach (27) ist dies nur méglich fiir 13 — 1? ; Gleichung 
(28) liefert dann bei c? 4 0 

sin 2 4, = sin 2 a. , 

Bei der Deutung dieser einfachen Gleichung ist zu be- 
achten, daB die Nullstelle des Sinus sich bereits beim 
Weiterschreiben des Argumentes um 7 wiederholt, wah- 
rend fiir sin2« ~0 die Gleichung 20,=20,+k2z, 
0, = 0) +a folgt. 

Die zweite Alternative liefert o,—a,—=+27; 
die Diagonalen der Rechtecke arten also entweder in 
unendlich lange Parallelstreifen aus oder in zwei 
aufeinanderliegende Strecken von begrenzter Lange, also 
Rechtecke von der Breite Null. ‘sti rinnung 

Dagegen liefert sin 2 «, = sin 24,= 0, z.B. «, = 0 cr a Deedee ip yetee teen os . 
und a,—27/2. Die Diagonalen bilden dann ebenfalls 
rechte Winkel; die Rechtecke sind Quadrate. Dabei sind die erzeugenden Wellenfronten 
den Koordinatenachsen parallel, die Quadratseiten liegen also zu diesen unter 45°. Es 
handelt sich um das in Abb. 9 dargestellte Quadrat. 

Andere Rechtecke gibt es bei c2?40 mit dem einfachen Ansatz (24) nicht. Auf den Fall 
c? = 0 wird weiter unten eingegangen. 

Gibt es schiefwinklige Rhomben ? In diesen stehen die Diagonalen senkrecht aufeinander, 
es ist also sin (x) — &,) = 1. Dies ist in vielfacher Weise zu realisieren; schiefwinklige Rhomben 
sind méglich. 
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Fiir einen allge meinen Uberblick tiber die in der angegebenen Weise erhaltlichen Parallelo- 
gramme ist von (28) auszugehen. Dies besagt, daB sich solche Wellenziige zu stehenden Schwin- 


gungen superponieren lassen, bei denen der Wert 
1 0 


See ene ese FS (31) 


Ge Ihe 
einen fiir die beiden Wellenziige konstanten Wert hat. 

Bei gegebenem Wert von C ist hierdurch eine Kurve /(a) definiert, auf welcher 
die Endpunkte der vom Nullpunkt gezogenen Strecken 1, und J, liegen miissen. Normal zu diesen 
Strecken liegen die Wellenfronten der erzeugenden Wellen, d.h. die Diagonalen der gesuchten 


Parallelogramme. In Abb. 10 sind derartige Kurven fiir C = 1 und verschiedene Werte von | 
) angegeben. Abb. 11 zeigt eine in dieser Weise erhaltene Rhombusform, Abb. 12 ein Parallelo- | 


gramm mit zwei zur y-Achse parallelen Seiten. Es ist nicht méglich, in der angegebenen | 


cinfachen Weise Rechtecke zu 
erhalten, bei denen die Seiten 
jeweils den Koordinatenachsen 
parallel sind. 


Abb. 12. Ein Parallelogramm mit zur 
Abb. 11. Ein Rhombus mit schrager Randbelastung. y-Achse parallelen Seiten. 


Gleichung (31) liefert reelle Figuren auch noch fir C—0. Die betreffenden Formen 
stellen dann Platten ‘dar, die sich bei den entsprechenden Randbelastungen gerade an der 
Grenze der Ausbeulung befinden; die Schubspannung hat fiir sie gerade den kritischen 
Wert. Gleichung (31) ergibt mit C= 0 nach Multiplikation mit [4 


Sar il 
Danae es e 
Osim 2x (32) 
Dies hat reelle Lésungen fiir 1 nur fiir 
ina oo =< We 


Die betreffende Kurve zerfallt in zwei hyperbelahnliche Zweige, die im zweiten und vierten 
Quadranten liegen. In Abb.13 ist diese Kurve fiir ? = 827? dargestellt; die Abb. 14 und 15 
zeigen Parallelogramme, die in dieser Weise erhalten wurden. 

Wie aus Abb. 15 hervorgeht, gehéren hierzu auch Rechtecke; hier liegt der Sonderfall c2? = 0 
vor, der friither bei der Frage der méglichen Rechtecke erwahnt wurde und fiir die Gleichung (31) 
noch andere Lésungen als Quadrate besitzt. Die jetzt erhaltenen Rechtecke liegen ebenfalls 
unter 45°, lassen sich aber nicht in Quadrate iiberfiihren; ihr Seitenverhaltnis nahert sich bei 


unbegrenzt zunehmender Seitenlange beliebig nahe der Einheit, ohne diesen Wert exakt zu 
erreichen. 


; : aa we ee one 
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Die der Druckspannung parallele Seite liBt sich beliebig klein machen, wie man aus der 
Konstruktion ohne weiteres erkennt. Es ist aber nicht moglich, daB diese Seite langer wird 
als die in Richtung der Zugspannung hiegende; in diesem Fall kann ein Ausbeulen moht in der 
vorausgesetzten Form eintreten. 

Fiir diese Rechtecke lat sich nun leicht die Formel fiir den zugehorigen kritischen Wert 
von ¥ aufstellen: Wenn man beriicksichtigt, da® 2 « der Winkel zwischen I und der Diagonalen 


4 


Abb. 13. Charakteristische Kurve zur Gewinnung von 
Plattenformen beim Ausbeulen. Abb. 14. Parallelogramm bei Ausbeulung durch schrage Randkrifte. 


ist, so kann man in (32) durch einfache Rechnung / und sin 2 « durch die Seitenlangen a und 6 
eliminieren, wobei a in Richtung der Zugspannung, 6 in Richtung der Druckspannung ange- 
nommen ist. Mit 6= b/a erhalt man so a 
8a? (1462) 
fa (10) 
Dies geht nach unendlich fiir 6 > 0 und 6—> 1. Man 
darf aber hieraus nicht den SchluB ziehen, daB die 
quadratische Platte unter der hier angenommenen Be- 
lastung iiberhaupt nicht ausbeult. Gleichung (33) und 
das Vorausgegangene bedeutet lediglich, daB die qua- 
dratische Platte unter dieser Belastung nicht als 
Ganzes, d.h. nicht ohne Knotenlinie, ausbeulen 
kann. Offenbar kann sie aber durch eine in Richtung 
der Zugspannung liegende Mittellinie in zwei Recht- 
ecke mit der einen Seite a und der anderen b = a/2 
unterteilt werden, fiir welches sich aus (33) mit d= 1/2 
der tatsachliche niedrigste Wert von #—0O fir 
das Quadrat unter gleichem Zug und Druck zu 

) Abb. 15. Rechteck bei Ausbeulung durch Zug 


Drrit = 3 ae \ 
t und Druck. 
ergibt. ‘ 
'  BeiderAnwendung dieserFormeln sind naturgem4B die hier benutzten Dimensionen zu beachten. 

Die héheren kritischen Werte erhalt man in einfacher Weise durch Unterteilung des Quadrats 
in mehr als zwei kongruente Rechtecke oder Quadrate. 

Fiir die anderen aus (32) entwickelten Parallelogramme mit einfachen Ausbeulungsformen 
lassen sich ebenfalls die zu (33) analogen Beziekungen aufstellen, die dann naturgemaB weniger 
einfach ausfallen. Da es sich um elementare Umrechnungen der einfachen Grundformel (32) 
handelt, kann die Durchfiithrung der Rechnung hier unterbleiben. 

Es sei noch vermerkt, da (33) fiir das unter 45° liegende Rechteck fiir beliebige, auch irratio- 


nale Seitenverhaltnisse 0 gilt. 


(33) 


a krit —— 


2 
200 a (34) 
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5. Ein allgemeiner Ansatz fiir stehende Wellen. Die entwickelten Parallelogramm- | 
formen gelten nur fiir den einen speziellen Wert von #, fiir den sie jeweils bestimmt wurden. Es 
ist nicht méglich, ihre Abmessungen und Winkel beizubehalten und # zu verandern. Eine 
gewisse Ausnahme hiervon machen nur die fiir ¢min = 0 angegebenen Formen, bei welchen die 
Winkel bei Anderung von @ erhalten bleiben und die Dimensionen sich so andern, da wieder 
die Ausbeulungsgrenze vorliegt. 

Die praktisch wichtige Aufgabe, die Lésung fiir eine bestimmte festbleibende Abbildung als 
Rand bei beliebig verinderten Werten von # zu finden, wird also durch die bisherigen Ergebnisse 
nicht gelést; die angewandte Methode der Ubereinanderlagerung von Wellen miiBte hierfitr — 
verallgemeinert und modifiziert werden. . 

Als erstes sei untersucht, was sich beim gleichzeitigen Auftreten ,,.unendlich vieler“ gerad- 
liniger Wellen mit stetig verandertem Winkel a ergibt. 

Diese sollen zu einer stehenden Schwingung fahren, ihre Wellenlangen und Winkel 
unterliegen also der Bedingung (28), welche die Gleichheit bei der stehenden Wellen auf- 
tretenden Schwingungsfrequenz ausdriickt. Da in der friiheren Mitteilung I gezeigt wurde, daB 
es Schwingungen der gesuchten Art mit harmonischer Zeitfunktion gibt, geniigt es, in dem 
Ausdruck (17) fiir eine geradlinige Einzelwelle die Koordinatenfunktion 


U(X, ¥,-0) asi - (xcosa + ysin a) (35) 


zu betrachten. Eine solche Welle habe die Intensitat A (a) da, der Allgemeinheit wegen sei sofort 
noch das Kosinusglied hinzugefiigt, so da man fiir eine allgemeine Ubereinanderlagerung 
zu dem Ausdruck 
2x 
U(x.) = [[4@) sin 77 
gelangt. Hierin kann nun bei der Integration / nicht, wie in vielen bekannten Problemen dieser 
Art, konstant gehalten werden; vielmehr ist nach (31) / eine Funktion von «a, fiir welche 
sich durch Auflésung von (31) nach | die Form ergibt 


1 fer a eae eas 
— === |} sin 24 Vs n22a4 k 
rE, | ysl 1 


(x- cosa + y sin a) + B(x) cos ae (xcosa+ysina)|da (36) 


| (37) 
orn et hae 6 
als neu auftretende Systemkonstante eingefiihrt ist. 
_ Der allgemeine Ansatz (36) geht hiermit iiber in 
2x / ‘ 
Us. 98) = f | Ae) sin [2 YR, eh ree _| 
eo | sin 20 + /sin? 20+ k, (38) 


da. 


+ B(a) cos (2212 k, j seen ) 
/sin 2a + ) sin? 2a + ki /_ 
Ob ein derartiger Ansatz mit den darin auftretenden Funktionen zu einem Erfolg fiihrt, mag 
hier dahingestellt bleiben. Nachdem erkannt worden ist, welche komplizierten Eigenschaften 
die schubbeanspruchte Platte in bezug auf die Fortpflanzung von Wellen besitzt, ist das 
Auftreten derartig schwieriger Formen wie (38) kaum verwunderlich; auch erklart sich durch 
diese Anschauungsweise von einem anderen Gesichtspunkt der verhaltnismakig groBe Aufwand, 
den die Untersuchungen in den Mitteilungen I und II erforderten; er liegt in der Natur des 
Problems. 
o Zusammenfassung. Als Ergebnisse dieser Mitteilung kénnen etwa folgende gebucht 
werden: 
a) Die Erkenntnis iiber die eigentiimliche Struktur der schubbeanspruchten Platte beziiglich 
der auf ihr méglichen Wellengebilde und ihrer Fortpflanzungsgeschwindigkeit. 
b) Die Gewinnung von Parallelogrammformen als einfache Liésungen fiir spezielle Werte der 
Schubbelastung und mit jeweils verschiedenen Randbelastungen. 
c) Die allgemeine Ausbeulungsformel fiir Rechtecke unter 45° und fiir andere Parallelo- 
gramme mit entsprechender Randbelastung. 
d) Ein Einblick in die grundsatzliche Schwierigkeit allgemeinerer Ansatze. 


(Eingegangen am 12. Juni 1950.) 
Anschrift des Verfassers: Professor W. Kucharski. Berlin-Charlottenburg, Hardenbergstr. 34. 
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Beitrag zur Berechnung der turbulenten Grenzschichten*. 
Von J. Rotta. 


1. Einleitung. Fiir die Beurteilung der Strémungsverhaltnisse um einen Kérper und ins- 
 besondere zur Abschatzung seines Strémungswiderstandes ist das Verhalten der an den Kérper 
angrenzenden Strémungsschicht, die entweder laminar oder turbulent sein kann, von sehr groBer 
_ Bedeutung. Wahrend bei den laminaren Grenzschichten die physikalischen Verhiltnisse geklart 
und auch die mathematischen Probleme soweit durchgearbeitet sind, daf fiir die Praxis brauch- 
| bare Rechenverfahren zur Verfiigung stehen, ist die exakte Berechnung der turbulenten Grenz- 
_ schichten noch nicht midglich. Die bekannten Naherungsverfahren! lassen sich verbessern, wenn 
| man beachtet, dai die kinematische Zahigkeit » und die geometrische Beschaffenheit der Wand 
) (Wandrauhigkeit) das Geschwindigkeitsprofil nur in einer gegen die Grenzschichtdicke 6 sehr 
diinnen Schicht 6,,in Wandnahe beeinflussen, und daf die StrémungsgréBen im iibrigen Bereich 
der Grenzschicht bei geeigneter Normierung als von der Zahigkeit und der Wandrauhigkeit fast 
unabhangig erscheinen. 
| 


2. Grenzschichtgleichungen. Die folgenden Ausfithrungen beschrinken sich auf im Mittel 
stationare Strémungen einer inkompressiblen Fliissigkeit. Die x-Achse liege parallel zur Wand, 
und y sei der senkrechte Abstand von ihr. Der zweidimensionalen Grundstrémung, welche durch 
die zeitlich gemittelten Geschwindigkeitskomponenten U und Vin x- und y-Richtung gebildet 
wird, iiberlagert sich die ungeordnete, zeitlich schwankende Turbulenzbewegung mit den Kom- 
» ponenten u,v und w, die stets dreidimensional ist. Die Integration der Grenzschichtgleichung 


fiir die ebene Wand 
U 0U pou Op OT 2] 
0 PPR: ae ee ray: 2} 


in welcher p den zeitlichen Mittelwert des statischen Druckes, 9 die Dichte des Strémungs- 
' mediums und 


GU ree : 
oS ae = uw) (2,2) 


die zeitlich gemittelte Schubspannung in der x y-Ebene sind, fithrt mit Benutzung der Kon- 
 tinuitatsgleichung 


0U OV 
a +5=0 (28 
| und den zugehérigen Randbedingungen 
as W=a05 ea OF aes v= w= 0, | (2,4) 
yo: U+U,; wd; vo, w2 > 0 | 
| auf die bekannte Impulsgleichung 
d ans dU ic 
ax Ui) aie Ufo, 7 : (2.5) 
- Hierbei sind 
Pe eee ay (2,6) 
1 U, 
q 0 ‘ 
die Verdrangungsdicke, 
| ea ele pees 22 


* Eine ausfihrlichere Darstellung findet sich in: J. Rotta, Uber die Theorie der turbulenten Grenz- 
| schichten. Mitteilungen aus dem Max-Planck-Institut fiir Strémungsforschung Nr. 1. Gottingen 1950. 

1 E, Gruschwitz, Ing.-Arch. 2 (1931), S. 321. — A. Kehl, Ing.-Arch. 13 (1943), S. 293. — A. E. v. Doen- 
_ hoff, N. Tetervin, Determination of General Relations for the Behavior of Turbulent Boundary Layers. 
NACA-Report No. 772 (1943). 
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die Impulsverlustdicke, U, die am AuSenrand vorgegebene Geschwindigkeit und t) die Wand- | 
schubsspannung. ae: : 
Der in (2,2) auftretende zeitliche Mittelwert des Produktes wv der zueinander senkrecht wir- | 
kenden Schwankungskomponenten ist als weitere, bei laminarer Strémung nicht auftretende 
Unbekannte anzusehen. Man muB sich nun nach weiteren Beziehungen umsehen, die den Zu- | 
sammenhang zwischen dieser Unbekannten uv und den iibrigen StrémungsgréBen vermitteln. | 
Kine solche Beziehung ist z. B. die Bilanz des zeitlichen Mittelwertes der auf die Masseneinheit — 
bezogenen kinetischen Turbulenzenergie 
u2 + v2 + w? 
E 7 is 
worin Uberstreichungen wieder zeitliche Mittelung bedeuten. Nach einer solchen Bilanz, die erst- 
malig von L. Prandtl! angegeben wurde, ist die in der Zeiteinheit erfolgende totale (substantielle) 
Anderung der Turbulenzenergie gleich der Differenz zwischen der der Grundstrémung ent- 
nommenen Energie (Arbeit der gemittelten Spannungen) und der durch Reibung in Warme um- 
gesetzten Energie (Dissipation S) einschlieBlich der mit benachbarten Orten ausgetauschten 
Energie (Energiediffusion Q). Mit den bekannten Grenzschichtvereinfachungen lautet diese | 


jae: (2,8) 


Energiebilanz 
dE 0E au aQ 
Oe ee SE yas “09 = 1 aaa (2,9) 
- Ox Oy Oy ¥ 
———————S=s —— — —_—. 
Totale Anderung Der Grund- Dissipation Energie- 
der Turbulenzenergie strémung ent- diffusion 


zogene Energie 


Hierbei ist die auf die Masseneinheit bezogene Dissipation 


dU \2 du\2 dv 
Ge HF 0 
und die Energiediffusion in Richtung der y-Achse 


Lig eh are aes 
lech 


2 dw 2 dw , dv\2 du Ow\2 
: [SS ee eco Lagi fes 
+2(32) oe ! a2} eee 


ae aaa 
oe nig ay) 


S=y (2,10) 


2 Puss Ape ee) SEALY 
aia al & insed mat (2,11) 
2 Q 
Mit p werden die zeitlichen Schwankungen des Druckes bezeichnet. Durch Integration iiber y 
1aBt sich hieraus ein Energie-Integralsatz entwickeln, der die Form 


co 


freaky d [ 
> dx (Ui5s) = D ee | UE dy (2,12) 
“Bnexgiestromvertast” Tnaiipation = : —~ 


Zunahme des 
Turbulenzenergiestromes 


der Grundstrémung 


hat mit der Energieverlustdicke 


(2,13) 


und der Dissipationsfunktion 
(2,14) 


Wahrend die Impulsgleichung (2,5) fiir turbulente Grenzschichten die gleiche Form wie fiir 
laminare hat, ist der Energiesatz (2,12) fiir das Verhalten turbulenter Grenzschichten bezeich- | 
nend. Die Energieverluste der Grundstrémung werden im wesentlichen zunichst in kinetische 
Turbulenzenergie umgesetzt, die dann ihrerseits durch Reibung in Warme ithergefithrt wird. 
In (2,12) kommt zum Ausdruck, daB die Umsetzung der Grundstrémungsenergie in Turbulenz- 
energie und die Umwandlung der Turbulenzenergie in Warme nicht am gleichen Ort stattzu- 
finden brauchen. La®t man die Schwankungsgeschwindigkeiten gegen Null gehen, so erhalt 
man aus (2,12) den von K. Wieghardt? fiir laminare Grenzschichten angegebenen Energiesatz 


d 


1 
2 


1 
2 


dx (Ui 


0 


aU\2 
(ay) dy. 


(2,12 a) 


1 ip Prarie Nachr. Akad. Wiss. Géttingen. Math.-Physikal. Klasse, 1945, S. 6. 
K. Wieghardt, Ing.-Arch, 16 (1948), S. 231. 
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3. Aufteilung der Grenzschicht in wandnahen und iuBeren Teil. In W andabstanden, die 
klein gegen die Grenzschichtdicke § sind, weicht die Schubspannung 7 nicht merklich von der 
Wandschubspannung t, ab, und die Strémungsverhiltnisse sind hier vom Druckgefiille dp/dx 
praktisch unabhangig. Wie man aus der Versuchserfahrung weib, iiben die Zahigkeit und die 
Wandrauhigkeit nur in einer an die Wand angrenzenden Schicht von der Dicke 6, einen un- 
mittelbaren EinfluB auf die Strémungsvorgange aus. Ist nun diese Dicke 6,, hinreichend klein, 
so gibt es sicher Wandabstinde y, die gréfer als 6,, sind, die aber andererseits doch sehr klein 
gegen die Grenzschichtdicke 6 sind, so daB in diesem Gebiet 5, < y <6 ein universelles Gesetz 
fiir die Grenzschichtstrémung gilt, nach welchem alle StrémungsgréBen durch nur zwei dimen- 
sionsbehaftete GréBen festgelegt sind, namlich durch die Schubspannungsgeschwindigkeit 


uh == We (3,1) 
Q 


und den absoluten Abstand y von einer Bezugsebene, die praktisch mit der Wandoberflache zu- 


sammenfallt. Der Geschwindigkeitsverlauf der Grundstrimung ist hier durch die bekannte Be- 
ziehung ! 


oU 
ge (3,2) 


gegeben. Darin ist x ~ 0,4 eine universelle Konstante. 
Die Existenz des universellen Strémungsgesetzes im Gebiet 6, < y <6 erméglicht die Zer- 
legung der Grenz- 
schichtstrémung in 
einen wandnahen Teil 
(0 <y <6) und einen 
_ auBeren Teil (y >6,.), 
Abb. 1, deren Stré- 
_mungenasymptotisch 
in die  universelle, 

turbulente Grenz- 

schichtstrémung 

iibergehen; die des 
_ ersteren mit wachsen- 
dem y, die des letz- 
teren mit abnehmen- 


dem y. Diese Auf- 


’ 
_teilung bietet den Vor- 
teil, dag man nach . 
hangig hangi 

Einfithrung des Para- Bereich Geschwindigkeitsgesetz ae peed ae vaeae 
meters v* die Stré- - aes 
-mungsvorgangein je- | mene nabs 0<y<6 | U=rrf a, ee tae U, (x) 
: eil v rauhigkeit 
dem Teil unter Re- U, (2) 
duktion der EinfluB- i] ae y, Wand- 

a i w Wy El reel (Sh wach von wate 
gréBen getrennt fiir AU | Buberer Teil Oe ei Neat | Pr renOs th, (s | rauhigkeit 


| v*/U,) 
sich (experimentell 3 v, Wand- 


dU “e v* 


universelles 


i OeEL Ow <x ¥ € 6| bee pe - | rauhigkeit, 
oder theoretisch) un Gesetz |” = ¥ <0 | eee See aaa 
tersuchen und beide 
Teile nach Bedarf zu- Abb. 1. Aufteilung der Grenzschicht in wandnahen und auferen Teil. 


-sammenfigen kann. 
Fiir den wandnahen Teil (0 < y <6) besteht ein Geschwindigkeitsgesetz von der Form 


co otg(2), (3,3) 


_wobei die Funktion f(v*y/v) im allgemeinen aufer von v* y/y noch von der Wandrauhigkeit ab- 
hangt. Die vorliegenden experimentellen Ergebnisse an glatten und rauhen Wanden lassen sich 

mit Hilfe des Mischungswegansatzes von L. Prandtl bis unmittelbar an die Wand verstehen und 
-formelmaBig darstellen®. Hier interessiert davon nur die asymptotische Form fiir y>06, , die 


1 L. Prandtl, Fihrer durch die Strémungslehre. S. 105ff. 3. Auflage. Braunschweig 1949. 
2 J. Rotta, Ing.-Arch. 18 (1950), S. 277. 
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sich aus (3,2) durch Integration ergibt: 
7 
Wee = ee Cc). (3,4) 


Die GréBe C ist eine Funktion der Wandrauhigkeit. 

Der auBere Teil (y>6,) ist von der Zahigkeit unabhangig und im wesentlichen durch zwei 
Funktionen festgelegt; namlich durch die Geschwindigkeitsverteilung U, (x) am AuBenrand und 
durch die Schubspannungsgeschwindigkeit v* (x). Wahrend die Geschwindigkeit U, (x) im all- 
gemeinen eine vorgegebene Funktion ist (z. B. Potentialstrémung um den Kérper), stellt v*(x) 
eine Funktion dar, die die Grenzschichtrechnung erst liefern soll. Zum Gliick ist aber das Ge- 
schwindigkeitsprofil des auBeren Teiles bei geeigneter Normierung nur verhaltnismabig wenig 
von v*/U, abhingig, wie noch gezeigt wird, so daB man die beschriebene Aufteilung der turbu- 
lenten Grenzschicht fiir eine rationelle Berechnung mit Erfolg anwenden kann. 
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Abb. 2. Geschwindigkeitsprofile bei konstantem Druck nach Messungen von F. Schultz-Grunow 
an glatten Wanden und von W. Tillmann an rauhen Wainden. 


Die Voraussetzungen fiir die Aufteilung in zwei voneinander unabhiangige Bereiche sind 
meistens in befriedigender Weise erfiillt. Dies ist jedoch keineswegs selbstverstandlich und be- 
darf im Einzelfall der Nachpriifung. Zu diesem Zweck sei angegeben: Die Dicke 6,, betragt bei 
glatten Wanden 6, ~ 40 »/v*; bei starker Wandrauhigkeit ist 6,, durch die Abmessungen der 
Rauhigkeiten bestimmt. Nach den Versuchen von J. Nikuradse! an sandrauhen Rohren ist 6,, 
etwa gleich der KorngréBe der Rauhigkeit, wobei die y-Werte von einer Bezugsebene aus ge- 
messen werden, in der U im Mittel verschwindet. Der Wert v*/U, und der drtliche Reibungs- 
beiwert cy stehen in folgender Beziehung miteinander: 


2 ey 3.5 
of as (3,5) 
4. Empirische Geschwindigkeitsprofile. Um nach dem Impulssatz (2.5) und dem Energie- 
satz (2,12) die Grenzschicht fiir eine vorgegebene Geschwindigkeit U, (x) berechnen zu kénnen, 
mu man entweder, wie es bei laminaren Grenzschichten geschieht, einen passenden analytischen 
Ansatz fiir das Geschwindigkeitsprofil machen, oder man muB® sich die ertorderlichen funktio- 
nalen Zusammenhange aus Versuchsergebnissen ermitteln. Dieser Aufgabe wollen wir uns jetzt 
zuwenden. Da neuerdings McSreihen vorliegen, bei denen die Wandschubspannung durch eine 
besondere Messung?? bestimmt wurde, kann eine Sichtung des vorliegenden Versuchsmateriales 
mit gr6Berem Erfolg betrieben werden, als es bisher méglich war. 


ae J. Nikuradse, Strémungsgesetze in rauhen Rohren. VDI-Forsch.-Heft 361 (1933). 


* H. Ludwieg, Ing.-Arch. 17 (1949), S. 207. 
* H. Ludwieg und W. Tillmann, Ing.-Arch. 17 (1949), S. 288. 
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eas Darstellung der duBeren Geschwindigkeitsprofile tragen wir den Wert (U, — U)/v* iiber 
yv*/d,U, auf. Die KoordinatenmaBstibe sind dabei so normiert. daB der Integralwert 


wird, wie der Vergleich mit (2,6) zeigt. Abb. 2 veranschaulicht eine solche Auftragung von Grenz- 
schichtprofilen an der ebenen Platte ohne Druckgefiille! 2, und Abb. 3 gibt Profile bei Druck- 
anstieg wieder ?®. 

Auf Grund der hydrodynamischen Differentialgleichungen laBt sich theoretisch zeigen, dah 
zwar bei turbulenter Rohrstrémung die Geschwindigkeitsverteilung? (U, — U)/v*, iiber y/r auf- 
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Abb. 3. Geschwindigkeitsprofile bei Druckanstieg nach Messungen von H. Ludwieg und W.Tillmann. Hy» = 64/09. 


getragen, das sog. Mittengesetz, vom Verhaltnis v*/U, unabhangig ist, das aber bei der Grenz- 
schicht an der ebenen Platte ohne Druckgefalle eine Abhangigkeit von v*/U, zu erwarten ist, 
auch wenn die Anlaufstrecke und die Reynoldszahl hinreichend groB sind. Die genaue Betrach- 
tung der Abb. 2 zeigt auch tatsichlich eine solche Abhangigkeit von v*/U; allerdings ist sie sehr 
gering. Nach Abb.3 wirken sich jedenfalls Anderungen im Druckgefalle wesentlich starker auf 
die Grenzschichtprofile aus. 

Fiir den duBeren Teil des Geschwindigkeitsprofiles gilt also das allgemeine Gesetz 


* 
U,—U= vt F(22), (4,1) 


in welchem die Funktion F (yv*/6,U,) hauptsachlich von U, (x) und in geringem MaBe von v*/U, 
abhangt. Fiir kleine Wandabstande y—> 6, haben die Profile die asymptotische Form 


1 yor 
Ee hed eers oe 
U,— U=v ( . In 3, van Kk) 4 (4,2) 


wobei die GriéBe K fiir jede Profilform einen anderen Betrag hat. 


1 F, Schultz-Grunow, Luftf.-Forschg. 17 (1940), S. 239. 

2 W. Tillmann, Untersuchungen iiber Besonderheiten bei turbulenten Reibungsschichten an Platten, 
Bericht UM 6627 Kais.-Wilh.-Inst.fiir Strémungsforschung. Géttingen 1945. 

3 H. Ludwieg und W. Tillmann, Ing.-Arch. 17 (1949), S. 288. 

4 U, ist in diesem Falle die Geschwindigkeit in Rohrmitte, r der Rohrhalbmesser. 


3* 
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Fir die Durchfithrung von Naherungsrechnungen geniigt es, iiber den Zusammenhang zwi- 
schen einzelnen KenngréBen (Verdrangungsdicke 6,, Impulsverlustdicke 6, u.a.) Bescheid zu 
wissen. Der nachste Schritt ist nun, diese Kenngréfen in eine Form zu bringen, die es gestattet, 
den Einflu& der Zahigkeit und der Wandrauhigkeit getrennt zu erfassen. Fiir die Impulsverlust- 
dicke 6, gilt nach (2,7) 


U U r U [ u\ 
02 ={20-%) Cs {( inlay | ( Ui) se 
0 


0 0 
dafiir 148t sich auch schreiben 
yee, 
6,=6 (1 Pa 4,3 
1 UG if ( ? ) 
Hierbei ist der Wert 
*: i poet Uf ay 
fie ae (4) 


unter den gemachten Voraussetzungen (6.< 6) praktisch unabhangig von dem Geschwindigkeits- 
verlauf bei y<6,. In gleicher Weise 
laBt sich die Energieverlustdicke 
a= nach (2,13) auf die Form 


* %*\2 
wa bp—aft a (af 


bringen, wobei 


700 


600 


f U,— U}3 /yv* 
l= | v* a(5. i) (4,6) 
0 


ebenfalls von der Profilform bei y<6,, 
praktisch unabhangig ist. Da das 
Geschwindigkeitsgesetz fiir y-Werte 
Ow Sy <6 sowohl (3,4) als auch (4,2) 
geniigt, soergibt sich durch Einsetzen 
von (3,4) in (4,2) 
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200 Ui _ 1m Re, +B, (4,7) 
x Schultz-Grunow, 8p./r=0 v x 
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Abb. 4. Zusammenhang zwischen den Profilparametern J, und I, nach (4,11) URGE UTILS aoe 
und Versuchen. B=C oh K (4,9) 
ist. 

Die GriBen I,, I, und K sind reine Formparameter, die ohne weiteres aus der Profilform her- 
geleitet werden kénnen und nicht von der Form des Wandgesetzes (3,3) abhangen, wenn die Be- 
dingung 0,.<6 hinreichend erfiillt ist. Wir machen nun die hauptsachlich durch die geringe Ab- 
hangigkeit von v*/U, gerechtfertigt erscheinende Annahme, da® zwischen diesen GriBen ein 
eindeutiger Zusammenhang besteht, den man empirisch aus vorliegenden Messungen ermitteln 
kann und der in Abb. 4 und 5 dargestelltist. Nahert man das Geschwindigkeitsprofil fiir 6> y>6 
durch den Ansatz ee 


U,—U A ¥y il 
so ae (1 Lin & (4,10) 
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an, sola ich di ai i 
. ssen sich die Zusammenhange zwischen J,, I, und B naherungsweise durch 


T 2 3 
1 Wakes 4 
zi 11 1 
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beschreiben, wobei A als Parameter auftritt. Hiermit ist es méglich, 09,63 und cp zu bestimmen 
wenn 0,, [,, Re, und die Wandrauhigkeit gegeben sind. i 


5. Turbulenzenergiestrom. Umdie x4 g 
GréBe des Turbulenzenergiestromes —_25|-39 : 
in (2,12) ermitteln zu kénnen, be- 
nétigt man das ,,Turbulenzprofil*‘, | 
| dasist die dimensionsloseAuftragung —,, | S| 
der Turbulenzenergie E/v*?  itber “J 
yv*/d,U,. Hierbei strebt E/v*? in 
Wandnahe y— 6, einem universell 
giltigen Wert zu und fallt fiir y <6,, 
sehr rasch auf Null ab. Obwohl es 
mit bekannten Hitzdrahtanordnun- 
gen grundsatzlich méglich ist, die 1 45 
quadratischen Mittelwerte aller drei 4 
Schwankungskomponenten und da- 
oe E experimentell zu bestimmen, ke | YL Esai Gearon On base 
iegen auswertbare Messungen an ot 

: = ; ? + Druckansteg 
Grenzschichten nur fiir die Kom- Lf Ee Dnchaateg 
ponente wu vor, die allerdings den re ss Druckabhal! 
wesentlichsten Beitrag zu E liefert. 
Die in Abb. 6 u.7 dargestellten Langs- 
schwankungsprofile J/u2/v* wurden 
von W. Tillmann* gemessen. Da fio 5 0 6 20 35 7) 
v*/ Uy nach Abschnitt 4 nur maBigen Abb. 5. Zusammenhang zwischen dem Verhiltnis v*/U,, der Reynoldszahl Re, 
EinfluB auf das aubere Geschwindig- und dem Profilparameter I, nach (4,11) und Versuchen an glatten Wanden. 
keitsprofil hat, so darf man auch nur 
einen geringen EinfluB8 von v*/U, auf das Turbulenzprofil erwarten. Abb. 6 bestatigt diese Ver- 
mutung fiir die Langsschwankungsprofile der Plattenstrémung ohne Druckgefalle bei glatter und 
rauher Oberflache, so daB der Gedanke nahe liegt, daB sich die mit Turbulenzprofilen zusammen- 
hangenden GréBen, der Turbulenzenergiestrom und die in (2,14) definierte Dissipationsfunktion 
D approximativ durch einen eindeutigen Zusammenhang mit dem Formparameter I, be- 
schreiben lassen. 

Die Anderung des Turbulenzenergiestromes hat bei zweidimensionaler Grenzschichtstrémung, 
soweit bisher festgestellt werden konnte, meistens keinen sehr grofen Anteil, so da eine etwas 
groBziigige Behandlung dieses Gliedes in (2,12) in der Regel statthaft 2u sein scheint. Die Gré- 
Benordnung der v- und w-Schwankungskomponenten laBt sich aus einigen alteren Messungen 
in Rohren und Kanilen? abschatzen. Nach dem vorliegenden Material kann man die Be- 


1 W. Tillmann, Uber die Wandschubspannung turbulenter Reibungsschichten bei Druckanstieg, 


Diplom-Arbeit. Gottingen 1947. 
2 H. Reichardt, Naturwiss. 26 (1938), S. 404. — H.C. H. Townend, Proc. Roy. Soc. A 145 (1934), 


S. 180. — A. Fage, Philos. Mag. 7, 21 (1936), S. 80. 


38 Rotta: Beitrag zur Berechnung der turbulenten Grenzschichten. 


Ingenieur-Archiv 


ziehung 
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Abb. 6. Liingsschwankungsprofile bei konstantem Druck nach Messungen von 
W. Tillmann an glatten und rauhen Wanden. 
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Abb. 7. Lingsschwankungsprofile bei Druckanstieg nach Messungen 
von W. Tillmann. 
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6. Dissipationsfunktion. 
Zur Bestimmung der Dissi- 
pationsfunktion D nach 
(2, 14) wollen wir ebenfalls 
vermittels der im Abschnitt 
3 dargelegten Ergebnisse 
den Einflu8 der Zahigkeit 
und der Wandrauhigkeit 
fiir sich erfassen. Den Bei- 
trag des wandnahen Ge- 
bietes (0< y <6) zu (2,14) 
kann man leicht ausrech- 
nen, da man hier t/9 = v** 
= konst. setzen und die mit 
der Grundstrémung fort- 
gefiihrte Turbulenzenergie 
vernachlassigen darf. Die 
in einem Streifen von der 
Breite y dissipierte Energie 
ist dann gleich der von der 
Grundstrémung geleisteten 
Arbeit zuziiglich der in 
diesen Streifen von aufen 
hinein diffundierenden ki- 
netischen Energie. Man er- 


halt also fiir y <6 
: 
S dy’ = v™ U(y) 


——O(y). 


Liegt die obere Integrations- 
grenze yim Gebiet der uni- 
versellen Grenzschichtstré- 
mung (0, < y <6), so kann 
man U (y) nach (3,4) ein- 
fiihren und Q(y) hat aus 
Griinden der fiir die Stré- 
mung geltenden Ahnlichkeit 
einen von y unabhangigen 
Wert Q = — c, v*, wobei 
cqeine universelle Konstante 
ist. Somit ergibt sich fiir 
QE? S10 


(6.1) 


Is dy’ = v*8 | 
0 


] * 
i. In aoe C+ c) : | 


Bei dem auBeren Grenz- 


schichtteil (y >6,,), bei dem 


(6,2) 
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die Reynoldszah] der Turbulenz groB ist, gilt fiir die Dissipation (2,10) die bekannte Beziehung 1! 


E3/2 
eatin (6,3) 


in der / eine die groBen Turbulenzelemente kennzeichnende Lange ist und c ein dimensionsloser 


Faktor, der nur von der Struktur der Turbulenz abhingt. Die Gleichung (6,3) folgt aus der Tat- 


sache, da8 die im wesentlichen in den gréBten Tur- QD 
bulenzelementen enthaltene kinetische Energie EF ,,| v*? _ 
laufend in immer kleinere Turbulenzelemente iiber. 
gefiihrt wird und schlieBlich in den kleinsten Ele- 
menten in Warme umgewandelt wird. Im Gebiet  23/—— 
der universellen Grenzschichtstrémung ist aus 
Ahnlichkeitsgriinden | proportional zu y, und E ist | 
konstant proportional v*2, Ferner ist hier die Dissi- “ Eb 
pation gleich der der Grundstrémung entzogenen 
Energie, so daB also fiir 6,, << y <6 a 
E323 : ; 
eat Cea ny (0,4) D=v"*{-06+5 75 log Rey) 
gilt. Damit ergibt sich dann as 
ae ; Bt ay =we(s,— 21 a 
| Mee ailcs ak oe Oy” re, 
yi i (6,5) 
wobei der Wert des Integralausdrucks oA | log Rey 
— | c 8) : (5, =) ro In fae (6.6) — Druckgefiille ta ‘Abheagigket voter Beal ee 
J ae eat sal MeBergebnissen von F. Schultz-Grunow. 
an 0, U, I’ 
760 7 2 
von der unteren Integrations- F-Lehvey § 
grenze y unabhangig ist, wenn sie yy le 
im Bereich 6,, < y <6 liegt. Die 160 
_ Addition von (6,2) und (6,5) 
liefert schlieBlich no} 
x» 10 
[Dy | S dy | ale D=v* (6+575 log Re,/) 
0 | (6,7) 20 = 
— y*3 (= In Re, + c] | Pi 
mit Mire 
=C+J,+¢,. (6,8) 
Da itber das Verhalten der Funk- a 20 € 
tionen / und ¢ aufer in der Ge- 
gend 6, <y <6 nichts bekannt A 
ist, kann J, auch bei bekanntem 60 : 
Turbulenzprofil nicht aus (6,6) | 
berechnet werden. Es bleibt da- 40 
her nur die Méglichkeit, durch 1) ies 
Einsetzen der versuchsmabig be- 
5 5 c 20 
stimmbaren Gréfen in die Energie- | 29 
gleichung (2,12) die Funktion D # 
durch Differentiation auszurech- 0 wate x L 
nen, ein Verfahren, das mit einigen a sai ee ee meee ee ey 
Unsicherheiten behaftet ist. Das W | 
in Abb. 8 dargestellte Auswer- “20 “up 


tungsergebnis der MeBreihe von Abb. 9. Zusammenhaung zwischen der Dissipationsfunktion nach (6,7), der Reynolds- 


- 2 Platte zahl Re, und dem Profilparameter J, nach Auswertung von Mefergebnissen an 
F.Schultz-Grunow? an der Plat 1 eee 


dU \? sai a tess 
1 Vel. FuBnote 1 von S. 32. Der Beitrag » ar der Grundstrémung ist hier zu vernachlassigen. 
2 Vel. FuBnote 1 von S. 35. 
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ohne Druckgefille bestatigt die Richtigkeit der Beziehung (6,7) hinsichtlich des Reynholdszan ci | 
flusses. Das Ergebnis der Auswertung mehrerer Versuchsreihen von H. Ludwiegund W. Tillmann ist 
in Abb. 9 dargestellt und gibt den Zusammenhang zwischen der Dissipationsfunktion und den 
Profilparameter I, an. Die Streuung der Versuchspunkte ist z. T. auf die in dem Auswertungs- 
verfahren begriindeten Unsicherheiten zuriickzufihren. Kine Umrechnung auf die Verhaltnisse von 
rauhen Wandenist an Hand von (6,8) durch Einsetzen des abgednderten Wertes Cvon (3,4) méglich. 

7. Ahnliche Lésungen. Fiir 
Ve die Weiterentwicklung ist 
; al Zi neben Versuchen zur Erhar- 
tung der angegebenen Zu- 
sammenhinge auch die theo- 
retische Durchdringung not- 
wendig. Mit der Grenzschicht- 
——— St t gleichung (2,1) und derEnergie- 
bilanz (2,9) ware unter Hinzu- 
nahme einiger Hypothesen die 
: | ; =) 7) Berechnung der turbulenten 
Grenzschicht —grundsatzlich 
denkbar. Man hatte dabei 
zwei partielle, nicht lineare 


3 
ae i Re,~10" Differentialgleichungen ane 
SS = = statt einer bei der laminaren 
Grenzschicht zu integrieren. 
Re,=10° m Daher hat man bei den turbu- 
-4¢ -02 oi 02 a4 06 08 10 12 144 lenten Grenzschichten auch in 


mathematischer Hinsicht gré- 
Bere Schwierigkeiten zu _ er- 
warten als bei den laminaren, 
bei denen die mathematischen 
Probleme keineswegs ‘einfach 
zu nennen sind. Man wird 
deshalb zunachst nach ein- 
fachen LéisungstypenAusschau 
halten, an denen sich verschie- 
dene Kinzelfragen  studieren 
lassen. In diesem Zusammen- 
hang ist es interessant, daB es 
unter gewissen Voraussetzun- 
Abb. 10. Abnliche Lésungen der Gleichungen fiir die turbulente Grenzschicht an genauch bei turbulenten Grenz- 
glatten Wanden. schichten sog. ahnliche Lésun- 
gen gibt, d. h. Lésungen, bei 
denen das Geschwindigkeitsprofil langs der Wand nur affin verzerrt wird; die partiellen Diffe- 
renzialgleichungen lassen sich dabei in gewéhnliche verwandeln. 

Zunachst laBt sich zeigen, daf fiir den auBeren Teil (y > 0,,) unter Vernachlassigung der Zahig- 
keit in den Grenzschichtgleichungen ahnliche Léisungen bestehen, wenn die Geschwindigkeits- 
verteilung am Aufenrande der Beziehung 

U,;=ax", (7,1) 
wobei aund m konstant sind, geniigt und der értliche Reibungsbeiwert c7 = konstant vorgegeben 
ist. Die geometrische Ahnlichkeit des Strémungsbildes bedingt, daB die Grenzschichtdicke hier- 
bei linear mit » anwachst. Das Geschwindigkeitsprofil laBt sich als (U,— U)/v* iiber y/x dar- 
stellen und ist von den beiden Parametern m und v*/U, abhangig. Fiir kleine y/x-Werte hat sie die 
asymptotische Form 


-0%  -02 0 02 04 96 08 10 12 a 


U,—U 1 y cmlee® 
i = - In a + K (m. | ; (7,2) 
wobei sich die Konstante K(m, v*/U,) bei gegebenen Werten m und v*/U, aus den Grenzschicht- 
gleichungen bestimmt. 


1 Vel. FuBnote 3 von S$. 34, sowie unverdffentlichte Versuchsreihen. 
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Die nur fiir Wandabstinde y>6,, giiltige Lésung mu durch das Wandprofil (3,3) erganzt 
werden, um daraus die vollstiindigen Geschwindigkeitsprofile zu erhalten. Die Bedingung fir 
den stetigen Anschlu8 des auBeren Teiles an (3,3) erhalt man durch Elimination von U/v* in 
(7,2) mit Hilfe von (3,4): 

U. * 

ee =n ye — K(m, | Bei 20, (7,3) 
Die diskutierten ahnlichen Lésungen des auBeren Teiles haben nur dann reale Bedeutung, wenn 
die Reynoldszahl U,x/» und die Wandrauhigkeit, deren Wirkung in der GriéBe C zum Ausdruck 
kommt, derart sind, daf (7,3) fiir alle x-Werte identisch erfiillt ist. Bei sehr groBen Reynolds- 
zahlen ist dies der Fall, wenn die die Rauhigkeit kennzeichnende Lange (z. B. KorngréBe bei 
Sandrauhigkeit) zu x proportional ist. Bei hydrodynamisch glatten Wanden und bei konstanter 
Rauhigkeit, bei denen C eine Konstante ist, kann die Bedingung (7,3) nicht fiir alle x-Werte 
streng erfiillt werden. Da x in (7,3) logarithmisch auftritt, darf man fiir hinreichend groBe x- 
Werte den Ausdruck auf der rechten Seite von (7,3) zumindest intervallweise als konstant und 
somit die Bedingung fiir die ahnlichen Lésungen mit praktisch ausreichender Genauigkeit als 
erfillt ansehen. Wichtig ist aber, und das gilt auch fiir die Entwicklungen der Abschnitte 4 
bis 6, daB 6,, so klein ist, daB die Funktion (U, — U)/v* an der Stelle y = 0, nicht merklich von 
der Asymptote (7.2) bzw. (4,2) abweicht. 

Eine Abschatzung der ahnlichen Lésungen ist vermittels des Impulssatzes (2,5), des Energie- 
satzes (2,12) und der empirischen Zusammenhange von Abschnitt 4 bis 6 méglich. In Abb. 10 
sind die Ergebnisse einer solchen Rechnung zusammengestellt, wobei die Verhaltnisse glatter 
Wande zugrundegelegt wurden. Interessant ist, da® nicht fiir alle m-Werte eine physikalisch 
sinnvolle Lésung existiert . Dies folgt aus dem Impulssatz (2,5), wenn man (7,1) einfiihrt. Wegen 
U/U, <1ist nach (2,6) und (2,7) stets 6, >6,. Negative Werte von 6, und tT, kénnen bei posi- 
tivem U, nur auftreten, wenn sich die Strémung von der Wand ablést. In diesem Falle sind aber 
die Voraussetzungen der Grenzschichttheorie nicht mehr erfiillt. Nach dem Impulssatz muf 


daher m> — 1/3 sein. Die Abb. 10 zeigt, daB die Ablésung praktisch etwa bei m = — 0,2 zu er- 
warten ist. Zum Vergleich sei daran erinnert, das bei den entsprechenden ahnlichen Lésungen 
der laminaren Grenzschichten die Ablésung bei m = — 0,091 erfolgt. Hierdurch wird die be- 


kannte Erfahrungstatsache bestatigt, da® turbulente Grenzschichten einen gré®eren Druck- 
anstieg als laminare iiberwinden kénnen. 


8. Zusammenfassung und Schlu8. Die turbulente Grenzschicht 1a8t sich vorteilhaft in zwei 
Gebiete aufspalten (Abb. 1), welche in einem gemeinsamen Zwischengebiet (6, < y<6) anein- 
ander anschlieBen. Im auBeren Gebiet (y >6,,) hangt das Geschwindigkeitsprofil im wesentlichen 
nur von der Geschwindigkeitsverteilung U, am Aufenrand ab und la8t sich naherungsweise als 
eine der Schubspannungsgeschwindigkeit v*¥ = \t/e proportionale Funktion darstellen, ist aber 
sonst von der Zahigkeit und den Vorgangen an der Wand unabhangig. Im wandnahen Teil 
(0 <y <6) ist das Geschwindigkeitsprofil durch die Zahigkeit und die Wandrauhigkeit bedingt, 
aber vom Geschwindigkeitsverlauf der AuSenstrémung nur indirekt iiber den Parameter v* 
abhangig. Im Zwischengebiet herrscht ein universelles Strémungsgesetz, das nur von v* ab- 
hangt, aber sowohl von der Wandrauhigkeit und Zahigkeit, als auch von der AuBenstr6mung un- 
beeinluBt ist. 

Unter Beriicksichtigung dieser Erkenntnis liefert die Auswertung vorliegenden Versuchs- 
materials empirisch die Zusammenhange zwischen den verschiedenen Gr6Ken, die mit den an- 
gegebenen Integralsatzen fiir das Impuls- und Energiegleichgewicht zur angenaherten Berech- 
nung turbulenter Grenzschichten dienen kénnen. Ein hierauf aufgebautes Berechnungsverfah- 
ren bietet gegeniiber den bisher bekannt gewordenen Verfahren den Vorteil einer zuverlassigeren 
Bestimmung des Reibungsbeiwertes in Abhangigkeit von der Reynoldszah| und dem vorgegebenen 
Geschwindigkeitsverlauf, sowie der Anwendbarkeit bei rauhen Wanden. Die Geschwindigkeits- 
profile sind dabei zweiparametrig, da an den einparametrig vorausgesetzten duBeren Teil (4,1) 
je nach der vorhandenen Reynoldszahl ein Wandgesetz (3,3) angepaBt wird. Fiir die praktische 
Durchfithrung von Grenzschichtrechnungen ist es notwendig, die gegebenen Zusammenhange in 
eine bequemere Darstellung zu bringen, was ohne Schwierigkeit méglich ist. 


Aus dem Max-Planck-Institut fiir Strémungsforschung Gottingen. 
(Eingegangen am 1. Juli 1950.) 
Anschrift des Verfassers: J. Rotta, (20b) Gottingen, Hann., Bottingerstrabe 6. 
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Uber die theoretische Berechnung der Strémungsverluste eines ebenen 
Schaufelgitters*. 


Von H. Schlichting und N. Scholz. 


1. Uhersicht. In der flugtechnischen Aerodynamik sind in den letzten beiden Jahrzehnten 
Verfahren zur theoretischen Berechnung des Widerstandes eines einzelnen Tragfliigelprofiles 
entwickelt worden. Diese Verfahren werden hier erstmalig auf das Profil in Gitteranordnung, 
das sog. Schaufelgitter, tibertragen, das ein wichtiges Bauelement jeder Strémungsmaschine 
bildet. Es wird damit méglich, die Strémungsverluste (Wirkungsgrad) eines Schaufelgitters in 
Abhingigkeit von seinen geometrischen und strémungsmechanischen Parametern (Schaufelform, 
Gitterteilung, Staffelungswinkel, Anstellwinkel, Umlenkungswinkel) systematisch theoretisch 
zu untersuchen. Die rechnerischen Grundlagen hierzu werden hier fiir das ebene Schaufelgitter 
in inkompressibler Strémung auf grenzschicht-theoretischer Grundlage entwickelt. 

Fiir die Durchfithrung der Rechnung ist nacheinander fiir das vorgelegte Schaufelgitter zu 
ermitteln: 

1. die potentialtheoretische Druckverteilung lings der Schaufelkontur, 

2. die (laminare oder turbulente) Reibungsschicht langs der Schaufelkontur, 

3. die Nachlaufstrémung hinter der Schaufel. 

Das im Gitter umgesetzte Druckgefiille ist von entscheidendem KinfluB auf die Ausbildung 
der Reibungsschicht und das Eintreten von Ablésung, und damit auf die GréBe der Strémungs- 
verluste. Beim Beschleunigungsgitter, welches Druck in Geschwindigkeit umsetzt (Turbine), 
bleibt die Reibungsschicht diinn, und ihre Ablésung kann meist vermieden werden, so daB die 
Verluste klein bleiben. Dagegen sind beim Verzégerungsgitter, welches Geschwindigkeit in Druck 
umsetzt (Kompressor), die Reibungsschichten wesentlich dicker. Infolgedessen ist die Ablésungs- 
gefahr hier gréBer, und die Verluste nehmen erheblich zu. — Das Verfahren wird erprobt an 
einigen Beispielen von Beschleunigungs- und Verzégerungsgittern mit sehr diinnen Schaufeln. 


2. Einleitung. Fiir axial durchstrémte Rader von Strémungsmaschinen kann die Strémung 
im Laufrad und Leitrad in bekannter Weise durch ein ebenes Schaufelgitter idealisiert werden. 
Dieses entsteht dadurch, daB man einen koaxialen Zylinderschnitt abgewickelt denkt. Bisher 
sind die Strémungsverluste solcher Gitter ausschlieBlich empirisch ermittelt worden. Wir stellen 
uns nun hier die Aufgabe, die Strémungsverluste eines ebenen Gitters rein theoretisch voraus- 
zuberechnen. Dabei legen wir fiir die Berechnung der Gitterstémung die tragfliigeltheoretische 
Auffassung! zugrunde, bei welcher die einzelne Schaufel als ein Tragfliigel aufgefaBt wird, und 
die gegenseitige Beeinflussung der Schaufeln eine wesentliche Rolle spielt. Die theoretische 
Ermittlung der Strémungsverluste eines Gitters wird bei dieser Auffassung erméglicht durch 
Rechenverfahren, welche fiir den Einzelfliigel in der flugtechnischen Aerodynamik in den letzten 
beiden Jahrzehnten zu groBer Vollkommenheit entwickelt worden sind, aber bisher bei der 
Berechnung von Strémungsmaschinen noch kaum Eingang gefunden haben. Ihre Anwendung 
diirfte hier jedoch ebenso erfolgreich sein wie in der Flugtechnik. 

Wahrend die bisher bekannten Rechenverfahren fiir Gitterstrémungen sich auf die Ermittlung 
der reibungslosen Potentialstrémung beschranken und deshalb keinerlei Aussagen itber auf- 
tretende Verluste machen kénnen, wollen wir fiir die Gitterstrémung die Fliissigkeitsreibung nach 
den Rechenverfahren der Prandtlschen Grenzschichttheorie beriicksichtigen. Das bedeutet, daB 
wir nur in einer diinnen wandnahen Schicht Reibungskrafte annehmen, wahrend das ibrige 
Strémungsfeld von der reibungslosen Potentialstrémung beherrscht wird. Die Berechtigung 
dieser Annahme fiir sehr groBe Reynoldssche Zahlen wird in der Grenzschichttheorie ein- 
gehend begriindet und kann als allgemein bekannt gelten 2. 

Rein empirisch ist seit langem bekannt, und auch physikalisch ist es plausibel, daB® die Stré- 
mungsverluste im Verzégerungsgitter (Umsetzung von Geschwindigkeit in Druck, Verdichter) 
erheblich gréBer sind als im Beschleunigungsgitter (Umsetzung von Druck in Geschwindigkeit, 


* Die vorliegenden Untersuchungen wurden durch die Notgemeinschaft der deutschen Wissenschaft 
unterstiitzt. 

: Im Gegensatz zur Auffassung der Gitterstrémung als Kanalstrémung. 

* Man vel. z. B. H. Schlichting, Grenzschicht-Theorie, erscheint demniichst bei G. Braun, Karlsruhe. 
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Turbine). Das hier entwickelte Verfahren gibt eine rationelle Theorie fiir diese wichtige Er- 
scheinung. Uberdies ergibt sich bei weiterem Ausbau der Theorie die Méglichkeit, den Einflu8 
aller wesentlichen Parameter eines Gitters (Schaufelform, Teilungsverhiltnis, Staffelung, An- 
stellwinkel, Umlenkungswinkel) auf die Strémungsverluste systematisch zu untersuchen. Damit 
wird es erméglicht, einmal ein wichtiges Ordnungsprinzip fiir systematische experimentelle Unter- 
suchungen zu geben, und zum anderen vielleicht spater einmal das fiir einen vorgegebenen Ver- 
wendungszweck giinstigste Gitter rein theoretisch vorauszubestimmen. 

Die Strémungsverluste in einem Schaufelgitter entstehen zum iiberwiegenden Teil in den 
diinnen Reibungsschichten lings der Schaufelkontur. Einige weitere Energieverluste entstehen 
aber auch in dem Raum hinter dem Gitter. Hier bilden die Grenzschichten der einzelnen Schau- 
feln ,,Dellen* in der Geschwindigkeitsverteilung (Abb. 7). Diese Geschwindigkeitsunterschiede 
gleichen sich durch turbulente Vermischung stromabwarts aus, wobei ein weiterer Energieverlust 
auftritt. 

Fiir die Durchfiihrung unseres Verfahrens sind drei str6mungsmechanische Teilaufgaben zu 
lésen: 

1. Ermittlung der potentialtheoretischen Druckverteilung langs der Schaufelkontur (reibungs- 
lose Strémung); 

2. Berechnung der (laminaren und turbulenten) Grenzschicht langs der Schaufelkontur; 

3. Berechnung der Nachlaufstrémung hinter der Schaufel. 

Fiir die Berechnung der potentialtheoretischen Druckverteilung eines Schaufelgitters exi- 
stieren grundsatzlich zwei Methoden: die konforme Abbildung und die Singularitatenmethode. 

Bei der konformen Abbildung wird die Strémung um die unendliche Reihe von kon- 
gruenten Schaufeln durch eine komplexe Abbildungsfunktion auf die bekannte Kreiszylinder- 
strémung zuriickgefiihrt. Wesentliche Beitrage zu dieser Methode sind von F. Weinig! geliefert 
worden. Diese Methode, die ihrem Wesen nach ausschlieBlich auf das ebene Problem beschrankt 
ist, ist zwar nach J. E. Garrick? und W. Traupel® fiir beliebige Schaufelformen grundsatzlich 
als gelést anzusehen. In der rechnerischen Durchfiihrung ist sie jedoch, selbst fiir einfache 
‘Schaufelformen, so kompliziert, daB sie fiir unsere Zwecke, wo gréBere Reihen von Schaufel- 
gittern zu rechnen sind, praktisch nicht in Frage kommt. 

Bei der Singularitatenmethode, die grundsatzlich auch auf den dreidimensionalen Fall 
anwendbar ist, wird die Schaufelkontur durch Singularitaten (Wirbel, Quellen und Senken) 
ersetzt. Dieses Verfahren ist zuerst von A. Betz*° fiir die Gitterstrémung angewendet worden. 
Bei der einfachsten Form dieses Verfahrens wird jede Schaufel lediglich durch einen Einzelwirbel 
ersetzt. Der nachst bessere Ersatz ist eine stetige linienhafte Zirkulationsverteilung iiber die 
Schaufeltiefe, wobei man .,,unendlich diinne“‘ gewolbte Schaufelkonturen erhalt. Unter Hinzu- 
nahme von Quellen und Senkenverteilungen erhalt man dann Schaufeln endlicher Dicke. Wesent- 
liche Beitrage zu diesem Singularitatenverfahren sind von M. Schilhansl®, J. Ackeret’, E. Pistolesi®, 
V. Lieblein® und R. A. Spurr und H. J. Allen’ geliefert worden. Die Singularitatenmethode 
bietet die Méglichkeit, in einfacher Weise vom Einzelfliigel auf die Schaufel im Gitterverband 
iiberzugehen. Da die Strémung um das Einzelprofil als gelést angesehen werden kann!) 1% 18:14, 
ist damit bereits der in mathematischer Hinsicht schwierigste Teil der Aufgabe erledigt. 

1 #F. Weinig, Die Strémung um Schaufeln von Turbomaschinen. Leipzig 1935. 
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Die Ermittlung der potentialtheoretischen Druckverteilung an der Schaufel soll im folgenden 
nur ganz kurz gestreift werden, da sie im Rahmen unserer Untersuchungen lediglich den Aus- 
gangspunkt fiir die Berechnung der Strémungsverluste bildet. Wegen naherer Ausfiihrungen zu 
diesem Punkt sei auf eine Arbeit des zweiten Verfassers’* verwiesen. 

Die Kenntnis der potentialtheoretischen Geschwindigkeitsverteilung ermoglicht die Berech- 
nung der Grenzschicht an der Schaufelkontur nach Methoden, die sich von den fiir Einzelfliigel 


entwickelten nicht unterscheiden. Diese Grenzschichtrechnungen liefern langs der Schaufel- | 


kontur die Geschwindigkeitsverteilung innerhalb der Reibungsschicht bzw. gewisse, das Ge- | 
schwindigkeitsprofil charakteri- | 


sierende GréBen (Impulsverlust- 


Die Ermittlung des Wider- 


Tops Te7 "a 5[ve der Hinterkante des Fliigels zu- 
riickgefiihrt werden*®*°, In ahn- 
licher Weise kann nun auch der 

| Widerstand der Gitterschaufel | 
TUL co 


in der Reibungsschicht an der 
Hinterkante der Schaufel aus- 
gedriickt werden. Die Nachlauf- 
ines | strémung hinter den Schaufeln 

daha wird dabei fiir das Gitter in 
| gleicher Weise wie beim Einzel- 


> 
8 
a 
Sey 


Gittertront| 


w/e ™y  niaherungsweise ermittelt. Beim 

b “| 7, unendlich langen Gitter tritt — 
aber auch eine Anderung der 
ae Strémung sehr weit hinter dem 
-----— Gitter auf, indem sich der Druck 


Oey und die Geschwindigkeit nach 
as at GréBe und Richtung gegeniiber 
eee TD: den potentialtheoretischen Wer- 
y ten andern. Diese Anderungen von 
Druck, Geschwindigkeitsbetrag 
und Geschwindigkeitsrichtung 
gegenittber den Werten bei rei- 


dy Lites bungsloser (verlustloser) Stré- 
mung stellen die ,,Strémungsver- 
Abb. 1. Strémung durch ein ebenes Schaufelgitter. luste“‘ dar und sollenim folgenden 


vorausberechnet werden. 

3. Krifte und Geschwindigkeiten beim ebenen Schaufelgitter nach dem Impulssatz. Wir 
betrachten in diesem und dem folgenden Abschnitt die reibungslose inkompressible ebene Stré- 
mung durch ein ebenes Schaufelgitter, welches aus lauter kongruenten Schaufeln mit dem Ab- 
stand t (= Gitterteilung) besteht (Abb. 1). Die geometrischen Verhaltnisse eines solchen Gitters 
sind gegeben durch das Teilungsverhaltnis t/l (oder auch t/b), durch den Staffelungswinkel /, und 
das ,,Schaufelprofil**. Weit vor und weit hinter dem Gitter herrscht eine homogene Strémung 
mit den Geschwindigkeiten w, bzw. w,, den Richtungen f,, 6, und den Driicken p, bzw. py. 


1 N. Scholz, Ein einfaches Singularitatenverfahren zur Erzeugung von Schaufelgittern. Vortrag GAMM 
1950, Z. angew. Math. Mech. 30 (1950), S. 262. 

2 N. Scholz, Ein Singularitatenverfahren zur Berechnung der Druckverteilung von ebenen Schaufel- 
gittern. Unveréffentlichter Bericht d. Inst. f. Stromungsmech., der T. H. Braunschweig. 

3 J. Pretsch, Zur theoretischen Berechnung des Profilwiderstandes. Jb. 1938 dtsch. Luftf.-Forschg. 
oh LU lly 

4 H. Helmbold, Ing.-Arch. 17 (1949), S. 273. 


> H. B. Squire u. A, D. Young, The calculation of the profile drag of aerofoils. Brit. ARC Report 1838 
(1938). 


dicke, Wandschubspannung u.a.). | 


standes von Einzelfliigeln konnte | 
auf diese GrenzschichtgréBen an | 


durch das Geschwindigkeitsprofil | 


| fliigel mit Hilfe des Impulssatzes _ 
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Der Zweck des Schaufelgitters ist, entweder Druckenergie in Geschwindigkeitsenergie oder 
Geschwindigkeitsenergie in Druckenergie umzusetzen. Im ersteren Fall ist Po<p, und w,>w,, 
und das zugehérige Gitter hei®t ein Beschleunigungs- oder Turbinengitter. Wird andererseits 
Geschwindigkeitsenergie in Druckenergie umgesetzt (p,> pj, W, <W,), so liegt ein Verzégerungs- 
oder Pumpengitter vor. In beiden Fallen wird angestrebt, daB die Energieumsetzung mit méglichst 
geringen Energieverlusten vor sich gehen soll. Den Grenzfall p,—p, und w,—w, bildet das 
»,Gleichdruckgitter‘*, bei welchem lediglich eine Umlenkung der Strémung aber keine Energie- 
umsetzung stattfindet. Auch diese reine Umlenkung soll mit méglichst geringen Verlusten verlaufen. 

Da bei der Durchstrémung des Gitters wegen der Kontinuitat die Geschwindigkeitskompo- 
nente senkrecht zur Gitterachse vor und hinter dem Gitter gleich groB ist, #1... Wee. == Web» 
folgt, da mit der Anderung des Geschwindkeitsbetrages beim Durchgang durch das Gitter auch 
stets eine Anderung der Geschwindigkeitsrichtung verbunden ist. Diese Verhiltnisse sind aus 
dem ,,Geschwindigkeitsdreieck** (Abb. 1) zu ersehen. 

Bei der Durchstrémung des Gitters werden von der Fliissigkeit Krafte auf die Gitterschaufeln 
tbertragen. Fiir die ideale Fliissigkeit 14Bt sich der Zusammenhang zwischen den Geschwindigkeiten 
weit vor und hinter dem Gitter und den Kraften an der Gitterschaufel in einfacher Weise aus dem 
Impulssatz und der Bernoullischen Gleichung angeben. Fiir den spiteren Gebrauch seien diese 
Beziehungen hier kurz zusammengestellt. Bedeuten P, und P, die x bzw. y-Komponente der 
Schaufelkraft (fiir eine Schaufel) pro Langeneinheit der Schaufel, so gilt mit den Bezeichnun- 
gen nach Abb. 1 (9 = Dichte) 


Pia ot = (Wey — Wyy), (1) 
P,= — 0t,,, (Wey — Wy). (2) 
Unter Einfithrung der Zirkulation der Schaufel 
P= t (wyy — py) (3) 
{aBt sich dies schreiben in der Form 
Pe, ee ree (4) 
Po Oph (5) 


Die resultierende Schaufelkraft A — ye + Py, die als Auftrieb bezeichnet werden mége, ist 
somit (Formel von Kutta-Joukowsky) : 

A= ow. I. (6) 
Dabei bedeutet w., das vektorielle Mittel der Geschwindigkeiten w, und w, (Abb.1) mit den Kom- 
exten Deu, a Die Richtung von A ist senkrecht zur Richtung von w,. 


Zwischen den Geschwindigkeitskomponenten und den Richtungswinkeln bestehen die Be- 
nehungen : 


F; Wx < WX x0 3 _ Wx 
sin Bj = — >, Bit hy == ee sin B=? | 
1 


(7) 
ctg Bi —= —", ote p.=—™, ctg B.. => (ctg By +etg fy). | 


Wx oo WX 50 
Die Zirkulation nach (3) 1a8t sich auch schreiben 
T= tw,,, (ctg B, — etg f,). (8) 


Fiihrt man fiir die resultierende Schaufelkraft A in gleicher Weise wie beim Einzeltragfligel 
sinen dimensionslosen Auftriebsbeiwert ein durch 
; A 
o, ear oy a? (9) 
> w, | 
30 gilt fiir diesen nach (6), (7), (8) 
C= 2 — sin B,, (ctg by — ctg f,). (10) 
Das im Gitter umgesetzte Druckgefalle Ap = p, — pz, welches aus der Bernoullischen Gleichung 
sxrmittelt werden kann, ]a4Bt sich durch die Richtungswinkel in folgender Weise ausdriicken: 
ee ee 1 =~ 1 : (11) 
Q sin?B, sin? B, 


SS ayy 
9 W x00 


+e 
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Im Sonderfall B, =a — B, wird py= po. Dies ist das sog. ,,Gleichdruckgitter“, bei welchem 


lediglich eine Umlenkung der Strémung, aber keine Energieumsetzung stattfindet. 


Die Formeln (4), (5) und (6) fiir die Schaufelkrafte sind véllig unabhangig von der individuellen | 


Schaufelform. Sie enthalten lediglich die Geschwindigkeiten weit vor und hinter dem Gitter. 
Die Besonderheiten des Strémungsfeldes in der Nachbarschaft der Schaufeln, die stark durch die 


individuelle Schaufelform (sowie durch Teilungsverhaltnis und Staffelungswinkel) bestimmt | 


werden, gehen bei den bisherigen Betrachtungen gar nicht ein. 


Das gleiche Formelsystem erhalt man auch aus der Eulerschen Turbinentheorie (Stromfaden- 


theorie), welche mit den Begriffen der ,,unendlichen Schaufelzahl* und der ,,schaufelkongruenten‘ | 


Strémung arbeitet. Auch diese Theorie vermag die Einfliisse der Schaufelform, der Teilung und | 


des Staffelungswinkels nicht zu beriicksichtigen. 


Das Ziel unserer Untersuchungen ist die Ermittlung der Strémungsverluste im Gitter in | 
Abhangigkeit von den soeben angegebenen geometrischen GréBen und von dem Zustrémwinkel. 
Da die Strémungsverluste haupt- | 


(2) are rae ; 5 
sichlich in den Reibungsschichten | 
b am iN an den Schaufeln entstehen, und | 
a (| \ : 
da deren Ausbildung stark von 
L y(a)- der Druckverteilung und damit 
beg oa roa von der Schaufelform abhangig 
T | ist, ist es erforderlich, sich fiir die 
rie) Y Ermittlung der potentialtheoreti- 
f° t = Am t schen Druckverteilung ein be- | 
aa Uy) | pa ee quemes Rechenverfahren zurecht 
7) 
+ P od zu legen, welches es gestattet, 
- den Einflu®B der geometrischen 
t Ros t Parameter des Gitters in ein- 
| Veta Aiea facher Weise zu erfassen. 
Dee eee ceeeeg le bo ee 2 4. Ermittlung der Druckver- 
Ire 2 I 2 2 2 £ 
= b 4 s+ b = teilung an der Schaufelkontur. 
Abb. 2. Ersatz eines ebenen Schaufelgitters. Bei Verwendung der Singulari- 
a) dureb eine unendliche Reihe von Einzelwirbeln der Starke J’. ch se ies ° 
b) durch eine Wirbelbelegung der Zirkulationsdichte y (x). tatenmethode 1aBt sich em Schau- 
+b/2 felgitter in erster Naherung durch 
1 = ie y(x) dx Gesamtzirkulation einer Schaufel. eine unendliche Reihe von Einzel- 
y(x) — Schaufelkontur; t = Gitterteilung; b = Gitterbreite. wirbeln langs der Gitterachse er- 


setzen (Abb. 2a). Jeder Wirbel 
entspricht nach Lage und Starke der Zirkulation einer Schaufel. Das Geschwindigkeitsfeld 
eines auf der y-Achse im Punkt vt gelegenen Wirbels ist gegeben durch die komplexe Ge- 


aes : : i 
schwindigkeit w, (2) = w,—iwy = r= ae Samtliche Wirbel der unendlichen Wirbelreihe 
y=0, +1,+2,...,+ © ergeben somit das Geschwindigkeitsfeld 
ir wo , 
CIC rere ps Zoe (ep | 
OAD 
t 


Die unendliche Summe ist die Partialbruchzerlegung des hyperbolischen Cotangens ctgh, so | 


da8 man nach Ausfiihrung der Summation als Geschwindigkeitsfeld des Schaufelgitters mit 
der Teilung t hat 


eee fC 
1o(2) = * etgh fk (13) 
Sehr weit vor und sehr weit hinter dem Gitter (x= —- oo) hat man wegen ctgh (= co) = 1 
oe if 
14 ——= = = 106) <5 W 24 — =e Ot (14) 


in Ubereinstimmung mit (3). Die Zirkulation der Schaufeln erzeugt weit vor und weit hinter 
dem Gitter eine Geschwindigkeit parallel zur Gitterfront, die vor und hinter dem Gitter gleichen 
Betrag, aber verschiedenes Vorzeichen hat. 

Das durch (13) gegebene Geschwindigkeitsfeld stellt eine brauchbare Naherung dar fiir maBige 
und groBe Abstande von den Schaufeln. Am Ort des Schaufelgitters wird es jedoch unbrauchbar 
wegen der unendlichen Geschwindigkeiten in den Wirbelpunkten z—=- iyt. 
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Kine zweite wesentlich bessere Naherung fiir das Geschwindigkeitsfeld des Schaufelgitters 
erhalt man, wenn man nach der Art der Birnbaumschen Theorie des diinnen Kinzelfliigels jede 
, Schaufel durch eine kontinuierliche Wirbelbelegung lings der Tiefe mit der Zirkulationsdichte y (x) 
ersetzt (Abb. 2b). Dabei gilt fiir die Gesamtzirkulation einer Schaufel 

+ b/2 
P= | V(t) ak. (15) 

% =—}/2 
_ Legt man die Wirbelbelegung auf die Geraden y=vt(yv=0, +1, +2, ...), jeweils von 
p x—= — b/2 bis x= + 6/2, so hat man ein ungestaffeltes Schaufelgitter (8,—= 7/2). Durch 
| Integration iiber die Schaufeltiefe erhalt man in Analogie zu (13) fiir die komplexe Geschwindigkeit 

. b/2 

L / 
w(z) = af y(x’) etgh 


x= — b/2 


~ (z — x’) 


dx’. (16)! 


In erster Linie interessiert die induzierte Geschwindigkeit am Ort der Schaufel selbst (Abwind- 

| verteilung), also fiir z= x. Hier wird w(z) rein imaginar, somit 
+ b/2 

1 a 


ied © i) Pe at | y (x’)ctgh ie (x— x')| dx’. (17) 


' Sehr weit vor bzw. hinter dem Gitter, wo der ctgh = + | ist, erhalt man aus (17) wegen (15) 
» wieder das gleiche Geschwindigkeitsfeld wie bei der Einzelwirbelreihe (14). Das Geschwindig- 
i keitsfeld nach (17) bleibt auch noch sinnvoll im Grenziibergang t—> 0, d.h. fiir unendliche 
* Schaufelzahl, wie hier nicht naher ausgefiithrt werden soll. Es liefert dann die ,,schaufelkon- 
} gruente*‘ Strémung der Eulerschen Stromfadentheorie. 
Um die Schaufelkontur zu berechnen, hat man in gleicher Weise wie beim Einzelfliigel dem 
. Geschwindigkeitsfeld der Wirbelverteilung eine Translationsstrémung zu iiberlagern. Im ein- 
 fachsten Fall (ungestaffeltes Gitter) ist diese senkrecht zur Gitterfront mit dem Betrag w,_*. 
| Die Schaufelkontur ergibt sich dann aus der kinematischen Bedingung, daf die Tangente an 
die Schaufelkontur parallel zur resultierendenStrémungsrichtung ist. Dies ergibt im vorliegenden 
| Fall fiir die Kontur y(x) der .,unendlich diinnen*: Schaufel 
dy wy (x) 

(ae _ Waice , a 
Bei den im Abschnitt 8 angegebenen Beispielen verwenden wir das Formelsystem (15) bis (18) 
in der Weise, dafi die Zirkulationsverteilung iiber die Schaufeltiefe y (x) vorgegeben, und nach (17) 
| die induzierte Geschwindigkeit und nach (18) die Schaufelkontur berechnet wird. 
Bei der Berechnung von Beispielen hat die Ausrechnung der induzierten Abwindverteilung 
/nach (17) numerisch zu erfolgen. Hierzu mége noch eine Bemerkung gemacht werden, die 
| gleichzeitig den Grenziibergang zum Einzelfliigel (t—> oo) erkennen 1a8t. Bei der numerischen 
| Auswertung von Gl. (17) ist die singulare Stelle des Integranden hei x’ — x stérend. Es ist zweck- 
) maBig, diese abzuspalten in der Form 

+ b/2 ; + b/2 Ss re ; ; 
w(a)=—se f per ax ay v(x) loth | (x 2) ae ae dx’. (19) 


x’ = —b/2 


| In (19) stellt das erste Glied die induzierte Geschwindigkeit des Hinzelfliigels dar, die nach der 
» Birnbaum-Glauertschen Theorie bekannt ist. Das zweite Integral stellt die induzierte Geschwin- 
| digkeit aller itbrigen Schaufeln dar. Deren Integrand ist jetzt bei x =x’ regular, da {-+-\ eine 
' Potenzreihe in a(x — x’)/t darstellt, die mit dem linearen Glied beginnt. Damit ist auch das 
-zweite Integral numerisch leicht auswerthar. Im Grenzitbergang t— co verschwindet das 
| zweite Integral iiberhaupt, so daB man dann das bekannte Ergebnis des Einzelfligels erhilt. 


1 Die Wirbelbelegung wird also statt auf der Schaufelkontur auf der «-Achse angeordnet, analog zur 
| Birnbaumschen Theorie des Einzelfliigels. Dies entspricht einer Naherung fiir schwache Schaufelwolbungen 
, und kleine Anstellwinkel (Anstellwinkel Null: stoSfreier Eintritt). Sie ist aber auch noch gut, wenn Wélbung 
und Anstellung zwar grof, das Teilungsverhaltnis t/l aber klein ist, wie an anderer Stelle gezeigt wird (vgl. 
' FuBnote 2 von S. 44). Mads 
* Wird auch eine Komponente wy, parallel zur Gitterfront iiberlagert, so bedeutet dies eine Drehung 
der Kontur um den Winkel ctg Bs = wy~/wx., wobei wir ein um den Winkel fs gestaffeltes Gitter er- 
halten. Hierbei erfolgt gleichzeitig eine Langenverzerrung der Kontur, da die Breite } erhalten bleibt. 
Obwohl hierbei unser Formelsystem nicht mehr ohne weiteres gilt, soll in einer weiteren Veréffentlichung 
_ des zweiten Verfassers (vgl. FuBnote 2 von S. 44) gezeigt werden, dafi diese Methode eine noch gute Nahe- 
rung des gestaffelten Gitters fiir schwache Staffelungen darstellt. 
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Die resultierende Geschwindigkeit an der Schaufelkontur w(x) erhalt man aus der Uber- 
lagerung der Translationsstrémung w,., wy, und der vom Wirbelfeld induzierten Geschwin- 


digkeit w,,;, w,, in der Form: 


co? 


rt xe Oa. | + Saugseite | 20 
10() = V(r a) + Cyn Fy)” | Denckseite |” a 
Dabei ist w,,—= y(x)/2, wihrend w,; durch (17) gegeben ist. Fiir die bequeme Ermittlung von 


w(x) ist in der in Fufnote 2 von S. 44 zitierten Arbeit ein Hodograph angegeben worden. 


5. Einige grundsiitzliche Betrachtungen iiber den Einflu8 der Reibung auf die Gitterstr6mung. 
Bei der Strémung um einen jeden Kérper bilden die beiden an der Hinterkante zusammen- 
flieBenden Grenzschichten hinter dem Kérper eine Nachlaufstrémung (auch Totwasser genannt). 
Bei einem Einzelkér per (imGegensatz zum Gitter) hat in einigem Abstand hinter dem Kérper 
der statische Druck wieder den ungestérten Wert p,, wie vor dem Kérper. Der vom Widerstand 
des Kérpers verursachte Energieverlust tritt hier nur noch als ein Verlust an kinetischer Energie 


in Erscheinung. Durch Anwendung des Impulssatzes! auf eine Kontrollflache K nach Abb.3 


Abb. 3. a) Zur Berechnung des Widerstandes aus dem Impulsverlust im Nachlauf; K = Kontrollflache. 
b) Zur Definition der Verdrangungsdicke 6* der Grenzschicht. Vgl. FuBnote 3, S. 49. 


erhalt man fiir den Widerstand W pro Langeneinheit des mit der Geschwindigkeit V_, ange- 
strémten zylindrischen Kérpers: 


W=e [ r(y)[V.—v(y)]ldy, (21) 


wo v(y) die Geschwindigkeitsverteilung im Nachlauf bedeutet. Dabei ist das Integral iiber einen 
Querschnitt zu nehmen, der so weit hinter dem Kérper liegt, daB dort der statische Druck den 
ungestérten Wert p,, wieder erreicht hat. Fithrt man den bei Tragfligelprofilen iiblichen Wider- 
standsbeiwert ¢C,, ein durch 


W=C12V2. (22) 
sowie eine Impulsverlustdicke # durch 
+ oo 
VERDI AGI Loken ONE (23) 
y=—»v 


so laBt sich (21) auch in der abgekiirzten Form schreiben 
= d 
Cw 2 aE S (24) 


Der einfachste Falleiner Gitterstrémung ist ein aus symmetrischen Schaufeln bestehendes 
Gitter, die senkrecht zur Gitterfront eingebaut sind, b;= 2/2, und mit einer Anstrémung senk- 
recht zur Gitterfront (Abb. 4).. In diesem Fall findet keine Ablenkung der Strémung durch 
das Gitter statt. In geniigendem Abstand hinter dem Gitter? haben sich die von den einzelnen 


1 Vel. z.B. H. Schlichting, Grenzschicht-Theorie. s. Zitat auf S. 42. 
Es sei hier vermerkt, dafi der Ausgleich der Geschwindigkeitsunterschiede hinter einem Gitter sich 
wesentlich rascher vollzieht als hinter einem Einzelkérper. Wie aus theoretischen Rechnungen bekannt ist, 
andert sich im Nachlauf (Windschatten) eines Einzelkérpers bei laminarer und turbulenter Stromung 


(H. Schlichting, Ing.-Arch. I [1930], S. 533) die Nachlaufgeschwindigkeit mit 1/|/x («= Abstand von der 


Hinterkante des Stérk6rpers), dagegen hinter einem Gitter (R. Gran Olson, Z. angew. Math. Mech. 16 [1936], 
S. 257) mit 1/x. 
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Schaufeln herrithrenden Geschwindigkeitsdellen ausgeglichen, so dafs in einem Schnitt parallel 
| zur Gitterfront eine gleichmafige Geschwindigkeit w, vorhanden ist, die wegen der Kontinuitit 
gleich ‘wy ist. Der durch den Widerstand der Schaufeln verursachte Energieverlust tritt oni 
in geniigend groBem Abstand hinter dem Gitter ausschlieBlich als Druckverlust in Eracheinung: 


Ap, = P1 — Pr. 
| Durch Anwendung des Impulssatzes auf die Kontrollfliche K ergibt sich (p, — p,)t=W, 


wo W den Widerstand einer Schaufel pro Langeneinheit bedeutet. Hieraus folgt fiir den Druck- 
| verlust infolge Reibung 2 


Alo gre (25) 


Fihrt man anstatt W den dimensionslosen Widerstandswert nach (22) ein, wobei V durch 
| w,— w, zu ersetzen ist, so hat man 


ips teeny L 
7 (26) 
Se 


| Diese Gleichung sagt aus, daB mit ab- 
»nehmendem Teilungsverhaltnis t/l der 
; auf den Staudruck der Zustrémung be- 
) zogene Reibungsdruckabfall stark zu- 
» nimmt, was ja physikalisch unmittelbar 


US) | 
| einleuchtet. fi a 
Um jetzt fiir den allgemeineren Fall ~7,* 


Pe 


hp 


‘eines Gitters mit f,==-7/2 und einer 
»beliebigen Zustrémungsrichtung das 
| Wesentliche des Reibungseinflusses zu : ee : = = - 
-erhalten, betrachten wir zunichst ein 
' gestaffeltes ,,Plattengitter‘‘nachAbb.5!. 
| Die Gitterschaufeln seien unendlich QW 
diinne ebene Platten. Die Zustrémungs- 
richtung sei zunachst parallel zu den 
| Platten gewahlt, 6,— /,. In reibungs- 
loser Fliissigkeit wiirde in diesem Fall Abb. 4. Zum Reibungseinflu8 bei der umlenkungsfreien Strémung durch 


5 . : = ein ebenes Schaufelgitter. 8, = B,;= 90°; w,=w,. Die Kontrollflache K 
| keinerlei Beeinflussung der Strémung _ befindet sich ,,sehr weit“ hinter dem Gitter, so daB dort homogene Stré- 


durch das Gitter eintreten, es gilt somit? mung herrecht. pj —p,= Ap, > 0= Druckverlust infolge Reibung. 
Boia= By Wd P2id— Pi: (27) 
| Durch den Reibungseinflu8 tritt jedoch eine Anderung der Strémungsgré8en hinter dem Gitter 
' gegeniiber diesen Werten ein, die im folgenden mit Hilfe des Impulssatzes berechnet werden soll. 
| Hierbei kénnen wir in den Schaufelkanalen die wirkliche reibungsbehaftete Strémung durch eine 
_gedachte reibungslose Strémung um verdickte Platten ersetzen, deren Dickenverteilung durch 
die Verdrangungsdicke 6* der Grenzschicht gegeben ist®. Die Potentialstrémung wird durch 
die Reibungsschicht um den Abstand 6* von der Wand abgedrangt. In den Schaufelkanalen 
wird also der Reibungseinflu8 lediglich durch die Verdrangungswirkung erfafit. Nach dieser 
| Auffassung entsteht nun an der Schaufelhinterkante fiir die aus den Schaufelkandlen austretende 
/ 


ee ff ee 


Strémung eine plétzliche Querschnittserweiterung mit einem Totwasser, welches sich strom- 
abwarts mit der umgebenden Strémung vermischt. Die hinter den Schaufeln auftretenden 
_Energieverluste kénnen als StoBverlust aus dem Impulssatz berechnet werden. 


Wir betrachten jetzt die Kontrollflache K, nach Abb. 5*. Die Geschwindigkeit w, in der 


1 Man vergleiche hierzu die analogen Betrachtungen von C. Pfleiderer (Die Kreiselpumpen. Berlin 1949) 
an einem Schaufelgitter mit endlich dicken Platten in reibungsloser Strémung. | a 
2 Mit dem Index,, 2 id‘‘ werden im folgenden die Werte hinter dem Gitterin reibungsloser(idealer) Fliissig- 
keit bezeichnet. 


3 i - 
Es gilt Vo = if (Vo — Lop; 


y=0 
wenn v(y) die Geschwindigkeitsverteilung innerhalb der Grenzschicht ist (Abb. 3b). _ : 
4 Die Austrittsseite dieser Kontrollflache liege so weit hinter dem Gitter, dafi dort eine ither die Teilung 
konstante Geschwindigkeit w, vorhanden ist. 
A 
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Austrittsebene des Gitters ist wegen der Verdrangungswirkung gréBer als w,, wahrend die iiber 
die Kanalbreite gemittelte Richtung von w, mit derjenigen von w, iibereinstimmt: 
/ 
B= By = Baia (w, > w;)- (28) 
; t 
Die Kontinuitatsgleichung liefert w= Wy, aan 


und andererseits auch 

Oi = Wyse (29) 
Weiter liefert der Impulssatz fiir 
die y-Richtung fiir die Kontrollflache 
K,, dain y-Richtung innerhalb dieser 
Kontrollflache keine Reaktionskrafte 
wirksam sind, 


y= Wy. (30) 


Aus den Beziehungen (27), (28), (29) 
und (30) ergibt sich der Geschwindig- 
keitsplan in Abb. 5. Hieraus ist er- 
sichtlich, daB 


Wy > Weiq— Wy 
und (31) 

Bo< Bria= By. 
Durch den ReibungseinfluB wird also 
hinter dem Gitter der Betrag der 
Geschwindigkeit gegeniiber dem rei- 
bungslosen Wert vergréBert, und der 
Strémungswinkel verkleinert, d. h. 
die Strémungsrichtung wird durch 
den ReibungseinfluB gegen die Um- 
fangsrichtung hin abgelenkt. Dieses 
Ergebnis gilt in gleicher Weise fiir 
ein Verzégerungs- und Beschleuni- 
gungsgitter. Die zusatzliche Ab- 
lenkung durch die Reibung 


(AB)r = Bria — Be 
ist in jedem Fall positiv. Eine quan- 
titative Formel zur Ermittlung von 
J. aus den iibrigen Daten der Stré- 
mung wird im nachsten Abschnitt 
angegeben werden. 
Der Druckverlust infolge Rei- 
: bung fiir das gestaffelte Platten- 
Abb. 5. Zur Abschétzung des Reibungseinflusses bei der Gitterstrémung. 


»,Plattengitter® mit 8, = Bs== 90°. Durch den ReibungseinfluB tritt hinter dem gitter nach Abb. 5 bei Anstrémung 
Gitter eine Ablenkung gegen die Gitterfront ein: (4$)R = Bojqg—f2> 0- parallel zu den Platten ergibt sich 


durch Anwendung des Impulssatzes 
fiir die x-Richtung auf die Kontrollflache K, zu (p,— p,) t= Wsin B, oder 


Wy 


a “Pe. ud 


__ Wsin Bs 
Bi Pers © 
Unter Einfiihrung von ¢,, nach (22) kommt 
A po eS 
= =¢w — sin Bs (32) 
ae | 


| Es sei hier vermerkt, dafB dieser Druckverlust infolge Reibung noch nicht den Energieverlust infolge 
Reibung darstellt. Hinzu kommt der Anteil der Anderung der kinetischen Energie infolge der Reibung, 
da |w,| += | wiq|. Der gesamte Energieverlust infolge Reibung wird in Abschnitt 6 berechnet werden. 
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was fiir B.= 2/2 in (26) tibergeht. Dieser zusatzliche Druckabfall infolge Reibung tritt auch 
beim Beschleunigungs- und Verzigerungsgitter auf. Infolgedessen erhiilt man beim Beschleu- 


mit Reibung erhalt hierdurch 
_gegeniiber A eine Drehung um den 


wie auch beim Einzelfliigel, in die 


_verlustlosen Strémung. Geschwindigkeitsplan 


Verzigerungsgitter : { (33) 
P2— Pi< Poa — P1- | 
In beiden Fallen ist 
Po< Praia . (33a) 


_sten Abschnitt angegeben werden. 


nigungsgitter in der wirklichen Strémung einen gréSeren Druckabfall als in der verlustlosen 


(idealen) Strémung, wahrend beim Verzigerungsgitter in der wirklichen Strémung der Druck- 
anstieg geringer ist als bei der 


Kraftteolan 


4 


Beschleunigungsgitter : 
P1— P2>Pi— Pid» | 


Eine quantitative Formel zur Er- 
mittlung von A p, = pra — py fiir 
Beschleunigungs- und Verzige- 
rungsgitter wird ebenfallsim nach- 


Der Krafteplan einer reibungs- 
losen Gitterstrémung (Abb. 6) ist 
dadurch gekennzeichnet, dafB die 
resultierende Schaufelkraft A 
senkrecht auf der Richtung von 
w,, steht, die das vektorielle Mittel 
von w, und wy;4 ist (Abb. 1). Der 
Auftrieb A kann in die Komponen- 
fen P,—U in Richtung der 
Gitterfront (Umfangskraft) und 
P,=S normal zur Gitterfront 
(Schub) zerlegt werden. Die Rei- 
bung verursacht nun auch nach 
(33a) einen zusatzlichen Schub S,, 
derimmer positiv, d. h.inRichtung 
von w,,, Wirkt. Die resultierende 
Kraft R bei der Gitterstrémung 


Gleitwinkel ¢. Sie kann wiederum, 


Komponenten A (Auftrieb) und 
W (Widerstand) zerlegt werden. 
Im Gegensatz zum Hinzelfliigel 
erfolgt beimBeschle Evugreittes Abb. 6. Geschwindigkeits- und Krafteplan der Gitterstromung mit Reibung fiir 
eine Auftriebssteigerung um As, a) Beschleunigungsgitter (Turbine): py > pojq3 bi > Baja: 
beim Verzégerungsgitter eine Auf- b) Gleichdruckgitter Wee: Pr ait By = 1 — Boja- 
. “ c) Verzégerungsgitter (Verdichter): p< pojg3 bi < Baid- 
ae omdernny ee Ase Beim In reibungsloser Strémung ist die resultierende Schaufelkraft A senkrecht zur 
Gleichdruckgitter bleibt der Auf- ,mittleren** Strémungsrichtung woo. Der ReibungseinfluB verursacht einen Druck- 
ieb d h die Reibung unverdn- verlust A p, = Pojq — Po der einen zusatzlichen Schub S,, liefert. Aus A und S,, 
a bay a oe bas . ergibt sich die resultierende Schaufelkraft R in der wirklichen (reibungsbehafteten) 
dert, der Schub wird gleich dem Strémung. Zerlegung von S,, normal und parallel zu wo ergibt 4, baw. 
Widerstand. Dieser Fall liegt beim W = Widerstand. In allen drei Fallen gilt W = S,,- sin Bx. 
Einzelfliigel (t > oo) immer vor. | 
Zwischen dem Widerstand W und dem Zusatzschub infolge Reibung S, besteht in allen drei 
Fallen (Beschleunigungs-, Verzégerungs- und Gleichdruckgitter) die Beziehung 
W= S, sin B..; (34) 


auf die wir spater zuriickkommen werden. 


1 Man beachte, da in Abb. 1 die Werte mit dem Index ..2*+ diejenigen in der reibungslosen Strémung 


bedeuten, die hier mit ,,2 id‘* bezeichnet sind. 
A* 


~ 
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6. Ermittlung der Energieverluste im Gitter aus den GrenzschichtgréBen an der Schaufel- 
hinterkante. Durch Weiterfithrung der im vorigen Abschnitt fiir das Plattengitter angestellten 
Betrachtungen soll jetzt in diesem Abschnitt fiir ein beliebiges Schaufelgitter die Berechnung der 
Energieverluste ausgefiihrt werden. Der gesamte Energieverlust laBt sich, wie schon frither 
angegeben, aufteilen in die Grenzschichtverluste und die Mischungsverluste im Nachlauf. Die 
Ermittlung der Grenzschichtverluste erfordert die Berechnung der (laminaren oder turbulenten) 
Grenzschicht von der Vorder- bis zur Hinterkante der Schaufel. Die Mischungsverluste werden 
ermittelt durch Anwendung des Impulssatzes von der Schaufelhinterkanten-Ebene bis ,,weit 
hinten*‘, wo die Nachlaufstrémung ausgeglichen ist. Bei der rechnerischen Ermittlung der Ver- 
luste werden wir jedoch nicht diese beiden Anteile getrennt berechnen, sondern nur den Gesamt- | 


Energieverlust. 


Fiir die Berechnung der Reibungsschicht zu der in Abschnitt 4 ermittelten potentialtheore- | 
tischen Geschwindigkeitsverteilung lings der Schaufelkontur haben wir bei laminarer Strémung | 
das Verfahren von K. Pohlhausen* mit der 
Verbesserung durch H. Holstein und 
T. Bohlen®, und fiir den turbulenten Fall 
das Verfahren von FE. Gruschwitz* ver- 
wendet. Die Rechnung verlauft in genau 
gleicher Weise wie beim Einzelfliigel, so 
daB es sich eriibrigt, hier naher darauf 
einzugehen. Im iibrigen sei fiir die Durch- 
fihrung der Rechnung auch auf das in 
FuBnote 1 von S. 42 zitierte Werk ver- 
wiesen. Fiir die Ermittlung der Energie- 
verluste ist es itbrigens nicht erforderlich, 
fiir die Grenzschicht an jeder Stelle die 
genaue Geschwindigkeitsverteilung zu 
kennen. Es geniigt die Ermittlung der sog. 
Grenzschicht-KenngréBen .Verdrangungs- 
dicke, Impulsverlustdicke und Wandschub- 
spannung. Insbesondere braucht man fiir 
das Folgende diese KenngréBen an der 
Schaufelhinterkante. Bedeutet v(n) die 
Geschwindigkeit in der Grenzschicht in 
Abhangigkeit vom Wandabstand 


b 


~ 
ae Sl 


/autoreite 


k——__Noch 
WAY AY \\\\ AY Kt \ WAAAY \\ 
\\ \\\ \\ \ 


Abb. 7. Nachlaufstrémung hinter einer Gitterschaufel (schematisch). a= sin B 
a) In der Hinterkantenebene: Geschwindigkeit in der Grenzschicht v(y) = 7 
auBerhalb der Grenzschicht wes Druck p’; Richtung fp; 4. 


re . . . . 
b) ,,Weit hinter dem Gitter in der wirklichen (reibungsbehafteten) und Ww, die potentialtheoretische Geschwin- 


Strémung: Geschwindigkeit w, Uronstast iiber Teilung t); Druck p,; digkeit auBerhalb der Grenzschicht an der 
la ee Schaufelhinterkante (Abb. 7), so sind die 
Verdrangungsdicke und Impulsverlustdicke an der Schaufelhinterkante 
Verdrangungsdicke : Or = i) L = We dn, (35a) 
6 Ws 
: f v(n) [ v(n) 
Impulsverlustdicke : Oy = [ — L a na dn. (35b) 
Onmtve 2 


Diese GréBen sehen wir auf Grund der Grenzschicht-Rechnung als bekannt an. 

Nach dem ZusammenflieBen der beiden Grenzschichten an der Schaufelhinterkante findet 
hinter dem Gitter ein Vermischungsvorgang statt, bis alle Geschwindigkeitsunterschiede aus- 
geglichen sind (Abb. 7). Die am Ende dieses Vorganges vorhandene Gesamtenergie der Strémung 
gegenitber der Potentialstrémung soll im folgenden berechnet werden. Der Strémungszustand 
weit hinter dem Gitter wird beschrieben durch Angabe der drei GréBen Druck, Geschwindigkeits- 
betrag und Geschwindigkeitsrichtung. Mit py, w,, 8, seien die Gré®en in der wirklichen (reibungs- 
behafteten) Strémung, und mit dem Index ,;, die GréBen in der verlustlosen Strémung weit 


1 K. Pohlhausen, Z. angew. Math. Med. 1, (1921). S. 235. 

2 H. Holstein und T. Bohlen, Ein einfaches Verfahren zur Berechnung laminarer Reibungsschichten, 
die dem Nahrungsansatz von K. Pohlhausen geniigen, Lilienthal Bericht, S. 10,5 (1940), 

3 E. Gruschwitz, Ing. Arch., 2. S. 321, (1931). 
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hinten bezeichnet, wahrend durch ps 
werden sollen. 


Fiir die Ermittlung der drei GréBen p,, w,, By aus denjenigen an der Hinterkante, die aus der 
Grenzschichtrechnung bekannt sind, stehen drei Gleichungen zur Verfiigung, namlich die Konti- 
nuitatsgleichung und die Impulsgleichungen fiir die x- und y-Richtung fiir die in Abb. 7 an- 
gegebene Kontrollflache K. 

Um diese Rechnung auszufiihren, geben wir den Ausdriicken fiir On und @yq nach (35a, b) 
noch eine etwas andere Gestalt, indem wir iiber die Gitterfrontrichtung y integrieren und ferner 
die Langen mit der Teilung t dimensionslos machen. Als Integrationsintervall wahlen wir jetzt 
die Teilung t, so da®B sich das Integral iiber die Grenzschichten an Ober- und Unterseite der 
Schaufel erstreckt. Die so erhaltenen dimensionslosen Verdrangungs- und Impulsdicken sind 


wegen dn= dy: sin fj. 


w, die Werte an der Schaufelhinterkante gekennzeichnet 


OF ob 1 OF unt ; v(y)| y ‘ 
— CP ae un = aw 5, ef 36 
= t-sin Bria J h w;, a( t ) Seed 
—~ = Onob + On es ; i een v(y) v(y) d Be (36b) 
t- sin B:id J altos t } ‘ 


Die potentialtheoretische Geschwindigkeit am Schaufelende w, ist wegen der Verdrangungs- 
wirkung der Grenzschicht gréBer als die Potentialgeschwindigkeit wy jq weit hinten, und zwar gilt 


“ W2 id 


= ye (37) 


Entsprechend ist der Druck p, an der Hinterkante kleiner als Pzia Weit hinten, und zwar gilt 


Ww 


nach der Bernoullischen Gleichung 


, Olmes 1 | 
Ps = Prid 2 W2 id la —A*)? =| . (38) 


Der Abstrémwinkel in der Hinterkanten-Ebene kann im Mittel gleich dem der Potentialstro6mung 


weit hinter dem Gitter angenommen werden, [) = [2 ia. 


1. Kontinuitat: Die Kontinuitatsgleichung liefert die Aussage, daB die x-Komponente 
der Geschwindigkeit vor und hinter dem Gitter gleich sein mu, sowohl fiir die verlustlose als 
auch fiir die wirkliche Strémung: 


ae : 
Wi, == W oid es — Woy = Weig 81M’ P2 ig = W, sin Py. 
Somit ist 
W, __ sin Paid 


Wid sin BP, 


(39) 
2. Impulssatzin y-Richtung: Dieser liefert 


1 
0 sin Psia Cos Pe ia [ v?(y)d = == 0 sin Py Cos BP, W, « 
0 
Nach Division durch gw.? und durch Hinzufiigen und Wiederabziehen von v(y)/w, und | auf 


der linken Seite ergibt sich 
: v(y) 
sin Boia * COS Boia | ae t ees 


2 2 
0 o 


wei 
t 


Hierdurch sind die GrenzschichtgréBen an der Hinterkante nach (36) eingefithrt, so das man hat 


: yo ws 
sin Bia CoS Big (1 — O A*) = sin f, cos Py En 


Durch Einfithrung von (39) und (37) und nach einigen Umformungen wird schlieBlich : 


Ope Oe Ae (40) 
tefia vd ae 


3. Impulssatzin x-Richtung: Dieser liefert 


1 
g sin® Baia [ (yd > +p, = 0 w; sin® B, + pp. 
0 
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Ganz analog zu dem Vorigen ergibt sich hieraus 
(Ore 
jee (Ee) 


P2—P2 _ 9 gin? Baia 
g Waid 

Die Druckdifferenz gegenitber dem potentialtheoretischen Wert betragt bei Einfithrung von (38) | 
und mit der Abkiirzung (40) 


SO *__ Are 
Paid Pe ae 7 — 9 (x 1) sin2 Boia : (41) | 


1 =P (y — 1) sin? Baia - 


a 


2 q 
SchlieBlich erhalten wir den gesuchten gesamten Energieverlust als den Gesamtdruckverlust | 
der Potentialstrémung zu 


2 
Wid 


/ @ ay? ey? 
81 — 82 = A8 = Pra Pg Miia (p. r3 wi). (42) 


\ 


Mit (39) und (41) ergibt sich aus (42) der Gesamtdruckverlust bezogen auf den potentialtheore- | 
tischen Austrittsstaudruck, zu 
Ag  24*—A® Bs 
ie c= 2 (x 1) sin? fs ia 
“9 Waid 
Hierbei rithren der erste und zweite Term von der Druckanderung, der letzte Term von der 


Geschwindigkeitsinderung her. Diese ergibt sich aus (39) in Verbindung mit (40) zu 
= /sin? Boia + 7? cos* Boia - (44) 


Die Anderung des Abstrémwinkels gegenitber dem potentialtheoretischen Wert ist durch (40) 
gegeben. Die Gleichungen (40), (41) und (44) geben die Lésung der oben gestellten Aufgabe, 
nimlich die Anderung von p,, , B, gegeniiber den potentialtheoretischen Werten. Dariiber 
hinaus gibt (43) noch den hauptsiachlich interessierenden gesamten Energie-Verlust. Dabei ist | 
die Ermittlung dieser GréBen zuriickgefiihrt auf die Grenzschicht-KenngréBen A* (= Verdran- 
| 
| 


(72 — 1) cos? Boia - (43) 


ee 


Wid 


gungsdichte) und 0 (= Impulsverlustdicke) an der Schaufelhinterkante. 


Im Grenzfall verschwindender Reibung, Re— co, wird 4* = O= 0, und damit aus (40) 
4=1, und somit aus (40) bis (44) py = pria, Wo = Weia y= Poias Ag = 9. 

Die GréBe 7 nach (40) ist in jedem Fall gréfer als eins, da das Grenzschichtendickenverhaltnis 
an der Hinterkante 64/0y >1 ist. Wir erhalten damit qualitativ wieder die Ergebnisse der vor- 
angegangenen groben Abschatzung fiir den EinfluB der Reibung: VergréBerung des Druckab- 
falles, VergréBerung der Abstrémgeschwindigkeit und Verkleinerung des Abstrémwinkels. 

Der Gesamtdruckverlust Ag enthalt sowohl die Reibungsverluste in der Grenzschicht als auch 
die Mischungsverluste im Nachlauf. Darin sind auch die StoBverluste infolge der plétzlichen 
Querschnittserweiterung an der Hinterkante enthalten (Carnotscher StoBverlust). Die vor- 
stehenden Gleichungen kénnen deshalb auch auf den Fall einer endlich dicken Hinterkante an- 
gewandt werden, indem man zur Verdrangungsdicke A* noch die GréBe sy/t+ sin 6, hinzuschlagt 
(sy — Dicke der Hinterkante). Der Wert der Impulsverlustdicke bleibt dadurch ungedndert. 

Gitter mit /,— 90°, ohne Umlenkung: Als Abschlu8 dieser Betrachtungen mige noch 
kurz wieder das Gitter von Abb. 4 behandelt werden, das aus symmetrisch profilierten Schaufeln 
mit B,—= 90° besteht und mit 6,= 90° angestrémt wird. Hierfiir gilt B,;a— 6, = 90°. Damit 
ergibt sich aus (42), wenn man wegen die héheren Potenzen von A* gegeniiber der ersten ver- 
nachlassigt, A4* <1 und da wegen wy iq = w,= w, noch Ag = pzia— po = Ap, gilt, die ein- 
fache Formel 

ar 263 
0 t 


2. 
7 


Fithrt man den Widerstandsbeiwert der Schaufel nach (24) ein durch 2 3/1 = ¢,,, so erhalt man 


A pv a l 
0 4 poet Ue on 
ce 


in Ubereinstimmung mit (26). 
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; 7. Definition von verschiedenen Gitterverlustheiwerten. Uber die geeignete Definitiou von 
Gitterverlustziffern und Gitterwirkungsgraden ist im Schrifttum mehrfach diskutiert worden » 23, 
Im Hinblick auf die nachfolgenden Beispiel-Rechnungen mégen hier einige Gitterverlustziffern 
und Wirkungsgrade charakterisiert und miteinander verglichen werden. Kine gewisse Schwierig- 

keit liegt darin, eine den verschiedenen Aufgaben eines Gitters angemessene Bezugsgrife fiir 
den Energieverlust im Gitter zu finden. Die drei Aufgaben eines Gitters, namlich: 1. kinetische 
Energie in Druckenergie (Pumpengitter), 2. Druckenergie in kinetische Energie (Turbinengitter) 
umzusetzen und 3. die Strémung umzulenken (Gleichdruckgitter), sind nicht immer klar von- 
einander getrennt. Zum anderen ist es notwendig, zum Vergleich der drei Gitterarten mit- 
einander einen Beiwert zu finden, der allen drei Arten gerecht wird. Hierzu scheint uns die von 
H. Marcinowski vorgeschlagene Gitterverlustziffer C, besonders geeignet, die den Gesamtdruck- 
verlust im Gitter auf den Durchsatz im Gitter bezieht. Mit G als Durchsatz in kg/s und F als 
freien Strémungsquerschnitt, gemessen in der Gitterebene, gilt fiir die Geschwindigkeitskom- 


ponente normal zur Gitterfront (y—spez. Gewicht) 
G 


Ws yF : (45) 
Bezieht man die Gitterverlustziffer auf den Staudruck von w,.. so ist 
e Ag 
Vie eg : (46) 
) Q  *o 
Sie ergibt sich aus (43) wegen wy. = wyiq sin By ja ZU 
oi Arr Genesee 
Sg he Buca (1 — A* 2 (%— 1) (4° — 1) ctg? Baia | - (47) 


Im Strémungsmaschinenbau ist es itblich, die Gitterverluste durch Einfiihrung einer Ge- 
_schwindigkeitsziffer wy zu beriicksichtigen, indem man statt des Gesamtdruckverlustes eine 
aquivalente Abminderung der Austrittsgeschwindigkeit w,;, der verlustlosen Strémung annimmt. 
Die Definitionsgleichung fiir y lautet somit 


aa ‘ 
P21 59 “2 = Paid > (pes ia)? : 


Fihrt man den Gesamtdruckverlust nach (42) ein, so folgt 


0 2 
Ag = § wha (lL—y’) - (48) 
Der Zusammenhang mit der Gitterverlustziffer €, nach (46) ist somit gegeben durch 
a oo yh yt 
6 yl — sin? Beiale > C.= are : (49) 


Eine weitere wichtige Grée ist der vor allem vom Einzelfliigel her gebrauchliche Wider- 
standsbeiwert ¢,,, der bereits in (22) eingefiihrt wurde. Es ergibt sich, wenn S, die Schubanderung 
infolge des Widerstandes ist, aus dem Krafteplan in Abb. 6 und aus (34): 

WS, sin 8, = Ap, tsin B.~ Ag t sin B.., 
wo Ap, = Pria— Po den Druckverlust infolge Reibung bedeutet, der in guter Naherung mit 
dem Gesamtdruckverlust 4g iibereinstimmt*. Damit wird ° 
pit aed =, > sin® B.,,. (50) 
Q we, -L 
Das Verhaltnis von Widerstand zu Auftrieb liefert die Profilgleitzahl ¢ = ¢,,/C,. Mit ¢, nach (10) 
ergibt sich somit aus (50): 


fw Sg Sin Bow 51 
Ca: (ctg Baia — etg By) 61) 


1 H. Marcinowski, Definition einer Gitterverlustziffer. Vortrag der Braunschweiger Gittertagung 1944 
(als Manuskript verbreitet). 
2 W. Traupel, Neue allgemeine Theorie der mehrstufigen achsialen Turbomaschinen. Ziirich 1942. 
3 C. Keller, Axialgebldse vom Standpunkt der Tragfliigeltheorie. Ziirich 1934. 
4 DaB der Widerstand W nicht exakt durch die Verlustziffer ¢, ausgedriickt werden kann, hingt damit 
zusammen, da W als die Kraft in Richtung von w,, definiert wurde. w, stellt aber das vektorielle Mittel 
aus den ,,reibungslosen‘‘ Geschwindigkeiten weit vor und hinter dem Gitter dar. 
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Einen Wirkungsgrad von Beschleunigungs- und Verzdgerungsgittern definiert A, Bets in 
der Weise, daB er den Energieverlust auf die im Gitter umgesetzte Energie bezieht. Dies gibt 
fiir Beschleunigungs- und Verzégerungsgitter verschieden definierte Wirkungsgrade, namlich 


a Ag 

Turbinengitter: 9 47 = i Sh = ; | 
Py — Po (By 19» (52) 

PO iit Pal he peru ‘| 

Pumpengitter: 7p ee aon 
Durch Einfithrung der Profilgleitzahl erhalt man hieraus 

Py eee te he ee etl Pon 53 
ea ratty es ee = wPOR Ti. © Co 


Beide Definitionen des Wirkungsgrades versagen jedoch im Falle des Gleichdruckgitters, wo die 
umgesetzte Energie gleich Null ist. f 
Einen fiir alle Falle brauchbaren Wirkungsgrad erhalt man, wenn man die vor und hinter 
dem Gitter vorhandenen Gesamtenergien ins Verhaltnis setzt, den Energieverlust also nicht 
auf die umgesetzte Energie, sondern auf die Gesamtenergie E, vor dem Gitter bezieht 
ee cle eT Ag rs 
oT Ey | Pisa ds Path 


Setzen wir das Energieniveau der Strémung fest durch 


E=0: q=0, p=P, 
messen wir also die Driicke als Differenzdriicke gegeniitber dem Druck p, hinter dem Gitter, 


dann ist wegen p, +4, —= Po + +Ag: 


ee ey 
We igre (54) 
Nach Einfiihrung der Gitterverlustziffer (46) wird 
iL 


Lan fae Cy - sin? Boia (55) 
Dieser Wirkungsgrad ermdglicht einen Vergleich von Pumpen-, Gleichdruck- und Turbinengitter. 

Grenzfalle sehr kleiner und sehr groBer Teilung. Fiir sehr kleine Teilungen t 
wachsen die auf die Teilung bezogene Verdrangungs- und Impulsverlustdicke 4* und @ nach 
Gl. (36a, b) stark an. Nach Gl. (47) wird der Gitterverlustbeiwert ¢, unendlich, wenn A*=1 
ist. Nach (36a) bedeutet A*=1, daB der Zwischenraum zwischen zwei Schaufeln an der Hinter- 
kante gerade durch die Verdrangungsdicke der beiden Schaufeln ausgefillt wird, so daB also 
nichts mehr durch das Gitter hindurchflieBen wiirde. Sollen die in Abschnitt 6 gemachten 
Annahmen noch sinnvoll bleiben, dann muf die Verdrangungsdicke klein gegen die Teilung 
bleiben, d. h. es muB A* <1 sein. Hierdurch ist fiir das Teilungsverhaltnis eine untere Grenze 
festgelegt, bis zu welcher unser Formelsystem anwendbar ist. Es muB etwa sein: 

tae OFF ob =F OF unt 


l 


(56) 


L- sin Bria 
Dem entspricht A* < 1/5. 

Fiir sehr gro®e Teilungen treffen zwar unsere Annahmen itber die Geschwindigkeits- 
verteilung in der Hinterkanten-Ebene immer weniger zu, da wir uns mehr und mehr den Ver- 
haltnissen an einer Hinzelschaufel naéhern. Trotzdem aber liefert der Grenziibergang t/l—> co 
ein durchaus sinnvolles Ergebnis. Wie hier nicht naher ausgefiithrt werden soll, erhalt man aus 
(50) fiir den Widerstandsbeiwert der Gitterschaufel mit ¢, nach (47) im Grenziitbergang t/l > co 
den Ausdruck fiir den Widerstandsbeiwert des Hinzelfliigels, wie er in (24) angegeben wurde mit 


8. Beispiele und Ergebnisse. Das vorstehend erlauterte Verfahren zur theoretischen Berech- 
nung der Gitterverlustbeiwerte soll jetzt an einigen Beispielen erprobt werden. Dabei kommt es 
uns besonders auf die charakteristische Gegeniiberstellung von Beschleunigungs- und Verzége- 
rungsgittern an, Um sinnvolle Vergleiche anstellen zu kénnen, wird fiir Beschleunigung und 
Verzégerung das gleiche Gitter zugrunde gelegt und lediglich die Strémungsrichtung umgekehrt. 
Die umgesetzte Geschwindigkeits- bzw, Druckhiéhe ist somit bei beiden Gittern gleich. Ferner 


1 Siehe FuBnote 4 von S, 43. 
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sollen Gitterteilung und Schaufelform derart variiert werden, daf fiir die Potentialstrémung 
weit vor und hinter den Gittern in allen Fallen die gleichen Strémungsverhaltnisse herrschen. 
Ks sollen also Gitter mit verschieden groBer (é) . 
Teilung potentialtheoretisch die gleiche Aufgabe i £ 

erfiillen. Um dies zu erreichen, mu mit wach- %0+-—72 
sender Teilung die Wélbung der Schaufeln zu- 
nehmen (Winkeliibertreibung). 

a) Potentialstrémung: Fiir die Berech- 96 
nung der zugehdérigen Potentialstrémung an den 
Schaufeln eignet sich das in Abschnitt 4 beschrie- 
bene Verfahren besonders gut. Dabei wird eine gz 
Zirkulationsverteilung iiber die Gitterbreite vor- 
gegeben, die sich fiir verschiedene Teilungsver- 27 0 Sy a7, 
haltnisse nur durch einen konstanten Faktor 
unterscheidet. Fiir gleiche Zu- und Abstrém- 
winkel muf dabei, wie aus (8) zu ersehen ist, die 
Zirkulation einer Schaufel /’ proportional dem 
Teilungsverhaltnis des Gitters sein. Ferner muB, 
da zunachst nur ,.unendlich diinne*: Profile nach 
dieser Methode behandelt werden sollen, stoB- 
freier Eintritt vorausgesetzt werden!. Es muB 
also die Zirkulationsverteilung iiber die Schaufel- 
tiefe y(x) am Eintritt und Austritt Null werden. 
Wir wahlen deshalb eine elliptische Zirkulations- 
verteilung iiber die Schaufeltiefe und setzen mit 


€ = 2 x/b (Abb. 8a) 


08 


OY 


Cyl (57) 
2Wxo— 

X00 Z 
Dabei ist C eine Konstante. die mit der ge- 
forderten Gesamtumlenkung zusammenhangt. 
Durch Integration von y(&) itber die Gitter- 
_breite nach (15) erhalt man die Gesamtzirkulation 


einer Schaufel in Verbindung mit (8) zu 


= ctg Pria—etgpi= C5. (58) 


Wx a0 t 


Hierdurch ist bei vorgegebenem Zu- und Ab- 
-strémwinkel die Konstante C festgelegt. Somit 
ergibt sich der Auftriebsbeiwert einer Schaufel 
nach (10) zu 


Ue Se 0 
C= Ca sim 6 i. (59) =7 0 +7 
Abb. 8. (a) Zirkulationsverteilung, (b) (c) Abwindverteilung, 
5 6 “1 (d) Kontur einer Schaufel im Gitterverband bei verschiedenen 
Schaufelform: Die Ermittlung der Schau- Teilungent/b. Die Zirkulationsverteilung ist elliptisch. 
felform erfordert zundchst die Ermittlung der Wy 7 — induzierte Abwartsgeschwindigkeit der Einzelschaufel; 
| Abwindverteilung wy (x) langs der Schaufel. ist die ,,Stromfaden-Theorie® (schaufelkongruente Strémung). 
Diese kann in bequemer Weise nach (19) aus-  Wegen C= 1 zugebériger Umlenkwinkel 46 — 75° fiir alle 
Pe Fe , . Schaufelkonturen nach Gl. (58). 
gefithrt werden, wenn man dort fiir y(x’) die 
Verteilung (57) einsetzt. Es ist zweckmafbig, hierbei den Anteil des Einzelfliigels [erster Term 
in (19)] und den EinfluB aller iibrigen Schaufeln [zweiter Term in (19)| getrennt zu berechnen. 
Der Anteil des Einzelfliigels liefert nach Birnbaum? 


t 
yg We ene eh E20 ,,..3 


Voulonigy 


Aw, = induzierter Abwind aller ibrigen Schaufeln; t/b — 0 


1 Bei unendlich diinnen Profilen mit nichtstoBfreiem Eintritt erhalt man an der Kintrittskante potential- 
theoretisch eine unendlich groBe Geschwindigkeit. Hiermit kann die Grenzschichtrechnung nicht in Gang 


gesetzt werden. 
2 Vel. FuBnote 11 von 5S. 43. 
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In Abb. 8b, c sind beide Anteile dargestellt fiir die Teilungen t/b = 0, 1/2, 3/4, 1,2. Ein Vergleich 
beider Anteile zeigt, daB z. B. fiir t/b =1 der Abwind aller iibrigen Schaufeln etwa halb so groB 
ist wie der Abwind der Kinzelschaufel. 

Die Schaufelkonturen ergeben sich aus der Abwindverteilung durch nochmalige Integration 
nach (18). Sie sind in Abb. 8d angegeben. Fiir die Einzelschaufel mit elliptischer Zirkulations- 
verteilung ist die Kontur bekanntlich ein Parabelbogen. Fiir die Schaufel im Gitter mit der 
gleichen Zirkulationsverteilung trifft dies nicht mehr zu; die Schaufel ist hier an den Enden 
etwas flacher. Die Winkelitbertreibungen an den Schaufelenden infolge der endlichen Schaufel- 
zahl sind in nebenstehender Tabelle angegeben (§* = Winkel der Schaufelkante gegen Gitter- 
front). 


Tabelle 1. Winkeliibertreibung bei ver- Fiir die weitere Rechnung miissen die Zu- und Ab- 


schiedenen Teilungsverhaltnissen und  .t+gmwinkel festgelegt werden. Hierfiir haben wir gewahlt: 
elliptischer Zirkulationsverteilung. F 3 Pree P 
| 2id co 
t/b Omec0.5 ym U0 ie20 1 
| | Gleichdruckgitter . . . . |102,5° ee on 
ctg p* | Beschleunigungsgitter, . . 60° | 35° | 4 ‘ 
ctep,  ~ : | Lew | NF eee Verzogerungsgitter . 35° | 60° | 45 


Dies entspricht beim Beschleuni- 
gungsgitter einem Geschwindigkeits- 
verhaltnis w,/w,—1,51, und einem 
auf den Durchsatz bezogenen Druck- 
gefalle nach (11) von 


Prt Ps_ I 172 
Gg sit fi esinepy 2s 
ge 


Damit ergibt sich fiir das Beschleuni- 
gungs- und das Verzégerungsgitter 
die Konstante C nach (58) zu C+ 
= 1,70, und fiir den Auftriebsbei- 
wert hat man aus (59): 


Pia 20 + = (35 - (60) 


Die nach dem in Abschnitt 4 ange- 
gebenen Verfahren ermittelte poten- 
tial-theoretische Geschwindigkeits- 


IG) 


46 


Abb. 9. Potentialtheoretische Geschwindigkeit an der Kontur der Gitter- 

schaufel fiir verschiedene Teilungsverhiltnisse t/b. Schaufelformen nach 

Abb. 13. Bei Strémung von links nach rechts, Beschleunigungsgitter mit 

By = 60°3 Bo;q = 35°; von rechts nach links, Verzégerungsgitter mit By = 35°; 
Baiq = 60°. 


U = Umschlagspunkt laminar-turbulent, angennommen im Druckminimum 
= Geschwindigkeitsmaximum. 

A = turbulenter Ablésungspunkt beim Verzigerungsgitter. 

t/b = 0 ist die ,,Stromfaden-Theorie® (schaufelkongruente Strémung). 


verteilung ist fiir Beschleunigungs- 
und Verzégerungsgitter in Abb. 9 
angegeben. Strémung von links nach 
rechts bedeutet Beschleunigungs- 
gitter, von rechts nach links Ver- 
zogerungsgitter. Die zugehdérigen 
Schaufelkonturen sind nicht vdllig 
identisch mit denen in Abb. 8, die 
fiir ein Gleichdruckgitter (Staffe- 


lungswinkel der Schaufelsehne gegen 
die Gitterfront 6,=90°) gelten. Durch die Hinzunahme der Geschwindigkeitskomponente wy, 
wird die Schaufelkontur gedreht (ctg B,=wy,/w,,), wobei eine Streckung der Sehne erfolgt, da 
die Gitterbreite b erhalten bleiben muB. Wélbung und Wolbungsriicklage erfahren hierdurch eine 
Anderung gegeniiber dem Fall des Gleichdruckgitters. Die geometrischen Parameter der unter- 
suchten Gitter sind in der Tabelle 2 zusammengestellt (x = 0: Schaufeleintrittskante). 
Fiir die Schaufelgitter mit t/b—3/4 und 1 ist die Schaufelform und die Druckverteilung auf 
der Schaufelkontur in Abb. 10 und 11 dargestellt. 


b) Grenzschichtrechnung: Der zweite Teil der Untersuchung besteht in der Ausfithrung 
der Grenzschichtrechnung fiir beide Seiten der Schaufelkontur. Fitr die laminare Grenzschicht 
wurde das Naherungsverfahren yon K. Pohlhausen und fiir die turbulente Grenzschicht das- 
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Tabelle 2. Geometrische Parameter der untersuchten Schaufelgitter vgl. Abb. 131, 
EE 


Gitter Teilung Wélbung Sabi Peete eae: aie Winkeliibertreibung 
5 wi megess ore 
= E. is 2 4irm meklage 4 bee winkel Eintritt Austritt 
Nr. des ae a some 
Gitters 
Art b | cp/ | 
ae Mt si Hh ap 7, =B,*— Bp, Y¥2=B,* — By 
i ey 
(1) Nex ell ONS 0,07 Orc OM LO 252 ole 77.52 90° B= Ne ° 
(2) | Pounce Lear iuseel oes |p sucl aut mse | 90° | — si | 4 sae 
Oman [10 | 1,0 | 0,09 | 0,50 | 1025°| 775° | 90° | —eo° | + g'o° 
(4) Beschleuni- f 0,5 | 0,35 0,07 0,43 60° 35° 45° ie 2,0° oe 18° 
(|) ung 10 | 0,71} 6,09 | 0,41 | 60° | 35° 45° = |) 414,92 *|) 214.9 
6 i g | 2.0 141-17 0,15° | 0.85 | 60> | 35° 45° | +40:7° | —24°5° 
ee rao leon | ose] ses. | foe diate | etait be 
: ; ; ; 35 60 45 61° ee 3.02 


\ 


Druckverteilung 
= iy 5 
+ 


WER 


Ape 


Abb. 10. Druckverteilung und Verlauf der Grenzschicht an der Schaufel eines 
a) Beschleunigungsgitters (Turbine). b) Verzégerungsgitters (Verdichter). 
Teilung t/b = 1; Bs = 45°; Reynolds Zahl Re = woo I/y = 105; stoffreier Eintritt. 
Beschleunigungsgitter: B, = 60°3 Bogqg = 35°; w,/w, = 1,51 
Verzégerungsgitter: f, = 35°; Byjq = 60°; w/w, = 1/1,51 = 0,66 
B= Bauchseite; R = Riickenscite. 
An der Schaufel des Beschleunigungsgitters ist die Reibungsschicht wegen des starken Druckabfalles sehr diinn; Ablésung 
tritt nicht ein. Beim Verzigerungsgitter gibt der starke Druckanstieg eine sehr dicke Reibungsschicht mit Ablésung auf 
beiden Seiten der Schaufel (A = Ablésungsstelle), und infolgedessen wesentlich gréfere Strémungsverluste. 


jenige von E. Gruschwitz verwendet2. Fiir simtliche Rechnungen wurde die Reynoldssche Zahl 


re"! og 


1 Der Winkel 6* der Schaufel am Eintritt bzw. Austritt errechnet sich fiir beliebige Anstro6mung aus 
ety BY = wy)t0ra, + cts Bo. 

2 Vel. FuBnote 3 von S.52. In beiden Fallen konnte von einer rechnerischen Vereinfachung von E. Trucken- 
brodt, Bericht 45/9 des Instituts fiir Stromungsmechanik der T. H. Braunschweig (noch nicht veroffent- 
licht), vorteilhaft Gebrauch gemacht werden. Vel. H. Schlichting, Grenzschicht-Theorie, Kap, 12 und 22. 
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gewihlt. Eine gewisse Unsicherheit bei der Durchfiihrung der Grenzschichtrechnung riihrt von 
der Wahl des Umschlagpunktes laminar-turbulent her. Wir haben die Annahme getroffen, 
da® der Umschlagspunkt im Geschwindigkeitsmaximum liegt, was nach Untersuchungen von 
H. Schlichting und A. Ulrich! fiir gréBere Re-Zahlen recht gut gutrifft. Damit ergeben sich 
die in Abb. 9 eingetragenen Umschlagspunkte. Die Grenzschichtrechnung liefert den Verlauf 
der GrenzschichtgréBen Impulsdicke ? und Verdrangungsdicke 0* lings Saug- und Druckseite 
der Schaufelkontur, ferner fiir jede Stelle langs der Wand einen Formparameter des Grenz- 
schicht-Geschwindigkeitsprofiles, welcher die Neigung zur Ablésung charakterisiert. 

Auf die Wiedergabe der Einzelheiten der 
Grenzschichtrechnung kann hier verzichtet 
werden. Fiir die anschlieBende Berechnung 
der Gitterverlustbeiwerte benétigt man aus 
der Grenzschichtrechnung lediglich die Werte 
von Impulsverlustdicke J, und Verdrangungs- 
dicke 6% an der Schaufelhinterkante. Die 
dimensionslosen Werte dieser GréBen [O und 
A* nach (36a,b)]| sindin Tabelle 3 angegeben. 

Samtliche untersuchten Beschleunigungs- 
und Gleichdruckgitter haben keine Ablésung, 
wahrend man in den Verzégerungsgittern mit 
t/b=1 und 3/4 vor allem auf der Druckseite 
Ablésung erhalt. Auf die Berechnung des Ver- 
zogerungsgitters fiir die Teilung t/b—=2 wurde 
deshalb verzichtet. In Abb. 10 ist der Grenz- 
schichtverlauf fiir das Beschleunigungs- und 
Verzégerungsgitter mit t/b—= 1 gegeniiber- 
gestellt. Die Grenze der Reibungsschicht ist 
an einer Schaufel mafstablich (nicht iiber- 
héht) eingetragen. Der starke Unterschied in 
den Reibungsschichtdicken beim Verzége- 
rungs- und Beschleunigungsgitter tritt hier- 
bei deutlich hervor. Die Ursache hierfiir ist 
im Druckverlauf an der Schaufel zu sehen. 
BeimBeschleunigungsgitter ist fast ausschlieB- 
lich Druckabfall vorhanden. Hierbei bleibt 
die Grenzschicht sehr diinn, und Ablésung 
tritt nicht auf. Beim Verzégerungsgitter muB 
die Grenzschicht fast ausschlieBlich im 
Druckanstieggebiet strémen. Dies gibt ein 
Abb. 11. Geschwindigkeitsverteilung in der Reibungsschicht und starkes Anwachsen der Grenzschichtdicke 
ee rt ae baptara ian a be und u. U. schlieBlich Ablésung. Abb. 11 soll 
[nfolge desstarkenDruckanstieges kommt die Stromung autbeden 2 ansChauliches Bild von der Geschwindig- 

Seiten der Schaufel zur Ablésung. (A = Ablésungsstelle.) keitsverteilung in der Grenzschicht an einer 

Schaufel am Verzégerungsgitter t/b — 3/4 
geben. Die verzégerte Potentialgeschwindigkeit mit starkem Druckanstieg bedingt ein sehr 
starkes Anwachsen der Grenzschichtdicke stromabwarts und fithrt auf beiden Seiten der 
Schaufel zur Ablésung und Riickstro6mung an der Wand?. 


1 H. Schlichting und A, Ulrich, Zur Berechnung des Umschlages laminar-turbulent. Jb. 1942 dtsch. 
Luftf.-Forschg., 5. I 8; vgl. auch H. Schlichting, Forschg. Ing.-Wes. 16 S. 65, (1950). 

2 In Abb, 10 und 11 ist auffallig, daB beim Verzégerungsgitter der Ablésungspunkt auf der Bauchseite 
(Druckseite) weiter vorn liegt als auf der Riickenseite (Saugseite), was mit der Erfahrung nicht recht iiberein- 
zustimmen scheint. Hierzu mége folgendes bemerkt werden: Der Systematik halber wurden Beschleuni- 
gungs- und Verzogerungsgitter mit dem gleichen Druckgefalle Ap/o/2-w?x.— 1,7 gerechnet. Fiir ein 
Beschleunigungsgitter ist dieser Betrag durchaus normal, jedoch fiir ein Verzégerungsgitter reichlich gro8. 
Die hier gewahlte elliptische Zirkulationsverteilung tiber die Tiefe gibt fiir die Vermeidung der Ablésung 
keineswegs eine besonders giinstige Schaufelform, wie man aus der Druckverteilung in Abb. 10 und 11] er- 
kennt. Der stoBfreie Kintritt liefert auf der Bauchseite ganz vorn sofort einen steilen Druckanstieg. Dies 
gibt schon vorn grofe Grenzschichtdicken, die dann bald zur Ablésung fiihren. Auf der Riickenseite 
dagegen hat man ganz vorn einen schwachen Druckabfall und danach zunachst einen sanften und spater 
einen steileren Druckanstieg. Infolgedessen bleibt hier die Grenzschichtdicke geringer als auf der Bauch- 
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Dieses charakteristische verschiedenartige Verhalten der Reibungsschicht im Beschleu- 
nigungs- und Verzégerungsgitter wird auch durch Strémungs- Pivioerapisen in Abb. 12 bestatigt. 
Auf diesen Strémungsaufnahmen nach dem Interferenzverfahren zeigen die Streifenverschie- 
bungen den Dichteunterschied 


gegeniiber der ungestérten NN 
Strémung an. Die Grenz- \S Wy i/ 
schicht wird durch ein sehr NS HI) 
plitzliches Abknicken der |, | 
Streifen angezeigt, da dort (qo | 
durch die innere Reibung gri- 

Bere Temperaturunterschiede : 

und damit gréere Dichte- . 

anderungen als in der iibrigen 


Strémung auftreten. Die Ab- HHH] 
B i ! PE 


lésungist alsstarke Verdickung * b 
der Grenzschicht zu erkennen. Abb. 12. Strémungsaufnahmen an Schaufeln im Gitterverband, nach der 
: : rR Inuerkreneemerenmethods (nach EF. Eckert). 
7 m = 
Bei der Schaufel i ‘Turbinen (a) Beschleunigungsgitter (Druckabfall): 6, = 48°; t/l= 0,524; Re L595 LO. 
gitter (Abb. iby a) tritt auf der (a) Verzégerung seitter (Druckanstieg): #B, = 20° : t/I = 0,64: Re = 1,97: 105. 
Die Pfeile geben ake Zustrémungsrichtung an. Der Rand der Grenzschicht ist durch 


Bauchseite itberhaupt keine ein plétzliches Abknicken der Streifen zu erkennen. Bei der Schaufel im Beschleuni- 
Ablésung auf, dagegen auf der gungsgitter, (a), (Turbine) auf der Bauchseite keine Ablésung, auf der Riickenseite 
es ° 5 ¥ i ecnerache Abloésung nahe der Hinterkante (gestrichelte Gerade). Bei der Schaufel des 
Riickseite in der Nahe der Hin- Verzégerung sgitte rs, (b), (Pumpe) auf der Riickenseite sehr starke Ablésung 
Serkante Se schiwache Abe unmittelbar an der Schaufelnase. 
lisung, die durch eine gestrichelte Linie angedeutet ist. Im Gegensatz dazu zeigt das Pumpen- 
5 5 I 


gitter (Abb.12b) schon in der Nahe der Schaufelnase auf der Riickenseite eine starke Ablésung, 


die ein Totwasser verursacht, das etwa den halben Schaufelkanal ausfiillt. 


Tabelle 3. Ergebnisse der Berechnung der Verlustbeiwerte fiir Gitter verschiedenen Teilungsverhiltnisses. 
|B, — Bzia| = 25°, vgl. Abb. 13 und 14. Re = 105, 


50 | 0,354| 0,02 | 0,426 |130,0| 40,7 | 0,170 | 0,934 | 2,13 | 57,4" 


ae eae |.2° | 0 
zs 3 = a3 | | \| | = = Boe 3 & An 
Belere | ee | ale | j= | 2) | al apis cE she) le 
Bol 22 | Ble | Sle rae esl ea es ee ae Een aie sd ae 
er] sis | 3B | Be | xis8 | le | Ble | Som | |e es Sie ae ie # 
Beis | le | ela] loin | + le Sle feu WM Se. hse aln Salem 
ly. ee: e{* | 3: ples. a T 4 
3 \ mS] onl” — Sn eae es i 

\ & 4 : Be Bese tess aie ee er dle eS ae 
A t/b t/l p/l c,  |10%-A* 10°. © - | yw  |100¢,| B ~| We | np, p | 1%, 
SE CS fide == DRS Sia ED 3 or oa ll 
0,028 | 0,028 | 0 0,025 | 199.2 | 139.3 0, 554 NO,080 15a" i730 16 — 0.654 
a) 10.5 | 0.5 1 0,443 | 10.6 7,41| 0,018 | 0,991 | 0,910 | 77,5° | 49 — 0,983 
(2)n; Oxton | 0,75..; 0 0.664 6.85, 4,68) 0.012 | 0,994 | 0,876 | Mee 716 ae} l 0,989 
ey | 1 | 1 0 0,887 5,36 3,55) 0,009 | 0,995 | 0,913 | 77,5° | 97 — 0,992 
(1) | 0,9 10,5 0,02 0.443 Diled 7,41| 0,019 | 0,991 | 0,973 | 76,6° A5 — 0.981 
Beschleunigungsgitter: 6, = 60°; Pyiq = 35°; w/w, = 1,51 
0,027 | 0,019| 0 0,023 | 200 138 | 1,504 |- 0,721 | 1,011 | 34,1°) 23 | 0,394 | 0,670 
(4) | 0,50 0,354 | 0 0.426 130 6.4 | 0,039 0,993 | 0,489 | 35,0° | 87 0.975 0,987 
(5) 11,00 | 0.207! 0 | 0.851 | 31] 27 | 0.017 | 0,996 | 0,425 | 350°} 200 | 0,990 | 0,995 
(6) | 2,00 1,414 | 0 1,702 I-74 TO, 007) 05999 | 073525)735,0 4:83 0,996 0,998 
(4) [0,50 | 0,354 | 0,02 | 0,426 | 77,0 | 6,4 0,055 0,991 0.690 | 33,0° 62 0.969 | 0,982 
Verzégerungsgitter: B, = 35°; By;q = 60°; w/w, = 0,66 
0,120 | 0,085 | 0 0,102 | 200 | 132 | 05685° | 0.7205 E1,91 | 58,9° Dee | 0,680 0,677 
(1) | 0,50 | 0,354; 0 | 0,426 | 76,7 | 40,7 | 0,136 | 0,948 | 1,70 | 59.1°| 25,0] 0,923 | 0,906 
(8) | 0.75 | 0,530 0 0,638 | 100.6 43,7 | 0,159 | 0,939 | 2,98 58,6° | 21,4) 0,913 | 0,894 
(7) 0.5 20,0 0,905 | 0,887 


seite, und es tritt deshalb die Ablésung dann erst weiter hinten ein. Das soeben Gesagte ergibt sich sehr 
deutlich aus dem Verlauf der Grenzschichtdicken langs der Schaufelkontur, der hier nic ht mit seteilt wurde. 

Durch leichte Abanderung der Schaufelkontur in ‘der Nahe der Eintrittskante (Verdickung “des Schaufel- 
kopfes) laBt sich leicht eine Schaufelform erhalten, die beziiglich der Ablésung giinstiger ist als die hier unter- 
suchte. Die ,,Ziichtung‘‘ solcher giinstiger Schaufelformen, wobei auch die Stofe mpfindlichkeit (Anderung 
der Verluste mit dem Zustrémwinkel) zu beriicksichtigen ist, war jedoch nicht Zweck der vorliegenden 
ersten orientierenden Rechnungen. Sie bleibt spateren Untersuchungen vorbehalten. 
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c) Die Gitterverlustbeiwerte. Der dritte und letzte Abschnitt der Untersuchung besteht 
nun in der Ermittlung der Verlustbeiwerte. Die Grenzschichtrechnung liefert die Impuls- 
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Abb. 13. Theoretisch ermittelte Verlustbeiwerte €, von ebenen Schaufelgittern 


mit verschiedenem Teilungsverhiltnis t/I. 
Ag 
Som 
5 Q 2 


ws. 
PEAS 


= Verlustbeiwert nach Gleichung (46). 


Ag = g1 — 8, = Differenz der Gesamtdriicke vor und weit hinter dem Gitter 
(Energieverlust infolge Reibung). 

Woo = achsiale Geschwindigkeitskomponente. P = Plattengitter. 

Samtliche Gitter sind so ausgelegt, da sie bei verlustloser Strémung gleiche 
Geschwindigkeitsdreiecke haben mit der Ablenkung 4f = 25°. 

Staffelungswinkel 6, = 45°; Reynolds-Zahl Re = wx I/v = 10°. 


Auftriebsbeiwert C, =-A/I 5 wr, = 1,204/L. 
Weitere Daten der Gitter nach Tabelle 2. 


verlustdicke #y und die Ver- 
drangungsdicke dj an der 
Schaufelhinterkante. Hieraus 
erhalt man die dimensionslosen 
Werte 9 und A* nach (36a, b). 
Damit kénnen dann die Ver- 
lustbeiwerte C,, Y, Cw» 7, P und 
nz nach (47), (49), (50), (53) 
und (55) ohne weiteres aus- 
gerechnet werden. In Tabelle3 
sind die Ergebnisse dieser 
Rechnung zusammengestellt. 
Fir das Teilungsverhaltnis 
t/b—=1/2 wurde neben der un- 
endlich  diimnen Schaufel 
(sy—=90) auch noch der Fall 
einer endlich dicken Hinter- 
kante mit sy/l—0,02 durch- 
gerechnet. Dies bedingt eine 
VergréBerung der Verdran- 
gungsdicke an der Hinterkante 
um sj, was neben héheren Ver- 
lusten vor allem verkleinerte 
Abstrémwinkel f, bringt. Um 
das Verhalten der Verlustbei- 
werte bei sehr kleiner Teilung 
zu erkennen, wurden fiir die 
drei Gitterarten noch je ein 
extrem kleines t/b gerechnet, 
und zwar wurde dafiir in (56) 
das Gleichheitszeichen gewahlt. 
In Abb. 13 ist der Verlust- 
beiwert C, iiber der Gittertei- 
lung t/l aufgetragen. Dabei 
mige bedacht werden, dah 
samtliche Schaufelgitter so 
gewahlt worden sind, daf sie 
gleichen Zustrémwinkel und 
gleiche Ablenkung mit 


AB =|B1 — Bo ial = 25° 


besitzen. Bei kleiner Teilung 
t/L ist deshalb die Auftriebs- 
belastung der einzelnen Schau- 
fel gering, dagegen bei groBer 
Teilung hoch. Die nach (60) 
errechneten ¢,-Werte sind in 
Abb. 13 mit angegeben. Der 
Verlauf des Gitterverlustbei- 
wertes ¢, in Abhangigkeit von 
der Teilung t/l zeigt einen 
charakteristischen Unterschied 
fiir die Beschleunigungs- und 
Verzégerungsgitter: Fiir sehr 
kleine t/l steigt ¢, fiir alle drei 
Arten von Gittern stark an, 
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was einfach daher riithrt, da bei sehr enger Schaufelstellung die Reibungsverluste wegen der 
groBen Schaufelzahl sehr gtoB werden. Fiir die Verzégerungsgitter sind die Verlustbeiwerte 
-durchweg wesentlich gréBer als fiir die Gleichdruck- und Beschleunigungsgitter, in guter Uber- 
einstimmung mit der Erfahrung. Dies riihrt daher, da der Druckanstieg im Verzigerungsgitter 
fiir die Grenzschicht wesentlich gréBere Verluste zur Folge hat als der Druckabfall im Beschleu- 
migungsgitter. Da® das Gleichdruckgitter noch etwas giinstiger als das Beschleunigungsgitter 
liegt, kommt daher, da8 der Durchsatz pro Langeneinheit des Gitters beim Gleichdruckgitter 
ein Maximum hat, wenn man Gitter mit gleichem w. vergleicht. Da ¢, ein auf den Durchsatz 
bezogener Beiwert ist, liegt das Gleichdruckgitter hierbei besonders giinstig. Beziiglich des Wir- 
kungsgrades 7, dagegen (Abb. 14) tritt dieser Vorteil des Gleichdruckgitters nicht in Er- 
scheinung. Mit wachsender Tei- 7 

lung tritt bei allen drei Arten ‘a oo 
von Gittern in Abb. 13 bei | 
gréBeren t/l ein Wiederanstieg A 
von ¢, ein, und zwar beim Ver- 
| zégerungsgitter schon bei wesent- 
lich kleinerer Teilung als beim 
Beschleunigungsgitter. Das Wie- 
deransteigen von¢, mit wachsen- 09 
dem t/l ist eine Folge davon, daB 
bei groBer Auftriebsbelastung an 
| der einzelnen Schaufel Ablésungs- 
erscheinungen auftreten, verur- 
_ sacht durch ortliche starke Druck- 
anstiege. Das giinstigste Teilungs- 


3 
T+—+2+—_}—3 + 
+ Beschleunigungsg/tter | 


is iP 
Gleichdruckgitver & 


Verzogerungsgitter c 


Deen 


verhaltnis t/1 (Minimum der ¢,- 

Kurve) ist somit fiir das Verzige- 
| rungsgitter wesentlich kleiner als 
fiir das Beschleunigungsgitter. 
- Daritber hinaus ist der optimale 
t/l-Bereich beim Verzégerungs- 49774 
_gitter wesentlich schmiiler als 
beimBeschleunigungsgitter. Beim 

Verzégerungsgitter kommt es des- 
| wegen viel mehr auf die richtige 
Wahl des Teilungsverhaltnisses 
_analsbeimBeschleunigungsgitter. | ; ae 
Zum Vergleich sind in Abb. 13 0 02 O4¥ 06 08 10 te 14 46 
' auch noch die Verlustbeiwerte Abb. 14. Theoretisch ermittelte Gitterwirkungsgrade fiir die ebenen Schaufelgitte: 
fiir ein aus ebenen Platten be- . ae 13. 
stehendes Gitter eingetragen, bei eng hee NESE 
-welchem der Zustrémwinkel /, 
gleich dem Staffelungswinkel f, ist, so daB in reibungsloser Strémung keine Ablenkung vor- 
-handen ist. Der Verlustbeiwert ¢, kann nach (50) aus dem Widerstandsbeiwert ¢,, ermittelt 
werden, der fiir die langsangestrémte ebene Platte bekannt ist. Er hangt nur von der Reynolds- 
schen Zahl ab, und zwar gilt mit Re=w, 1/y: 


nach Gleichung (55). 


- ,66 
laminar: i — ; (61a) 
0,148 
turbulent : Cy = == (61b) 


Damit wird der Verlustbeiwert des Plattengitters mit §,— 3, nach Gl. (50): 


laminar: ¢, = eer 2 5 (62a) 
esin® B ,, 
; 0,148 l 
turbulent: ¢, = (62b) 


j/Re sin® B ~ th 


Die hiernach fiir B,,—=45° und Re —10° errechneten Kurven sind in Abb. 13 mit eingetragen. 
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Die Beschleunigungsgitter und Gleichdruckgitter liegen, wie zu erwarten, zwischen demlaminaren 
und turbulenten Plattengitter, wihrend das Verzégerungsgitter gréBere Verlustbeiwerte hat 
als das turbulente Plattengitter. 

Abb. 14 gibt den Verlauf des Gitterwirkungsgrades 7,, Wie er nach (54) eingefiihrt wurde, in 
Abhangigkeit vom Teilungsverhiltnis t/l. Wie wir bereits erwahnten, ist 7, fiir den Vergleich 
der drei charakteristischen Gitterarten Beschleunigung, Gleichdruck, Verziégerung besonders 
veeignet, was sich durch diese Auftragung bestatigt. Das Gleichdruckgitter liegt mit seinem 
Wirkungsgrad zwischen Beschleunigungs- und Verzégerungsgitter. 

Kin Vergleich der hier theoretisch ermittelten Verlustzahlen, insbesondere von yr und 7, 
fiir Beschleunigungsgitter (Tabelle 3), mit den in der Praxis tiblicherweise angenommenen mag 
die theoretischen Werte als allzu giinstig erscheinen lassen. Dazu mige bedacht werden, da} 
es sich hier um die reinen Reibungsverluste des ebenen Gitters handelt, wahrend in der Praxis 
noch mancherlei andere, den Wirkungsgrad verschlechternde Einfliisse vorhanden sind. Dies 
sind einmal Verluste an den seitlichen Begrenzungen (dreidimensionale Strémung) sowie Einfliisse 
der Machzahl, die hier nicht beriicksichtigt wurden. Ferner handelt es sich hier um unendlich 
diinne Profile bei stoBfreiem Eintritt, die beziiglich der Verluste besonders giinstig liegen diirften. 
Auch wird in der Praxis vielfach die Schaufeloberflache nicht als hydraulisch glatt anzusprechen 
sein, sondern die Rauhigkeit der Schaufeloberflache einen betrachtlichen zusatzlichen Widerstand 
verursachen. 

Ferner bleibt es einer experimentellen Nachpriifung vorbehalten, ob auch in Strémungs- 
maschinen mit ihrem teilweise sehr turbulenten Strémungszustand die Grenzschicht an der 
Schaufel bis zam Druckminimum laminar strémen kann, wie sie es beim Einzeltragfliigel tut. 
Die vorherrschend beschleunigte Strémung in einem Turbinengitter begiinstigt an sich die 
Laminarhaltung sehr. Bei mehrstufigen Maschinen arbeitet aber die von den vorhergehenden 
Stufen verursachte starke Turbulenz dem entgegen. Einen beachtlichen Beitrag zu dieser Frage 
hat kiirzlich W. T. Sawyer! gegeben. Aber auch wenn die hier gemachten Annahmen beziiglich 
der Lage des Umschlagpunktes aufgegeben werden und etwa die Schaufel voll turbulent ge- 
rechnet wird, wird der charakteristische Verlauf der Verlustbeiwerte, wie ihn Abb. 13 und 14 
zeigen, erhalten bleiben. 


9. Ausblick. Wie bereits aus den wenigen hier berechneten Beispielen zu ersehen ist, liefert 
die Theorie auSerordentlich wichtige Anhaltspunkte beziiglich der Wahl eines ,.giinstigsten 
Gitters“ fiir eine vorgelegte Aufgabe, die das Gitter erfiillen soll. Die grofen Erfahrungen, die 
in der Flugtechnik bei der Auswahl giinstiger Profilformen gesammelt wurden, kénnen hierbei 
mit verwendet werden. Die hier entwickelte Theorie stellt einen ersten Anfang dar, der noch erheb- 
lich weiter ausgebaut werden muf. Um insbesondere den Hinflu8 einer Anderung des Zustrém- 
winkels (StoBempfindlichkeit) theoretisch zu untersuchen, wird es nétig sein, bequeme und zu- 
verlassige Verfahren zur Ermittlung der potentialtheoretischen Geschwindigkeit an einer endlich 
dicken Schaufel im Gitterverband auszuarbeiten, damit auch der Fall des ,,nichtstoBfreien** Ein- 
tritts behandelt werden kann. Die Beriicksichtigung des Einflusses der Machzahl auf die potential- 
theoretische Druckverteilung sowie auf die Grenzschicht bleibt ebenfalls noch offen. Systema- 
tische Rechnungen von der hier geschilderten Art, bei denen die geometrischen Parameter des 
Gitters: Profilform, Teilungsverhaltnis und Staffelungswinkel, sowie der Anstellwinkel und der 
Umlenkungswinkel zu variieren sind, liefern dann das fiir einen vorgegebenen Verwendungszweck 
»gunstigste Gitter‘*. Die Ubertragung des Verfahrens auf das Radialgitter diirfte ebenfalls 
méglich sein, wobei die Hauptschwierigkeiten immer in der zuverlassigen Ermittlung der 
Potentialgeschwindigkeit liegen. Weitere wichtige Einfliisse ergeben sich bei einem axial durch- 
stromten rotationssymmetrischen ruhenden Gitter (Leitrad) aus dem radialen Druckgefalle 


sowie daritber hinaus beim umlaufenden Rad aus der Wirkung der Zentrifugalkrafte auf die 
Grenzschicht. 


10. Zusammenfassung. Es wird ein Rechenverfahren angegeben fiir die theoretische Ermitt- 
lung der Strémungsverluste eines ebenen Schaufelgitters in inkompressibler Strémung. Diese 
Aufgabe zerfallt in die folgenden drei Teilaufgaben: 

1. Ermittlung der Druckverteilung an der Schaufel in reibungsloser Strémung; 2. Berechnung 


der Grenzschicht langs der Schaufelkontur; 3. Berechnung der Nachlaufstrémung hinter der 
Schaufel. 


' W. T. Sawyer, Experimental Investigation of a Stationary Cascade of Aerodynamic Profiles. Mitt. 
des Institutes fiir Aerodynamik der E. T, H. Ziirich, Heft Nr. 17 (1949). 
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Fir die Lésung der ersten Teilaufgabe wird ein bequemes Verfahren mitgeteilt, das eine 
Verallgemeinerung der Glauertschen Theorie der diinnen Profile darstellt. Der Arbeitsaufwand 
zur Ermittlung der Geschwindigkeitsverteilung und der Schaufelkontur betragt fiir ein 
Gitter und fiir einen Zustrémwinkel etwa drei Stunden. 

Die Kenntnis der potentialtheoretischen Geschwindigkeitsverteilung erméglicht die Durch- 
fiihrung einer Grenzschichtrechnung langs der Schaufelkontur, wobei die GrenzschichtgréBen an 
der Hinterkante der Schaufel erhalten werden. Die Berechung des Energieverlustes im Gitter 
sowie die Anderung der Austrittsstrémung gegeniiber der reibungsfreien Strémung wird auf die 
GrenzschichtgréBen an der Schaufelhinterkante zuriickgefiihrt. Der Arbeitsaufwand fiir diese 
Ermittlung einschlieBlich der Grenzschichtrechnung betraigt ebenfalls etwa drei Stunden pro 
Gitter und pro Anstellwinkel. 

An einigen Beispielen von Verzégerungs-, Beschleunigungs- und Gleichdruckgittern wird 
die Ermittlung der Geschwindigkeitsverteilung und der VerlustgréBen durchgefithrt. Dabei zeigt 
sich der entscheidende Vorteil des Beschleunigungsgitters gegeniiber dem Verzégerungsgitter in 
wesentlich geringeren Verlustbeiwerten. 


[Aus dem Institut fiir Strémungsmechanik der Technischen Hochschule Braunschweig. | 


(Eingegangen am 1. Juli 1950.) 


Anschrift der Verfasser: Prof. H. Schlichting und Dr. N. Scholz, Braunschweig, Institut fiir Stromungs- 
mechanik der Technischen Hochschule. 
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Strenge Lésung der Integralgleichung fiir die Str6mung 
durch ein senkrechtes Fliigelgitter. 


Von J. Dorr. 


1. Aufgabenstellung. Ein regelmaBiges Gitter sehr diinner Fliigel befinde sich in einem 
stromenden Medium. Das Profil aller Fliigel sei gleich, und seine Skelettlinie weiche nur sehr 
wenig von der Profilsehne ab. Die Strémungsrichtung weit vor und weit hinter dem Gitter weiche 
nur sehr wenig von der Richtung der Profilsehnen ab. Mit der in der Birnbaumschen Theorie 
itblichen Linearisierung findet man dann fiir die Wirbelbelegung J/'(x) und den Abwind w (x) 

an der Skelettlinie der Fliigel die Integralgleichung 


fe 
w(x) = | ré Re} 5}, 6tg eee 
Si 


mit h=h’+ ih’. 
Wir werden zeigen, daB man diese Integralgleichung 
auf eine Form bringen kann, die in zwei Sonder- 
fallen unmittelbar die Lésung durch den von Betz 
fiir den Einzelfliigel angegebenen liésenden Kern 
erlaubt. 


2. Transformation und Lésung einer verwandten Integralgleichung. Wir betrachten vorerst 
die Integralgleichung 


flaj=sy | sy) Gig@@ May. (1) 


Es ist 
Sf 


ree TUX Jt x See ty jo itG 
a(x—y) SMa Cof> + Sinz Cof 
Stgq oT, —== = es ae er —. (2) 
Cine eine UY 
Sut Wes om h 
Fir die Funktionen g und f ist immer folgende Zerlegung méglich: 


g(y) = 8109) +e:(y) mit | 820) =— &(—9)- (3a) 


f(x) == f(x) +e (*) mit 


Dann kann man aber auf Grund der Kigenschaften des Ketnes (2) die Integralgleichung (1) in 
folgende zwei Integralgleichungen aufspalten: 


) 
ee = a Wo A) (3b) 


+ 1 


82(Y) 
— 1g axe aK Se oe eye 
fila) = 9; Sin shee ef Sin? = — Gin? ae 4 


= 
+1 


7 6 
; l &i(y) Sin > Coj = 
Ja(*) = 55 ime eras dy. 
CU eek 


Dafitr kann man auf Grund der Beziehungen (3a) und (3b) auch schreiben 


+1 


fy (2) = 5,60) * i Olay, (4a) 
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-] _ my 
; & (y) Go} | 
fy (*) = Th Bs Tt 2 : es dy. (4b) 
na Cm 
ai 


Wir betrachten nun als erstes den Sonderfall h’’—= 0 und setzen 


Gin == en, Gin “X= cé 
h h 
ne | my 
MITES Si h ; 


Dah wegen h’’ = 0 reell ist, bildet sich durch diese Transformation die reelle Achse der y-Ebene 
auf die reelle Achse der &-Ebene ab. 


Es sei ferner 


fle())=e™, sly )=v). (6a) 
Dann ist 
+1 
th pee If Wo (E fs 
@, ()) = an V lf en cae dé, (6b) 
a4) 2 
41 
Li CS) se 
ral) =, | eS ae. (6c) 
=i 


Damit haben unsere Integralgleichungen genau die Form, fiir die Betz den lésenden Kern an- 
gegeben hat. Man findet also als Lésungen von (6b) und (6c) unter der Voraussetzung, daB h 
reell ist, 


tJ pals AE~-.- 
P2 (S) as k pee? | (ONES dy 
Vie-e@2 xf1I—& ' anfi—&) fit &7? (n— 
—1 
rd ee 
Z p(n) V1—7 
i gee eee aH 
nV 1— & 1 S 


Ersetzt man jetzt wieder 7 und € nach (5) durch x und y, so erhalt man 


ot fl ; : _ 
> [paseo Te ee tses 
( ) k Sin 9 fils) | Cin? — Sin? —— 
E2\Y ates “ paine eee dx, 
my i : | TG % wy 
wv ; J (Ls ate Bn, MOPS rae, AOS 
Co} ~ He \/ Sin’, — Cin? = h | Sin? — Sine | Sin | Sin 
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td ie ee Se 
© f(x) | Sine ™ — Sin? =* 
: : : pee PC dx 
1 at = Pir eye Ty Olt 
h | Gin? — Sin? 7 Cm on 
h oa 
Die letzte Relation kann man auch wie folgt schreiben: 
Sey, 
Sec ey eee Yt 
fule)(Gin E+ Gin eV Eel 
== : Soi —d 
8119) = ae ay oj, dx 
¥ Cit — Citta. 
i Nee ae x 
eK os UV || or oT po UH 
; Ff (%) oj —— Sin = Ce, ; 
ae ee re xe ee wey, 
hy/ Gin? ~ —Gin? = Cit Cie 
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Entsprechend folgt 


k Sin - Coj 2% 
0) =e 


Wane ee Mane gt) 
Cin? —- —Gin? —— 
x | on h Ks) h 


4${}— i 
ea pee awe = Se 1% 
| ws 9 TE 19 i12 —— — Gin? 
h |/ Gin? = — Sin? 42 Civ > Sin 
h h 
=i 
Die beiden Ausdriicke fiir g, und g, kann man jetzt aber zu einem einzigen zusammenziehen. 
Mit g = g, +, erhalt man 


Ly ss : 
k Gin + Gof =~ f(x) |/ Sine — Sine a 
8(y) — = = ie Se ES — ae =e Fm) Se 


1/ 3 “i oe te Pasa EM Bu 
n | Sin? + — Cin? “~ h | Sin® os Gir? =~ Sin h 


Damit ist die allgemeine Lésung von (1) gefunden unter der Voraussetzung, dab h reell ist. Die 
Lésung ist nicht eindeutig, da k eine beliebige Konstante ist. Diese wird aber im Falle der 
Gitterstrémung durch die AbfluBbedingung festgelegt. 

Man kann nun zeigen, daB (7) auch dann die Liésung von (1) ist, wenn h komplex ist, und 
zwar kann in h = h’ -+-ih’” die GréBe h” beliebig sein, sofern h’ >0 ist. Sofern h’= 0 ist, mub 
hingegen —1> h’’> -++1 sein. Der Beweis hierzu erfordert einigen Aufwand und soll deshalb 
nicht gebracht werden. 


3. Losung der Integralgleichung des cenkrechten Fliigelgitters. Auf Grund der Tatsache, 
da (7) die allgemeine Lésung der Integralgleichung (1) ist, 1aBt sich sofort fiir die Integralgleichung 


w(x) = [re Re a, Sta nae) : dé 


die Lésung angeben, wenn entweder h’’=0 oder h’=0 mit —1>h’’>1 ist. Der letzte Fall 
kennzeichnet ein waagerechtes Fliigelgitter, das man als Grenzfall eines sehr schiefen Gitters 
auffassen kann. Dieser Grenzfall hat keine groBe technische Bedeutung. Von groBer Wichtig- 
keit ist hingegen der Fall h’’—0, der ein senkrecht zur Strémungsrichtung stehendes Fliigel- 
gitter kennzeichnet. Man erhalt mit 


hh’ = 0 
4 +1 - : 
k Sin Gof" : w(x) |/ Sin’ Bere 7 
Box ats Si : : dx. 
| Ce RS eat ee ee SS Sin — (x — 
| Sin? = — Cin? | h| Sin? | — Sin? = Sin + é) 


—1 


Die Konstante k wird dabei durch die AbfluBbedingung 
I'(+1)=0 


festgelegt. Der Abwind w ist zu beziehen auf eine Gerade, die senkrecht steht zu der Winkel- 
halbierenden des Winkels, der aus der Strémungsrichtung weit vor und der weit hinter dem Gitter 
gebildet wird. Diese Vorschrift riithrt von der Tatsache her, daB ein regelmaBiges lineares Wirbel- 
gitter sowohl im Unendlichen vor als auch im Unendlichen hinter dem Gitter eine endliche Rich- 
tungsanderung der Strémung hervorruft und zwar sind diese Richtungsanderungen weit vor 
und weit hinter dem Gitter entgegengesetzt gleich. 


(Eingegangen am 12. Juli 1950.) 


Anschrift des Verfassers: Dr.-Ing. Johannes Dérr, Darmstadt, Technische Hochschule, Institut fiir praktische 
Mathematik, 


XIX. Band 1951. Poschl: Anwendung der Matrizenrechnung auf die Theorie der Fachwerke, 69 


Uber eine Anwendung der Matrizenrechnung auf die Theorie der Fachwerke. 


Von Th. Poschl. 


1. Gegeben sei ein statisch-bestimmtes Fachwerk mit n Gelenken oder Knoten K; und 
s== 2n— 3 Staben. Die Koordinaten der Gelenkpunkte in irgendeinem Cartesischen Achsen- 
system seien (x;, y;) (‘= 1,2, ...n), die Stablangen /;, — V(x = s,s (yi— yn)? Die in den 
Gelenken wirkenden Krafte seien, da ihre Richtungen beliebig sein kénnen, durch ihre Kom- 
ponenten (X;, Y;) (= 1, 2,...n) und die Stabkrafte in den Staben (i, k) durch S,, bezeichnet. 
Dann sind diese Stabkrafte bekanntlich durch die 2n Gleichungen 


= %,— %), 5 oy . ‘ 
Noe ye a 9 YD Suet (= 132), 7) (1) 
k lik k Lik 


bestimmt, wobei die Summen itber die Gelenkpunkte («;,, y,) zu erstrecken sind, mit denen das 
Gelenk (x;, y;) durch die Stabe verbunden ist. Da die in den Gelenken wirkenden Krifte den drei 
Gleichgewichtsbedingungen 


Sa—0, SY =—0, SF (een X;)=0 (2) 
71 t=1 2 


zu geniigen haben, so sind von den 2n Gleichungen (1) nur 2n—3 voneinander unabhangig. Die 
Stabkrafte S;, sind dann durch Aufliésung dieser Gleichungen bestimmt; wir nehmen an, daBR 
diese Auflésung eindeutig méglich, das Fachwerk also statisch-bestimmt sei. 

Wir stellen uns die Aufgabe, die Verschiebungen (u;, v;) der Gelenkpunkte K; in den Kraft- 
richtungen zu bestimmen, eine Aufgabe, die durch Anwendung der Matrizenrechnung in beson- 
ders einfacher und iibersichtlicher Weise gelést werden kann. 

Da die Gleichungen (1) linear sind, so sind auch ihre Lisungen S,, als lineare Funktionen der 
X,und Y; darstellbar. Der Einfachheit halber denken wir uns weiterhin die s-Stabe mit 1, 2,... s, 
die Stabkrafte mit (S,, S,.....S,) und die n Gelenke mit 1, 2,...n beziffert; ferner betrachten 
wir S,, S,,...5S, als die s Elemente einer einreihigen Spaltenmatrix in einem R, und die X; und 
Y; je als n Elemente in einem darin eingebetteten R,, und schreiben — was sich fiir die folgenden 
Ausfithrungen als zweckmabig erweisen wird — die Auflésung der Gleichungen (1) in der Matrizen- 


form* S=BX+1C9. (3) 
Ausfiihrlicher geschrieben haben also ©, X, ¥) die Form 
S xX, f Y, 
S, oe Y, 
el ere ees ea yee eee | 
Ss is vay 
wahrend % und © rechteckige Matrizen mit s-Zeilen und n-Spalten bedeuten: 
ina Wing oo ie, Cine Crowes Cre 
gcc | Pst» ane ++ Ban | scl emerson 
on do6 Ubi Bad: Co 0 0 OIE 
Die Gleichung (3) lautet also in ausfiihrlicherer Schreibweise : 
st Me Y, 
S, S ae i X, iss cere "") Wo : 
‘eS; Se lize ao ott Ne a CED alti \ ye , (4) 
Es ist klar, daB diese Darstellungsart — in entsprechender Ausgestaltung — auch fiir sta- 


tisch-unbestimmte Fachwerke, und auch fiir Raumfachwerke méglich ist, sowie auch fiir alle 
anderen Tragwerke, deren Theorie durch lineare Beziehungen beherrscht wird. 

2. Fiir die folgenden Anwendungen ist es von Bedeutung, daf die GroBen b;;, und C;, nicht 
durch Aufliésung der Gleichungen (1) bestimmt werden miissen, sondern auf direkte Weise durch 


1 Im folgenden werden im wesentlichen die Bezeichnungen nach R, Zurmiihl, Matrizen, Berlin-Got- 
tingen-Heidelberg 1950, verwendet. Man vergleiche auch die dort aufgefiihrten Schriften tber diesen 
Gegenstand, 


70 Péschl: Anwendung der Matrizenrechnung auf die Theorie der Fachwerke. _ Ingenieur-Archiv 


die elementaren Methoden der Fachwerkstatik erhalten werden kénnen; und zwar kann dies so- 
wohl durch Rechnung (mittels Ritterscher Schnitte), als auch durch Zeichnung (mittels der Me- 
thode der Krafteplane) unter Aufbringung geeigneter Belastungen geschehen. Es ist hierbei zu 
beachten, da® fiir wirkliche Bestimmung der Verschiebungen die GréBen 6;, und ¢;, einzeln — 
gebraucht werden. 


Man denke sich hierzu zunachst nur einen Knoten — etwa K; — mit einer Last X; = 1 be- 
lastet, wahrend alle anderen Knoten unbelastet und alle Y; = null sind. Die in den Staben wir- 
kenden Krafte ergeben dann unmittelbar die GréBen 


buss Ons es Ore 


Ebenso werden durch Belastung des Knotens K; mit 
der Last Y;—= 1, wahrend alle tibrigen Y und alle X; 
Null sind, die GréBen 


Ciget Cok nanan nce 


Durch 2 n-malige An- 
wendung dieses Vorganges 
sind die beiden Matrizen $ 


¢ und & vollstandig bestimmt. 
Sees VY — Die GréBen 6b;;, und c;;, wer- 
\ den als EHEinfluBzahlen 


(s. Zurmiihl, $.33), und die 

Matrizen 8 und € als 

4G (innere)EinfluBmatrizen 

z aa, eae ae bezeichnet. 

Abb. 1. Beispiel: rah Sale tet hoe der HinfluB matrix a bya) Me pera aieeee 

Vorgang an dem folgenden 

einfachen Beispiel: Gegeben sei ein Dreiecksfachwerk nach Abb. la mit n== 4 Gelenken und 

s==5 Staben, das durch zwei Krafte P,,(X,, Y,) und P, (X,, Y,) belastet ist. Man bestimme 
die Stabkrafte nach der angegebenen Methode. 

Zur Ermittlung der EinfluBzahlen 6;;, und ¢;, bendtigt man vier Krafteplane, fiir welche die 
Belastungen der Reihe nach durch X, = 1, X,= 1, Y,; = 1, Y, = 1 gegeben sind, wahrend alle 
anderen Lastkomponenten jeweils gleich Null sind. Die Krafteplane sind in Abb. 1b bis d ge- 
zeichnet (fiir X, = 1 reduziert sich die EinfluBmatrix auf b,,— 1). Als ,,Kinheit‘* wird eine be- 
liebige Lange gewahlt, in der dann alle anderen Komponenten gemessen werden. Man erhalt 


VAC, C =—(), 16, 0738 
0,66, I scat, Wall’, 
oe | 0-66, 0Oul J =| —0,38, +0,38 
—0,35, 0 —0719> — 0738 
+0,66, 0 Sek 5 1G 
und die Stabkrafte sind gegeben durch 
S,= 0,66 X,-+* — 0,76Y,— 0,38 Y, , 
S, = 0,66 X, + X, + 0,38 Y,+ 0,19 Y,, | 
S; = — 0,66 X,-+-* — 0,38 Y, + 0,38 Y, , (5) 
| S,= — 0,33 X,+* — 0,19Y, — 0,38Y, , 
[S50 bax ae 430 38¥e— 0,76, . 


4. Sobald die KinfluBmatrizen $ und © bestimmt sind, kann die Ermittlung der Verschie- 
bungen der einzelnen Knotenpunkte des statisch-bestimmten Fachwerks mittels des Castigliano- 
schen Satzes itber die Abgeleiteten der Formanderungsarbeit A erfolgen. Wir bestimmen zu- 
nachst die Verschiebungen der mit den eingepragten Kraften P; belasteten Knoten. 

Hierzu hat man die Ableitungen von A nach den einzelnen Lastkomponenten X,,Y;, zu bil- 
den, die unmittelbar die Verschiebungen in Richtung der betreffenden Lasten ergeben. Dazu 
ist erforderlich, das die Formanderungsarbeit A durch diese Lasten ausgedriickt wird. Dies ge- 
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ehisht durch die Formel 


x 1 aber. : 
A(S)= + >) Sir, (6) 


wobci in bekannter Weise die abkitrzenden Bezeichnungen 1; == 1;/E;F; eingefithrt sind, in denen 


jetzt J; die Stablangen, I’; die Querschnitte und E; die Elastizitatsmoduli der einzelnen Stabe be- 
deuten. 


Man hat nur die S vermége der Gleichungen (4), die wir hierzu in der Form schreiben 
S; aay (bj % X;, -|- Cik Y;) (J — ie 2 Pureeowoee S ) $ (4) 
c=] 


in die Gleichung (6) einzufiihren und erhalt die Formanderungsarbeit A, in den Belastungen 


(X, Y) ausgedriickt, 


lly l s n : 
A(XY) == >" D>) (Be Xe + op Ys)? - (6) 
j=! k=l 


Fiir die Bildung der Abgeleiteten von A nach den X,; und Y; brauchen wir die GréBen 
aS; aS; 


Ge Be Thora @) 


aus denen ebenfalls die Bedeutung der Matrizenelemente 6; und ¢;; hervorgeht. 


Die Komponenten der Verschiebungen (u;, v;) der Knotenpunkte K; nach den Achsen (x, y) 
erhalt man sodann nach dem Satz von Castigliano in der Form 


aa 8. 8S; son 
u; = Bac Sed S; Be i —— = a (bis Xi}, jel i) bi: rj 
j= | 
= STD (jx B55 7) Xn a ji] Yu) = DS) (ie Xe + Bin Ve), 
j=1 k=1 k=) 


v; = CaS eu = 5 sy (b; , X_ + Gn Yu) Gi 1; 
Ney Ee ay : : : 
=) > (dix, Ci Vj Aye Chk CFi Tj Y;) = (Vir ie OLE) 
7k = 


wobei gesetzt ist 


Ok = a eee io oe 6; 7: ; : 
k = je Oi t > Pir _ jh Ofe TF | (tea Zoran, 172) | (9) 
: s Te Eee ly iat) 
Vik = D>) Oye Git; > Ore = 2 Ga Gi Tj» | 

a= ie 


Bezeichnet man noch mit St die einreihige Spaltenmatrix mit den Elementen r;, und mit ll und % 
die einreihigen Spaltenmatrizen mit den Elementen uj, v;, also 


Y uy v1 
u v 

ae Pe a Ba| % |, 
Ys s Un . Un 


so lassen sich die Gleichungen (8) in der einfachen Form schreiben 
U=BSWVRXACH RY, | 
V= BO RKCCO RY, | 


we ; ' eg! 
wobei %’, ©’ die transponierten Matrizen zu %, © und die jeweils aus den drei ersten Teilen be- 
stehenden Matrizen die auBeren Einflu®Bmatrizen oder Verschiebungsmatrizen be- 


deuten: 


(10) 


OF | Ol = BID ILS B* = || By |= © BR, | 


s (11) 
P*=\\y_|| =BUR, A* =||6,,/= CUR. | 
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Die Elemente dieser Matrizen sind die Koeffizienten der einzelnen Glieder in dem Ausdruck fiir 
die Forminderungsarbeit A, wenn diese durch die Lasten (X;, Y;) ausgedriickt wird. 

Man erkennt ferner, daB wegen des Maxwellschen Satzes von der Gegenseitigkeit der Verschie- 
bungen die Grifen « bis 6 gewisse Symmetrieeigenschaften erfiillen miissen; und zwar lauten 
diese Beziehungen fiir den vorliegenden Fall 

Op == Ons Om = Oni fiir k=1, | 


(12) 

pix = Vki e Bri eS Vik fiir k a= L und k = 1 5 | 
Der vollstindige Ausdruck fiir die Formanderungsarbeit A, in den Lasten (X;, Yj) ausgedriickt, 
lautet daher 


1 n F 1 n nm 
A ed 9 Pa, (X; u; +Y; v;) == = = ie, (ix X; AG +2 Bix X; Ne + Oi, Y; Y;) . 


In dem oben behandelten Beispiel ergeben sich fiir die 4uBeren EinfluBzahlen a, B, y,0 die 
folgenden (abgerundeten) Werte: 
he =o 0,66 ry . 
Sse 0p, = 0,29 7, 24,07 7, = O14 7, 10/07 r, — 0,29 r, ; 
Bin = Yor = 0,25 7, + 0,125 7, — 0,25 r, +013 r, — 0,5 7 


5°? 


ba = Yo = 0,38 7, , 
Bir = Yu = —0,5 r, +0,25 1, + 0,25 r, + 0,06 r, +0,25 rz, 
Boo = Yon = 0,19 7. 


wahrend die iibrigen Elemente die Werte haben 
6&1, = 0,435 (r, re +7, + rs) +0,1257,, O5.=—= Tos 
3, = 0,58 r, 420,144 (5, Er.) 20,036 7, , 
Oo9 = 0,144 (r, +r, + 71,) + 0,036 r, + 0,585 r;. 


5. Die Bestimmung der Durchsenkung eines verschieblichen Auflagers (wie z. B. des waag- 
recht verschieblichen Knotens 3 im obigen Beispiel Abb. la) erfolgt dadurch, daB im Knoten 3 
eine fiktive Kraft X, vom Betrage | in Richtung der méglichen Verschiebung angebracht wird, 
die dann geradeso behandelt wird, wie die iibrigen Lasten. Es wird die Matrix X um ein Element 
X; erweitert und die zugehérige Spalte der KinfluBmatrix mittels eines neuen Krafteplanes be- 
stimmt. Nach Bildung der partiellen Ableitung 0A/0X, wird dann X, = 0 gesetzt, um die ge- 
suchte Verschiebung u, zu erhalten. Auflager in festen Gelenken brauchen nicht beriick- 
sichtigt zu werden. 

Die Glieder in den Gleichungen (8), die die waagrechte Verschiebung uz im verschieblichen 
Auflager 3 ergeben, sind dann von folgender Form: 


us = = = (bjx Bj3 7) Xn + Cin bjs Tj Yu) = a (%34 Xn + Bsn Yu)» (13) 
Dae = 


wobei ahnlich wie oben gesetzt ist 


s 


Mee = 2 bin Bets Bar & Ge Bis H- (14) 
r= 


Fiir das oben behandelte Beispiel ergeben sich daher fiir diese EinfluBzahlen a, und fz, die 
folgenden Werte. 
Zunichst ergibt der Krafteplan fiir die Kraft X,. der in Abb. le) gezeichnet ist, fiir die Ein- 
fluBmatrix B eine dritte Spalte mit den folgenden Werten: 
bi, = 0,66, b3— 0,66, 63, —0,66, b,,= 0,66, b;,— 0,66, 
und damit folgen die zusatzlichen Glieder der Einflu®Bmatrix: 
X3, = 0,435 (7, + 72-7, +75) — 0,25 r,. 
Sa, = 0,607... 
Bs1 = — 0,5 r, +0,25 r, + 0,25 rz — 0,125 7,1 0,25 7, , 
Bae == — 055 ry 0,125 7, 0,257, = 02517, = Oro 
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Fiir die Verschiebungen der vier Knotenpunkte 1, 2, 3, 4 des gegebenen Fachwerks erhalt mz 
sodann die folgenden Werte: ; a 
BO Ate Map Xe Pa Ya Pie Ys 
Wr its Yan 2 san at Oi Les 
Up = Og, Xy + M9 Xs + Bor Ya +fap Yo, 
% = Yor X1 + Yoo Xz +001 Yy +-Oe2 Yo» 
Ug == Oy) X + Og, X2-+ 831 Yi tf Yo, = 0, 


6. Als zweites Beispiel fiir die Anwendung des Verfahrens sei ein nicht-einfaches, statisch 
bestimmtes Fachwerk gewahlt, d.i. ein solches, das sich nicht aus Dreiecken zusammensetzt 
das also nicht aus Zweischlagen auf- und abbaubar ist. Es hat jedoch die Eigenschaft een, 
da8B sich die Stabkrafte durch Rechnung oder Zeichnung als lineare Funktionen der an den Knoten 
angreifenden Lasten darstellen lassen., Das Fachwerk be- 
steht nach Abb.2a aus n= 6 Knoten und s=2n—3=9 
Staben und ist in geometrischer und mechanischer Hin- 
sicht in bezug auf die Lotrechte symmetrisch. Die Be- 
lastung besteht aus den zwei gleichen Kraften X,, X,' an 
den beiden Knoten 1, 1’, aus der lotrechten Last Y, am 
Knoten 3 und — soferne auch die waagerechte Verschie- 
bung der Auflagerpunkte 2, 2’ ermittelt werden soll, — 
aus den beiden fiktiven Kraften X,, X,’ an diesen Punk- 


Abb, 2. Nichteinfaches Fachwerk. Ermittlung der EinfluBmatrix durch Krifteplane. 


ten 2,2’. Die Matrix Wt der elastischen Konstanten sei gleichfalls symmetrisch, fiir ihre Elemente 
r,; = 1,;/E; F; und die Groen der Lasten mégen keine besonderen Annahmen gemacht werden, so 
daB sie als Parameter in den Gleichungen stehen bleiben. 

Zur Bestimmung der inneren EinfluBmatrizen B und © fiir die Lasten X,, X, und Y; brauchen 
wir je einen Krafteplan, die in bekannter Weise (durch je zwei Rittersche Schnitte) erhalten wer- 
den kénnen und aus denen die Elemente der EinfluBmatrizen unmittelbar abgelesen werden 


kénnen. Die Krafteplane sind in Abb. 2b), c), d) dargestellt. Fiir die Krafte X,, X,, Y, erhalt 
man die folgenden EinfluBzahlen: 


nou; = 132, °b, "3,67, bay 52, Og, = — 2,45, 6,,=— 3.4, 
orev == — 1,32, b,== 2,0, bgp = 1,7, bgp = 2,3, 95. = 3,3, 
ee —— a Le Ce —= 0, Cg —— —" 42, C  — Dol, C53==— 9,0, 


wahrend alle anderen b;, und ¢;, Null sind. 
Aus diesen GréBen erhalt man sodann die Elemente a, 6, y,6 der auBeren EinfluBmatrizen 


oder der Verschiebungsmatrizen in rein schematischer Weise nach den Gleichungen 9; von diesen 
sind in dem gewahlten Beispiel nur die Gréfen 
O11> S10» Sor» %9> Biz» Pos» Yar» Vie» Yerr Yoo» Yair Ya2° d13> O33 
yon Null verschieden. Die gesuchten Verschiebungen folgen dann durch die Gleichungen 
Uy = 4, X,+Pis Y3= v1 
7 ila ast Sas / 
X, +613 Y3= % » 


y= Vu 
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Uy = Oy, X, + fos Y3= Uy 

Vy = Yor (X1 — Xj) +20 (X2 — X3) +0:3Y,;=0 (da d,;— 9), 
Us == Og, (Xy — Xj) + op (Xp — X53) + Pag Yu, = 0 (da fxg = 9), 
U3 = 31 X1 +053 Y3- 

7. Was die Anwendung dieser Darstellung auf statisch-unbestimmte Stabwerke betrifft, so 
erfolgt diese in derselben Weise wie fiir statisch-bestimmte Fille. Voraussetzung hierfiir ist, dal 
die Ermittlung der Stabkrifte durch eine der iiblichen Methoden erfolgt ist, und da die Stab- 
krifte ahnlich wie durch Gleichung (3) als lineare Funktionen der in den Knoten angreifenden 
Kriafte gefunden sind; da diese Krafte in der Ebene beliebige Richtungen haben kénnen, so miissen 
ihre Komponenten (X,, Y;) nach zwei Achsen eingefiihrt werden, so da fiir die Stabkrafte 
dieselbe Darstellung gilt wie fiir statisch-bestimmte Stabwerke in der Form 

O= Br Oe (3’) 
wobei die EinfluBmatrizen § und € nunmehr als bekannt angesehen werden. Fiir die Verschie- 
bungsmatrizen ll und & gelten dieselben Formeln (10) und (8) wie zuvor. Wenn es sich um die 
Verschiebung eines Knotens in Richtung einer dort angreifenden Last handelt, so ergibt sich die 
Verschiebung unmittelbar durch Bildung der Ableitung von A nach dieser Last, in allen anderen 
Fallen durch Einfithrung einer Hilfskraft, die nach Ausfiithrung der Ableitung Null gesetzt wird. 

8. SchluBbemerkung. In dieser Arbeit wird eine Anwendung der Matrizenrechnung auf die 
Ermittlung der Verschiebungen der Knotenpunkte eines Fachwerks gegeben, die in geschlossener 
Form durch die Einflu®matrizen des Fachwerks, durch die Matrix St der elastischen Konstanten, 
sowie durch die Lastmatrizen X, ¥) ausgedriickt werden. Das Verfahren ist auf alle Falle an- 
wendbar, die eine ahnliche Darstellung mittels linearer Funktionen zulassen. 


(Eingegangen am 7. September 1950.) 
Anschrift des Verfassers: Professor Dr. Th. Péschl, Karlsruhe, Wendtstrafbe 5. 


Nachtrag 


zu meinem Aufsatzin Bd. VIII, 5.106 des Ingenieur-Archivs 
»Lur Krmittlung der HinfluBflachen und inneren Krafte 


umfangsgelagerter Rechteckplatten“. 
Von W. Koepcke. 


Bezugnehmend auf die FuBnote auf S.132 in meinem obigen Aufsatz teilte mir Herr Prof. 
Dr.-Ing. habil. G. Worch mit, da er das darin beschriebene Verfahren zur Berechnung rechteckiger 
Platten bereits erstmalig 1939 in seinen Vorlesungen gebracht hat, es aber wegen dringen- 
derer Arbeiten erst 1943 in der August-Hertwig-Festschrift (Heft 6 der Forschungshefte aus dem 
Gebiete des Stahlbaues) unter dem Titel ,,Beitrag zur Berechnung der allseitig eingespannten 
Rechteckplatte“ veréffentlichen konnte. Somit beruht der Inhalt des Aufsatzes von Herrn Prof. 
Dr.-Ing. Worch auf eigenen Arbeiten und fubt nicht auf meiner Arbeit aus dem Jahre 1940. 


(Eingegangen am 25. September 1950.) 
Privatdozent Dr.-Ing. Werner Koepcke, Berlin-Charlottenburg, Hardenbergstr. 35. 
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LANDOLT-BORNSTEIN 


Zahlenwerte und Funktionen 
aus Physik . Chemie - Astronomie - Geophysik und Technik 


Sechste Auflage der ,,Physikalisch-Chemischen Tabellen‘® 


In Gemeinschaft mit J. Bartels, P. ten Bruggencate, K. H. Hellwege, KI. Schafer, E. Schmidt 
und unter vorbereitender Mitwirkung von J. D’Ans, G. Joos, W. A. Roth + 


herausgegeben von 


Arnold Eucken + 


In vier Binden (Erster Band in finf Teilen) 


1. Band: Atom- und Molekularphysik 


Im Mérz 1951 wird erscheinen: 


=~ 


ae : 2. Teil: Molekeln J. (Kerngeriist) 


fury 
5s 


Bearbeitet von 


P. Debye jr., E. U. Franck, F. Kerkhof, W. Maier, R. Mecke, H. Pajenkamp, H. Seidel, 
H. Stuart, E. Wicke 


Vorbereitet von Georg Joos 


f 


Herausgegeben von 


A. Eucken} und K. H. Hellwege 


Mit zahlreichen Abbildungen. Etwa 580 Seiten. 1951. 
In Moleskin gebunden DMark 168.— 


jt. y Inhaltstibersicht: 
Kerngeriist (geometrische und dynamische Higenschaften:) Atomabstinde und Strukturen, — Valenz-Energien chemischer 
Bindungen. — Trennungs-Energien chemischer Bindungen: Dissoziations-Energien zweiatomiger Molekeln. Trennungs- 


Energien mehratomiger Molekeln. — Schwingungen und Rotationen der. Molekeln. Ubersicht und Gebrauchsanweisung. 
Substanzverzeichnis. Literaturverzeichnis. Symmetrie der Molekeln und Higenschwingungen.  Eigenschwingungen (auch 
Tragheitsmomente und Kernabstande) einfacher Molekeln aus Raman- und Ultrarotspektren. Ultrarotspektren der 
einfachsten Molekeln. Ultrarotspektren weiterer ausgewahlter Molekeln, Ramanspektren weiterer ausgewahlter Molekeln. 
Mikrowellenspektren. — Die Hemmung der inneren Rotation einiger Molekeln im Gaszustand. 

Kein Physiker, Chemiker oder Ingenieur kommt gegenwartig ohne ein Nachschlagewerk aus, in dem die exakt durch 
Zahlenwerte angebbaren Ergebnisse der bisherigen Forschung iibersichtlich zusammengestellt sind. Nur wenige von 
ihnen kénnen sich mit den relativ kurzen Tabellenwerken begniigen, die eine gedrangte Auswahl aus dem Gesamtmaterial 
bringen. Die meisten Forscher und Praktiker sind auf ein ausfiihrliches Werk angewiesen, das im Prinzip das gesamte 


- in der Weltliteratur veréffentlichte Zahlenmaterial beriicksichtigt, und in dem man auch die Zitate der Originalarbeiten 


& 


findet, so daB“man sich tiber den Ursprung eines jeden Zahlenwertes unterrichten und seine Zuverlassigkeit von Fall zu 
Fall kritisch nachprifen kann. 


Die urspriinglich yon H. Landolt und R. Bornstein herausgegebenen Physikalisch-Chemischen Tabellen sind 
das erste Werk, das diese Forderung erfiillte. 


‘Die jetzt erscheinende 6. Auflege bedurfte im Hinblick auf die anschwellende Originalliteratur und mit Ricksicht auf die 
stetig steigenden Anforderungen, die an ein derartiges Werk gestellt werden, einer durchgreifenden Neugestaltung, was 
u. a. auch unmittelbar durch die Titelanderung zum Ausdruck kommt. Neu geschaffen wurde beispielsweise ein beson- 
derer Band, der die wichtigsten astronomischen und geophysikalischen Zahlenwerte und Funktionen enthilt, sowie ein 
solcher, der speziell den Bediirfnissen der praktisch tatigen Physiker, Chemiker und Ingenieure gerecht wird. 


Bisher erschien: 


I. Band, 1. Teil: Atome und Ionen. Mit 248 Abbildungen. XII, 441 Seiten. 1950. In 
Moleskin gebunden DMark 126.— 


Die weiteren Bande werden folgende Gebiete behandeln: 


I. Band, 3. Teil: Molekeln II: (Elektronenhiille), Etwa 500 Seiten. Erscheint etwa im 
Sommer 1951. 


4, Teil: Kristalle. Etwa 320 Seiten. Erscheint etwa 1951. | 
5. Teil: Atomkerne. Etwa 400 Seiten. Erscheint etwa 1951. 


II. Band: Makrophysik und Chemie. In Vorbereitung. 
III. Band: Astronomie und Geophysik. Etwa 400 Seiten. Erscheint etwa 1951. 
IV. Band: Technik. In etwa 3 Tcilen. Erster Teil erscheint etwa Ende 1951. 


Jeder Band und Bandteil ist einzeln kduflich 
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Der Kreisel 
Seine Theorie und seine Anwendungen 


Von 
Dr. R. Grammel 


o. Professor an der Technischen Hochschule Stuttgart 


Zweite, neubearbeitete Auflage 


Erster Band: 
Die Theorie des Kreisels 


Mit 137 Abbildungen, XI, 281 Seiten, 1950. DMark 30.—; Ganzleinen DMark 33.— 


Inhaltstibersicht: 


Einleitung. — Grundlagen: Grundlagen der Vektorrechnung. — Grundlagen der Mechanik. — Der Tragheitstensor. — 
Der syatiietrieche Kreisel: Der kriiftefreie symmetrische Kreisel. — Die gefiihrte Bewegung des syaumctidanen Ieee _ 
Der symmetrische Kreisel unter Zwang und Sto. — Der schwere symmetrische Kreisel. — Der Finflu8 der Rei pact 
Der wnsymmetrische Kreisel: Der kriftefreie unsymmetrische Kreisel. — Der schwere unsymmetrische Kreisel. — Be- 
sondere Probleme: Kreisel im erweiterten Sinne. — Gyroskopische Systeme: — Anhang: Darstellung der Kreiselbewegungen 
durch Thetatunktionen, — Namen- und Sachverzeichnis. ; 


Zweiter Band: 


Die Anwendungen des Kreisels 


c 
Mit 133 Abbildungen. VI, 268 Seiten. 1950. + 
DMark 30.—; Ganzleinen DMark 33.— 


Inhaltsibersicht: | ‘ 
Kreiselwirkungen bei Radsdtzen: Kollermiithlen. — Kritische Drehzahlen von Rotoren. —~ Fahrzenge. — Flugzeuge. — 
Kreiselgerate: Gerite mit KompaBkreiseln. — Der Kreiselkompa8. — Kiinstliche Horizonte mit Pendelkreisen. — Wende- 
und Lagekreisel. - Sonstige Kreiselgerite. — Unmittelbare Stabilisatoren: Richtkreisel. — Stiitzkreisel. — Dimpfkreisel. 
— Namen- und Sachverzeichnis. 
* 


Praktische Funktionenlehre 
Von 
Professor Dr.-Ing. Friedrich Télke 


Erster Band: 


Elementare und elementare transzendente Funktionen 


Zweite, stark erweiterte Auflage 


Mit 178 Abbildungen, 50 durchgerechneten Beispielen und einer Ausschlagtafel. 
XI, 440 Seiten. 1950. Ganzleinen DMark 39,— 


Inhaltsiihersicht: 


I, Abschnitt: Definierende Differential- und Integralgleichungen, Fundamentaleigenschaften und gegenseitige Be- 
ziehungen der elementaren und elementaren transzendenten Funktionen. — II. Abschnitt: Durch elementare und elemen- 
tare transzendente Funktionen ausdriickbare Integrale.—JII. Abschnitt: Funktionentafeln der elementaren Transzen- 
denten. — IV. Abschnitt: Anwendungen im Bereich der partiellen Differential- und Integralgleichungen. — V. Ab- 
sehnitt: Summierung von Reihenentwicklungen. — VI. Abschnitt: Tafel der zonalen Kugelfunktion ‘sowie ihrer Ab- 
teilungen und Integrale. 


' Auf zahlreichen Gebieten der Technik zeigt sich ein immer fiihlbarer werdendes Bedirfnis nach einer breit angelegten 
Darstellung der Praktischen Funktionenlehre. Bei der Inangriffnahme neuer Probleme ist man heute gezwungen, sich 
mit der Lésung von Integralen, Differential- und Integralgleichungen zu befassen, die langst technisches Allgemeingut 
sein kénnten. Um diese Verhaltnisse zu bessern, hat sich der Verfasser entschlossen, ein den heutigen technischen Bediirf- . 
nissen angepafites mehrbindiges Lehr- und Nachschlagebuch der Praktischen. Funktionenlehre zu schaffen. Die, vor- 
liegende, zweite Auflage des I. Bandes bringt cine Erweiterung um rund 180 Seiten. Uber den in der ersten Auflage 
behandelten Stoff hinaus enthiilt sie véllig neu die Abschnitte IV, V und VI (siche oben Inhaltsiibersicht) und wird daber 
auch viele Benutzer der ersten Auflage interessieren. on 
———— 
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